4.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]