A=[0 0 0];b=[0];Aeq=[0 0 0];beq=[0];lb=-1000*[1 1 1]';ub=1000*[1 1 1 ]';
x0=[1 1 1]';
[x,fval,exitflag,output,lambda,grad,hessian] = fmincon(@objfun,x0,A,b,Aeq,beq,lb,ub,@nonlcon)
%[x,fval,exitflag,output,lambda,grad,hessian] = fmincon(@objfun,x0,@nonlcon)
