Question

用matalb编写程序，用追赶法求解三对角线性方程组：

A为：对角元为 十个 4 的十阶方阵 两侧对角元 全部是 -1
等式右端十个数为 b=(7 5 -13 2 6 -12 14 -4 5 -5)
求解 Ax=b

说明尽量详细些，谢谢。如果要分，说声，追加有的是！

Answer 1

function x=Trid(a,b,c,d)
% 追赶法求解三对角的线性方程组 Ax=d
% b为主对角线元素，a，c分别为次对角线元素，d为右端项
% A=[ b1 c1
% a2 b2 c2
% ......
% a_(n-1) b_(n-1) c_(n-1)
% a_(n) b_(n) ]
% b=[b1...b_(n)]
% a=[0 a2...a_(n)]
% c=[c1...c_(n-1)]

n=length(b);

u(1)=b(1);
for i=2:n
l(i)=a(i)/u(i-1);
u(i)=b(i)-l(i)*c(i-1);
end

y(1)=d(1);
for i=2:n
y(i)=d(i)-l(i)*y(i-1);
end

x(n)=y(n)/u(n);
for i=n-1:-1:1
x(i)=(y(i)-c(i)*x(i+1))/u(i);
end