The jacobi_linsolve command finds the solution of a linear system of equations using the Jacobi iteration method.
Examples.
A:=[[100,2],[2,100]]; |
jacobi_linsolve(A,[0,1],1e-12); |
⎡ ⎣ | −0.000200080032,0.0100040016006 | ⎤ ⎦ |
⎡ ⎣ | −0.000200080032013,0.0100040016006 | ⎤ ⎦ |