给我一个共轭梯度法法解矛盾线性方程组的matlab代码

以下是使用共轭梯度法解矛盾线性方程组的matlab代码：

function [x,flag,relres] = mycg(A,b,x0,maxIter,tol)
% A: 系数矩阵
% b: 常数向量
% x0: 初始解向量
% maxIter: 最大迭代次数
% tol: 计算精度

n = size(A,1);
x = x0;
r0 = b - A*x;
p = r0;
for i = 1:maxIter
Ap = A*p;
alpha = (r0'*r0) / (p'*Ap);
x = x + alpha*p;
r1 = r0 - alpha*Ap;
relres = norm(r1) / norm(b);
if relres < tol
flag = 0;
return
end
beta = (r1'*r1) / (r0'*r0);
p = r1 + beta*p;
r0 = r1;
end
flag = 1;
end

这个matlab代码使用共轭梯度法来解矛盾线性方程组，其中A为系数矩阵，b为常数向量，x0为初始解向量，maxIter为最大迭代次数，tol为计算精度。