function ret=o(xu)
  ret = [ 2 .* xu(1) .* xu(1) + xu(2) - 15; -5 .* xu(1) .* xu(1) + 8 * xu(2) - 1 = 0 ];
endfunction

[x2,spr] = fsolve([-1;1], o);
[x3,spr] = fsolve([1;3], o);