MATLAB中解三对角线性方程组追赶法程序(共1页).docx

精选优质文档-倾情为你奉上MATLAB中解三对角线性方程组追赶法程序M文件函数：function x=followup(a,b,c,d)n=length(b);u(1)=c(1)/b(1);l(1)=b(1);for i=2:nl(i)=b(i)-a(i)*u(i-1);u(i)=c(i)/l(i);endy(1)=d(1)/l(1);for i=2:ny(i)=(d(i)-y(i-1)*a(i)/l(i);endx(n)=y(n);for i=n-1:-1:1x(i)=y(i)-u(i)*x(i+1)endx,end例子： a=0 -1 -1 ;b=4 4 4;c=-1 -1 0;d=2 4 10 ; x=followup(a,b,c,d)x = 1.0000 2.0000 3.0000专心-专注-专业