fdm.m

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

function x=fdm(A,b,x0,tol)
%x=fdm(A,b,x0,tol)
%线性方程组求解程序。
%时间：二零零七年五月二十八日
%地点：南华大学
%作者：王礼广
%A为系数矩阵，b为右向量,x0为迭代初值,tol为误差容许
%


rk=b-A*x0;
xk=x0;
d=norm(rk);
while d>tol
    pk=A*rk;
    ak=(rk'*rk)/(pk'*rk);
    xk1=xk+ak*rk;
    xk=xk1;
    rk=b-A*xk;
    d=norm(rk);
end
x=xk;
return;