Added new Gauss-Seidel routine - not working yet.

This is a naive implementation. More sophisticated approach to follow.