x=2;
function y=f(x);
y(1)=x(1)^2-x(2)^2+x(1)+1;
y(2)=x(1)*x(2)+x(2);
endfunction
fsolve("f",[1,2])