Implement solutions for the system `Ax = b`

* implement back substitution
* implement solutions for `Ax = b` by gaussian elimination
* implement solutions for `Ax = b` by gauss jordan elimination
* Add tests
* Correct the bug with fraction free gauss jordan elimination: Now the diagoal entries
  are equal to determinant of the matrix
Refactor code