4.5.4
fsolve
with the option
newtonj_solver
Input :
fsolve([x
^
2+y-2,x+y
^
2-2],[x,y],[0,0],newtonj_solver)
Output :
[1.0,1.0]