Last 2 dimensions of the array must be square

這個報錯是由於咱們在求解行列式的值的時候使用了：

np.linalg.det(D)

可是D必須是方陣才能夠進行運算，不是方陣則會報錯，咱們把以前的行列式更改成方陣就不會再報錯了，固然這也是numpy自身計算數值的一個缺陷。