library nlsys;
nlset;

proc F(x);
  local f1,f2,e;
  f1 = 0.5*sin(x[1]*x[2]) - x[2]/(4*pi) - x[1]/2;
  f2 = (1 - 1/(4*pi))*(exp(2*x[1] - 1) - 1) + x[2]/pi - 2*x[1];
  retp(f1|f2);
endp;

x0 = { 0.2, 0.3 };
output file = nl6.out reset;
call nlprt(nlsys(&f,x0));
output off;