ex3.m

来自「重点内容是计算中的误差、函数方程求根、插值呆逼近、数值积分和微分、线性代数方程组」· M 代码 · 共 26 行

n=10;
a=zeros(n-1,1);b=zeros(n,1);c=zeros(n-1,1);d=zeros(n,1);
a=linspace(1,1,n-1);b=linspace(-4,-4,n);c=linspace(1,1,n-1);
for i=1:n;
d(i)=-2;d(1)=-3;d(n)=-3;
end
A=diag(a,-1)+diag(b)+diag(c,1);
for w=0.9:-0.1:0.1
x=zeros(n,1);
m=0;det=0;ep=1.0000e-006;  
while(1)
    for i=1:n;
        s(i)=x(i);
        t=0;
        for j=1:n;
            t=t+A(i,j)*x(j);
        end
        x(i)=(d(i)-t)*w/A(i,i)+s(i);
        det(i)=abs(s(i)-x(i));
    end
    if norm(det,inf)<ep break;end
    m=m+1;
end
k(m)=m;
end
k