এসসিপি fsolve

from scipy import optimize

def equation(x):
    return [
        x[0] ** 2 + x[1] - 5,
        x[0] + x[1] ** 2 - 5
    ]

result = optimize.fsolve(
    equation,
    [1, 1]
)

print(result)