请问python怎么求解方程组？numpy好像只能求方阵的情况

刚学习python，用来做一些科研方面的计算，目前碰到到一个问题，要解一个方程组，A*B=C，其中A和C是已知的一些数，A是4行9216列，C是4行1列，B是待求的举证，为9216行1列。按照线性代数来说，我求出A的逆矩阵，然后左乘到方程上就能得到B矩阵，但是目前了解到的是numpy只能求行列相等的方阵的逆，这种行列不等的矩阵用pinv求逆矩阵然后乘上去，结果验证过不对。请问对于非方阵的情况下，怎么用python解方程组？这种问题在平时科研中经常会碰到。

确定有解吗？ 我记得pinv是能给出假解的