《MATLAB中解三对角线性方程组追赶法程序.doc》由会员分享,可在线阅读,更多相关《MATLAB中解三对角线性方程组追赶法程序.doc(1页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
MATLAB中解三对角线性方程组追赶法程序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