V = 70 %m/s
alpha = 5;
beta = 0;
p = 0;
q = 0;
r = 0;
psi = 0;
theta = 0;
phi = 0;
xe = 0;
ye = 0;
H = 1000;

x = [V alpha beta p q r psi theta phi xe ye H]';

n = 4800;
throttle = 65; %percentuale di manetta
uprop = [n throttle];

rho = 1.1117;
ps = 898.76;
T = 281.65;
mu = 1.76e-5;
g = 9.81;
yatm = [rho ps T mu g];

a = 336.44;
M = V/a;
qdyn = 0.5*rho*V^2;
yad1 = [a M qdyn];

eta = 0.7;
rhoSL = 1.225;
TSL = 288.15;
psi = (yatm(1)/rhoSL)*sqrt(TSL/yatm(3));

