数值分析 追赶法 matlab编程
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
clear;
A=[0 -1 -1 -1 -1];
B=[2 2 2 2 2];
C=[-1 -1 -1 -1 0];
f=[1 0 0 0 0];
n=length(B);
C(1)=C(1)/B(1);
for i=2:n-1
C(i)=C(i)/(B(i)-A(i)*C(i-1));
end
f(1)=f(1)/B(1);
for i=2:n
f(i)=(f(i)-A(i)*f(i-1))/(B(i)-A(i)*C(i-1)); end
x(n)=f(n);
for i=n-1:-1:1
x(i)=f(i)-C(i)*x(i+1);
end
x