clear all
close all

rho = 1025;  %kg/m^3

U = 1;
T = 10;
b = .25;  
alpha0 = 0*10;  phi0 = 90;  %Degrees
h0 = 1*1.3;  phi2 = 0;

motion = FF_Motion(U,alpha0*pi/180,phi0*pi/180,h0,phi2*pi/180,T,b);
clear U T b alpha0 phi0 h0 phi2;

AOA = FF_AngleOfAttack(motion);


figure
plot(motion.t,imag(motion.h));
xlabel('time (s)');
ylabel('h (m)');

figure
plot(motion.t,imag(motion.alpha)*180/pi);
xlabel('time (s)');
ylabel('alpha (deg)');

figure
plot(motion.t,AOA.Ubar);

figure
plot(motion.t,AOA.gamma*180/pi);
xlabel('time (s)');
ylabel('gamma (deg)');

figure
plot(motion.t,AOA.beta*180/pi);
xlabel('time (s)');
ylabel('beta (deg)');
