we already have rref so it shoud be easy to implement to solve matrix as equations system (RHS verctor is parameter)