Added reduced row echelon form, rref
Added function solve_linear for solving matrix of linear equation system