close all
clear all

C.rho = 1025;

C.lt = .5;  %Tail Length
C.I = 51;   %Total Inertia
C.cw = 12.22*.0254;
C.Lw = 12*2.5;
C.Lt = 1.5;
C.ct = 12.22*.0254;
C.it = 10*pi/180;


tspan = [0 40];
alpha0 = [0 0];
[t,alpha] = ode45(@(t,alpha) RHS(t,alpha,C), tspan, alpha0);


for i = 1:length(t)
   [th(i) U(i)] = IncFlow(t(i)); 
end

figure
plot(t,alpha(:,1)*180/pi,t,th*180/pi,t,(alpha(:,1)-th')*180/pi);