function x=cggm(A,b,x0,tol)
%x=cggm(A,b,x0,tol)
%線性方程組共軛梯度法求解程序。
%時(shí)間：二零零七年五月二十八日
%地點(diǎn)：南華大學(xué)
%作者：王禮廣
%A為系數(shù)矩陣，b為右向量,x0為迭代初值,tol為誤差容許

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