function f=TrussNonLinF(x)
y=x(1);x1=x(2);x2=x(3);
f=x1*sqrt(16+y^2)+x2*sqrt(1+y^2);
function [C,Ceq]=TrussNonLinCon(x,sigma)
y=x(1);x1=x(2);x2=x(3);
C(1)=20*sqrt(16+y^2)-sigma*y*x1;
C(2)=80*sqrt(1+y^2)-sigma*y*x2;
Ceq=[];
x=[1 1 1];sigma=10^5;
LB=[1,0,0];UB=[3,inf,inf];
[x,f]=fmincon(@(x)TrussNonLinF(x),x0,[],[],[],[],LB,UB,@(x)TrussNonLinCon(x,sigma),[],sigma)