febs.m

matlab编程的三个线性方程组求解子程序：cggm.m,fdm.m,febs.m.其中cggm.m为共轭梯度法求解子程序,fdm.m为最速下降法求解自程序,febs.m为三对角追赶法求解子程序.

function xx=febs(a,b,c,f)
%x=febs(a,b,c,f)
% 对角占优的三对角线性方程组追赶法求解程序。
%时间：二零零七年五月二十八日
%地点：南华大学
%作者：王礼广
%a,b,c为系数，f为右向量

n=size(f,1);
for k=1:n
    if b(k)==0
         error('There is a zeros B element %d in the data!',k);
    end
end
for k=1:n-1
    if c(k)==0
         error('There is a zeros C element %d in the data!',k);
    end
end
for k=1:n-1
    if a(k)==0
         error('There is a zeros A element %d in the data!',k);
    end
end
if abs(b(1))<=abs(c(1)) || abs(b(n))<=abs(a(n-1))
    error('There is no mask element is domainated!');
end
for k=2:n-1
    if abs(b(k))<abs(a(k-1))+abs(c(k))
        error('There is no mask element is domainated!');
    end
end

s=zeros(n-1);
s(1)=c(1)/b(1);
for k=2:n-1
    s(k)=c(k)/(b(k)-a(k-1)*s(k-1));
end
y(1)=f(1)/b(1);
for k=2:n
    y(k)=(f(k)-a(k-1)*y(k-1))/(b(k)-a(k-1)*s(k-1));
end
x(n)=y(n);
for k=n-1:-1:1
    x(k)=y(k)-s(k)*x(k+1);
end
xx=x';
return;