Basic Jacobi solver works
Some optimisations
Separate functions
Three types of solvers
WS4