function [n,xv,fvec,iflag]=c05fj(n,xv,fvev,iflag)
x1 = xv(1)   ;
x2 = xv(2)   ;
  fvec = [x1*x2-x2^3-1, x1^2*x2+x2-5]   ;
