数值分析上机试验.pdfVIP

数值分析上机试验

一、解线性方程组干脆法(教材49页14题)

追逐法程序如下：

functionx=follovmp(A,b)

n=rank(A);

For(i=1:n)

if(A(i,i)==0)

disp(Error:对角有元素为0,);

return;

end

end;

d=ones(n,1);

a=ones(n1,1);

c=ones(n1);

for(i=l：n1)

a(i,1)=A(i+l,i);

c(il)=A(i,i+l);

z

d(il)=A(ii);

zz

end

d(n,1)=A(n,n);

for(i=2:n)

d(i,l)=d(i,l)(a(il,l)/d(il,l))*c(il,l)；

b(il)=b(i,l)(a(il,l)/d(ill))*b(il,l)；

zz