zuigan.m

来自「关于一的微分方程数值解,一维波动方程初边值问题Crank-Nicholson格式」· M 代码 · 共 19 行

function [x]=zuigan(a,b,c,f)
%三对角方程的追赶法；
%a,b,c为A的三个对角，

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