Initially an Linear Equation Solver by Gauss Seidel method, but it can be used as a library for main matrix operations in Fortran