zhuigan.m

来自「在matlab环境下利用追赶法解线性方程组的算法实现」· M 代码 · 共 12 行

function x=zhuigan(a,b,c,f,m) 
for i=2:m
    a(i)=a(i)/b(i-1);
    b(i)=b(i)-c(i-1)*a(i);
    f(i)=f(i)-a(i)*f(i-1);
end
f(m)=f(m)/b(m);
for i=m-1:-1:1
   f(i)=(f(i)-c(i)*f(i+1))/b(i);
   x=f;
end