From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 548/CH6/EX6.4.a/6_4a.JPG | Bin 0 -> 21079 bytes
 548/CH6/EX6.4.a/6_4a.sce |  22 ++++++++++++++++++++++
 2 files changed, 22 insertions(+)
 create mode 100755 548/CH6/EX6.4.a/6_4a.JPG
 create mode 100755 548/CH6/EX6.4.a/6_4a.sce

(limited to '548/CH6/EX6.4.a')

diff --git a/548/CH6/EX6.4.a/6_4a.JPG b/548/CH6/EX6.4.a/6_4a.JPG
new file mode 100755
index 000000000..7cfdbea8e
Binary files /dev/null and b/548/CH6/EX6.4.a/6_4a.JPG differ
diff --git a/548/CH6/EX6.4.a/6_4a.sce b/548/CH6/EX6.4.a/6_4a.sce
new file mode 100755
index 000000000..9830fdde0
--- /dev/null
+++ b/548/CH6/EX6.4.a/6_4a.sce
@@ -0,0 +1,22 @@
+pathname=get_absolute_file_path('6_4a.sce')
+filename=pathname+filesep()+'6_4a_data.sci'
+exec(filename)
+clf();
+V=linspace(20,300,500);
+i = 1;
+Cl = 0;Cd = 0;Cl_Cd =0;Thrust = 0;
+while(i<=length(V)) 
+    Cl(i) = 2*W/(D*S*V(i)^2);
+    Cd(i) = Cdo + Cl(i)^2/(%pi*e*AR);
+    Cl_Cd(i) = Cl(i)/Cd(i);
+    Thrust(i) = W/Cl_Cd(i)/1000;
+    Power(i)=Thrust(i)*V(i);
+    Pa(i)=P*Pf*746/1000;
+   i = i+1; 
+end 
+xlabel("Velocity (m/s)");
+ylabel("Power (KN-m/s)");
+plot2d(V,Power,4);
+plot2d(V,Pa,5);
+disp("the intersection of both graph shows maximum velocity of CP-1 at sea level which is arround 81m/s.")
+
-- 
cgit