预处理共轭梯度法求线性方程组Ax=b的解

function [xn]= preconjgrad （Abx0Meps）
if nargin == 4
    eps = 1.0e-6;
end
r1 = b-A*x0;
iM = inv（M）;
z1 = iM*r1; 
p  = z1;
n  = 0;
tol= 1;
while tol>=eps
    alpha  = dot（r1z1）/dot（pA*p）;     
    x  = x0 +