function [M] = TotalMoment(alpha,theta,V)


[rho S Cp C ] = Constants();
C4 = C/4;

qinf = .5*rho*V.^2;
M0 = S*qinf*C^2;
F0 = S*qinf*C;

aoa = alpha-theta;

[cl cd cm]=Coef0015(aoa*pi/180);

cl = 0.75*cl;  %Rough adjustment for finite span.

L = cl*F0;
D = cd*F0;
M = cm*M0;

M_wing = M+L*(Cp-C4).*cos(aoa*pi/180)+D*(Cp-C4).*sin(aoa*pi/180);

M_spring = GangedSpring(alpha);  %Spring restoring force estimate from March2017 deployment

M = M_wing+M_spring;

%disp(['alpha = ' num2str(alpha)]);
end
