A = [];
B = [];
Aeq = [];
Beq = [];
x0 =[1.5; 3];
[x,fval]=fmincon(@myfun4,x0,A,B,Aeq,Beq,[],[],@nonlcon4)
