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 | ⎤ ⎦ |