习题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")

结果如下