The conjugate_gradient command takes two mandatory arguments
and two optional arguments. The mandatory arguments are an n× n positive
definite symmetric matrix A and a vector y of length n. The
optional arguments are a vector x0 of length n and a positive
number є.
conjugate_gradient uses the conjugate gradient algorithm to
return the solution to Ax = y to within є or
epsilon. The vector x0 is an optional initial approximation.
Input:
Output:
Input:
Output:
Input:
Output: