x0 = [-1,1];     % Make a starting guess at the solution
options = optimset('LargeScale','off');
[x, fval] = ... 
fmincon(@objfun2,x0,[],[],[],[],[],[],@confun2,options)
