1 files changed, 28 insertions, 0 deletions
diff --git a/Prof_Arya/Trim_Cessna_Body.m b/Prof_Arya/Trim_Cessna_Body.m
new file mode 100644
index 0000000..79be1d0
--- /dev/null
+++ b/Prof_Arya/Trim_Cessna_Body.m
@@ -0,0 +1,28 @@
+function F = Trim_Cessna_Body(x)
+
+m = 1043.26;
+S_ref = 16.1651;
+d = 1.225;
+V = 60;
+g = 9.8;
+
+C_L_0 = 0.25;
+C_L_alpha = 4.47;
+C_L_delta_e = 0.3476;
+        
+C_D_0 = 0.036;
+C_D_alpha = 0.13;
+C_D_delta_e = 0.06;
+        
+C_M_0 = -0.02;
+C_M_alpha = -1.8;
+C_M_delta_e = -1.28;
+
+C_L = C_L_0 + C_L_alpha*x(1) + C_L_delta_e*x(2);
+% C_D = C_D_0 + C_D_alpha*abs(x(1)) + C_D_delta_e*abs(x(2));
+C_D = C_D_0 + 0.3*C_L^2;
+C_M = C_M_0 + C_M_alpha*x(1) + C_M_delta_e*x(2);
+
+F = [C_M;
+    0.5*d*V*V*S_ref*(-C_D*cos(x(1)) + C_L*sin(x(1))) - m*g*sin(x(1)) + x(3);
+    0.5*d*V*V*S_ref*(-C_D*sin(x(1)) - C_L*cos(x(1))) + m*g*cos(x(1))];
+\ No newline at end of file