zhuigan.m

来自「用matlab编写的追赶法程序」· M 代码 · 共 16 行

function[x]=zhuigan(a,b,c,d,n)
l=zeros(n,1);u=zeros(n,1);y=zeros(n,1);
u(1)=b(1);
for i=2:n
    l(i)=a(i-1)/u(i-1);
    u(i)=b(i)-l(i)*c(i-1);
end
y(1)=d(1);
for i=2:n
    y(i)=d(i)-l(i)*y(i-1);
end
x(n)=y(n)/u(n);
for i=n-1:-1:1
    x(i)=(y(i)-c(i)*x(i+1))/u(i);
end