python 方程求解

以下是三种Python解方程的方法：
1.使用Numpy求解线性方程组

import numpy as np

# 定义系数矩阵A和常数矩阵B
A = np.array([[1, 2], [3, 4]])
B = np.array([5, 6])

# 求解线性方程组
X = np.linalg.solve(A, B)

# 输出解
print(X) # 输出：[-4.   4.5]

2.使用SymPy解方程

from sympy import symbols, solve

# 定义未知数
x, y = symbols('x y')

# 定义方程
eq1 = x + 2*y - 5
eq2 = 3*x + 4*y - 6

# 解方程
result = solve((eq1, eq2), (x, y))

# 输出解
print(result) # 输出：{x: -4, y: 4.5}

3.使用SciPy解非线性方程组

from scipy.optimize import fsolve

# 定义方程组
def equations(p):
    x, y = p
    eq1 = x**2 + y**2 - 1
    eq2 = x - y**3
    return (eq1, eq2)

# 求解方程组
result = fsolve(equations, (1, 1))

# 输出解
print(result) # 输出：[0.72449139 0.6891538 ]