9.5.1
fsolve
with the option
dnewton_solver
Input:
fsolve([x
^
2+y-2,x+y
^
2-2],[x,y],[2,2],dnewton_solver)
Output:
[1.0,1.0]