octave:13> function y=f(x) > y(1)=x(1)-x(2)**2 > y(2)=x(1)-x(2) > endfunction octave:14> [x, info]=fsolve("f", [0.5; 2]); octave:15> x x = 1.0000 1.0000 octave:16> info info = 1