MATLAB中解三对角线性方程组追赶法程序
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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:n
l(i)=b(i)-a(i)*u(i-1);
u(i)=c(i)/l(i);
end
y(1)=d(1)/l(1);
for i=2:n
y(i)=(d(i)-y(i-1)*a(i))/l(i);
end
x(n)=y(n);
for i=n-1:-1:1
x(i)=y(i)-u(i)*x(i+1)
end
x,
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