ineq2_ex3.m

good code for matlab by mili , i than you

function retval = ineq2_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;
S = FLAGW + (FD * L);
Q = (X1.*X1 + X1.*X2 + X2.*X2)/6.0;


retval = S.*Q./INERTIA;