10.5.4
fsolve
with the option
newtonj_solver
Example.
Input:
fsolve([x^2+y-2,x+y^2-2],[x,y],[0,0],newtonj_solver)
Output:
⎡
⎣
1.0,1.0
⎤
⎦