From 32b6a6e80d22a6d6e8449821d3117252de3e89ba Mon Sep 17 00:00:00 2001
From: Marius Peter <blendoit@gmail.com>
Date: Mon, 15 Jul 2019 13:29:36 -0700
Subject: commence class-ification of gui

---
 wing_scripts/get_z.m | 34 ++++++++++++++++++++++++++++++++++
 1 file changed, 34 insertions(+)
 create mode 100644 wing_scripts/get_z.m

(limited to 'wing_scripts/get_z.m')

diff --git a/wing_scripts/get_z.m b/wing_scripts/get_z.m
new file mode 100644
index 0000000..5387b52
--- /dev/null
+++ b/wing_scripts/get_z.m
@@ -0,0 +1,34 @@
+function z = get_z(x,u)
+
+
+
+if (x < 0 )
+    disp('invalid X')
+end
+
+M = 0.02;
+P = 0.4;
+T = 0.12;
+a0 = 0.2969;
+a1 = -0.126;
+a2 = -0.3516;
+a3 = 0.2843;
+a4 = -0.1015;
+
+if x <P
+    z_camber = M/P^2*(2*P*x - x^2);
+else
+    z_camber = (M/(1-P)^2)*(1 - 2*P +2*P*x - x^2);
+end
+
+%z_camber = M/P^2*(2*P*x - x^2);
+z_thickness = (T/0.2)*(a0*x^.5+a1*x+a2*x^2+a3*x^3+a4*x^4);
+
+if u==1
+    z = z_camber + z_thickness;
+else
+    z = z_camber - z_thickness;
+end
+
+
+
-- 
cgit v1.2.3