习题2.9
import sympy as sp
# 定义变量
x, y = sp.symbols('x y')
# 定义方程组
equation1 = sp.Eq(x**2 - y - x, 3)
equation2 = sp.Eq(x + 3*y, 2)
# 解方程组
solutions = sp.solve((equation1, equation2), (x, y), dict=True)
print("符号解：")
for sol in solutions:
print(sol)
from scipy.optimize import fsolve
import numpy as np
# 定义方程组函数
def equations(vars):
x, y = vars
return [x**2 - y - x - 3, x + 3*y - 2]
# 初始猜测值
initial_guess = [1, 0]
# 解方程组
solution = fsolve(equations, initial_guess)
print("数值解：")
print(f"x = {solution[0]}, y = {solution[1]}")
print("学号：3008")结果如下
