2-9
from sympy import symbols, Eq, solve
x, y = symbols('x y')
equations = (Eq(x2 - y - x - 3, 0), Eq(x + 3*y - 2, 0))
symbolic_solution = solve(equations, (x, y))
print("符号解:", symbolic_solution)
from scipy.optimize import fsolve
def equations(variables):
x, y = variables
return [x2 - y - x - 3, x + 3*y - 2]
initial_guess = [1, 1]
numeric_solution = fsolve(equations, initial_guess)
print("数值解:", numeric_solution)
print('3006')