📅  最后修改于: 2023-12-03 14:55:58.284000             🧑  作者: Mango
本程序旨在解决以下四次方程:
(x2 – 43)(x2 – 4x – 5) = 5(x + 1) – 3(x-5)
import numpy as np
def solve_equation():
# 获取参数值
a = 1
b = -4
c = -48
d = -50
e = 118
# 构造系数矩阵A和常数向量b
A = np.array([[a, 0, b, 0, c], [0, a, b, c, 0], [1, 0, -4, -1, -5], [0, 1, -4, -5, 0], [5, -3, 0, 0, -5]])
b = np.array([-e, -e, -6, 15, -4])
# 调用numpy.linalg.solve方法求解线性方程组
x = np.linalg.solve(A, b)
# 返回解析结果
return x
# 输出解析结果
print(f"方程的实数解为:{solve_equation()}")
首先,将式子化简:
x4 - 4x3 - 38x2 + 176x - 215 = 0
使用numpy库中的linalg.solve方法求解线性方程组,将方程的系数和常数向量构造成矩阵方程,得到方程的实数解。最后返回解析结果。