function [n,xv,fvec,iflag]=c05fj(n,xv,fvev,iflag)
x = xv(1)   ;
y = xv(2)   ;
z = xv(3)   ;
  fvec = [z^2-x^2-y^2+2*x+2*z-1, y*z-y-x+1, -2+x+y]   ;
