2.9

from sympy import symbols, Eq, solve

x, y = symbols('x y')

equations = (Eq(x**2 - 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 [x**2 - y - x - 3, x + 3*y - 2]

initial_guess = [1, 1]

numeric_solution = fsolve(equations, initial_guess)
print("数值解:", numeric_solution)

3022