Last 2 dimensions of the array must be square

这个报错是由于咱们在求解行列式的值的时候使用了：

np.linalg.det(D)

可是D必须是方阵才能够进行运算，不是方阵则会报错，咱们把以前的行列式更改成方阵就不会再报错了，固然这也是numpy自身计算数值的一个缺陷。