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

>>A=[-3 2 6;10 -7 0;5 -1 5]; >> b=[4 7 6]'; >> x0=[1 0 0]'; >> M=inv(diag(diag(A)))*A; >> [x,n]= preconjgrad (A,b,x0,M,eps) 未定义与 'double' 类型的输入参数相对应的函数 'preconjgrad' 您好，请帮看一下是哪里输错了吗？

有点看不懂