Matrix splitting (source code)

= Matrix splitting
{wiki}

A splitting $A=M-N$ with invertible $M$ produces
$$
x^{(k+1)}=M^{-1}Nx^{(k)}+M^{-1}b.
$$