````python
import numpy as np
代码4-3
方法一:使用det函数求行列式
arr = np.array([[2, 1, -5,1],[1, -3, 0,-6],[0,2, -1,2],[1,4,-7,6]]) # 创建分母的三维数组
arr1 = np.array([[8, 1, -5,1],[9, -3, 0,-6],[-5, 2, -1,2],[0,4,-7,6]]) # 创建分子的三维数组
arr2=np.array([[2,8,-5,1],[1,9,0,-6],[0,-5,-1,2],[1,0,-7,6]])
arr3 = np.array([[2,1, 8,1],[1, -3, 9,-6],[0, 2, -5,2],[1,4,0,6]]) # 创建分子的三维数组
arr4 = np.array([[2, 1, -5,8],[1, -3, 0,9],[0,2, -1,-5 ],[1,4,-7,0]]) # 创建分子的三维数组
求解行列式
D = np.linalg.det(arr)
D1 = np.linalg.det(arr1)
D2 = np.linalg.det(arr2)
D3 = np.linalg.det(arr3)
D4=np.linalg.det(arr4)
print(‘方程组的解x1为:’, D1 / D)
print(‘方程组的解x2为:’, D2 / D)
print(‘方程组的解x3为:’, D3 / D)
print(‘方程组的解x4为:’, D4 / D)
方法二:使用solve函数解线性方程组
D = np.array([[2, 1, -5,1],[1, -3, 0,-6],[0, 2, -1,2],[1,4,-7,6]]) # 创建系数行列式
arr = np.array([8, 9, -5,0])
x = np.linalg.solve(D, arr)
print(‘方程组的解为:’, x)
![在这里插入图片描述](https://img-blog.csdnimg.cn/20200318151010467.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2x4bG92ZXdlbmppbmc=,size_16,color_FFFFFF,t_70)
![在这里插入图片描述](https://img-blog.csdnimg.cn/20200318151017219.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2x4bG92ZXdlbmppbmc=,size_16,color_FFFFFF,t_70)