<<VilGebra`

LinearSolveCompleteSP[(({{-4, -1, -1, 2}, {4, -1, 0, 0}, {-3, 0, -1, 3}, {-3, -1, 1, -3}})),(({{1}, {-1}, {-1}, {4}}))]