<<VilGebra`

LinearSolveCompleteSP[{{-5,2,1},{2,-1,3}},{{0},{0}}]