For linear systems, see Matrix.solv().
Prints the form of the matrix defined by eqn statements. Each nonzero element is printed as an "*".
Does one iteration of the non-linear system defined by eqn statements. Returns the linear norm of the difference between left and right hand sides plus the change in the dependent variables.