cgm(A,b)

Solves Ax=b using the conjugate gradient method.

Solves Ax=b using the conjugate gradient method. The matrix A should be symmetric and positive definite.
Usage (examples provided at end of m-file):
x=cgm(A,b);

If the matrix is sparse, it is worth trying:
x=cgm(sparse(A),b);