我正在尝试解决具有多个解决方案的线性系统。这是我的代码的一部分：
from scipy.sparse.linalg import lsqr
solution = lsqr(M, b)[0]
现在，如果矩阵M为：
[1 1 1 1]
[0 0 0 0]
[0 0 0 0]
[0 0 0 0]
[0 0 0 0]
[0 0 0 0]
b是这个
[1 0 0 0 0 0]
此代码给出的解决方案是：
[1/4 1/4 1/4 1/4]
但是，出于我的目的，我想获得一个尽可能多的零向量，因此在这种情况下，应该是这样：
[1 0 0 0]
有没有办法做到这一点？我也可以使用scipy.sparse.linalg以外的软件包。谢谢！