?? jacobi.txt
字號:
% A為方程組的增廣矩陣
% e為誤差容量
function x=myJacobi(A,ep)
if nargin<2
ep=1.0e-5;
end
[n,m]=size(A);
x=zeros(n,1);
y=zeros(n,1);
k=0;
while 1
for i=1:1:n
s=0.0;
for j=1:1:n
if j~=i
s=s+A(i,j)*x(j);
end
y(i)=(A(i,n+1)-s)/A(i,i);
end
end
xmx=0.0;
for i=1:1:n
xmx=max(xmx,abs(x(i)-y(i)));
end
if xmx<=ep
for i=1:1:n
x(i)=y(i);
end
return;
end
for i=1:1:n
x(i)=y(i);
y(i)=0.0;
end
k=k+1;
if k>30
error('failure');
return;
end
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -