A = [];
B = [];
Aeq = [];
Beq = [];
x0 =[1.5; 3];
x=fmincon(@myfun3,x0,A,B,Aeq,Beq,[],[],@nonlcon3)
