function  [c,ceq,gc,gceq] = myExpNLCon(x)
c=[];

ceq(1)=1000*sin(-x(3) - 0.25) + 1000*sin(-x(4)-0.25) -x(1)+1000;
ceq(2)=1000*sin(x(3) - 0.25) + 1000*sin(x(3)-x(4)-0.25) -x(2)+1000;
ceq(3)=1000*sin(x(4) - 0.25) + 1000*sin(x(4)-x(3)-0.25)+1000;

eig(rand(200));    % code to take up time

if nargout > 2
gc=[];
gceq=[
  -1,  0, -1000*cos(-x(3)-0.25), -1000*cos(-x(4)-0.25)
   0, -1, -1000*cos(x(3)-0.25)+1000*cos(x(3)-x(4)-0.25), -1000*cos(x(3)-x(4)-0.25)
   0,  0, -1000*cos(x(4)-x(3)-0.25), 1000*cos(x(4)-0.25)+1000*cos(x(4)-x(3)-0.25) 
  ];
end
