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;