ineq1.m

good code for matlab by mili , i than you

function retval = ineq1_ex3(X1,X2)
% global statement is used to share same information between 
% various m-files
global ELAS SIGALL TAUALL GAM FS GRAV
global RHO CD FLAGW SPEED LP L DELT

AREA = 0.25* pi*(X1.^2 - X2.^2);
INERTIA = pi*(X1.^4 - X2.^4)/64;
FD = 0.5*RHO*SPEED*SPEED*CD.*X1;
MW = 0.5*FD*L*L;
MF = FLAGW * LP*LP*LP/(3.0*ELAS*I);
SIGBEND = 0.5*(MW + MF).*X1./INERTIA;
SIGW = GAM*GRAV*L;

retval = SIGBEND + SIGW;