1 files changed, 22 insertions, 0 deletions
diff --git a/51/CH1/EX1.6/1_6.sce b/51/CH1/EX1.6/1_6.sce
new file mode 100755
index 000000000..aecc14046
--- /dev/null
+++ b/51/CH1/EX1.6/1_6.sce
@@ -0,0 +1,22 @@
+clc;
+clear;
+p1=14.7;//psi(abs)
+V1=1;//ft^3
+V2=0.5;//ft^3
+//for isentropic compression, (p1(d1^k))=(p2/(d2^k))
+//volume*density=constant(mass)
+ratd=V1/V2;
+p2=((ratd)^1.66)*p1;//psi(abs)
+disp("psi(abs)",p2,"final pressure p2=")
+
+i=1;
+ratV=0.01:0.01:1.0;
+
+for j=0.01:0.01:1.0
+    pres(i)=p1/((j)^1.66);
+    i=i+1;
+    
+end
+
+plot2d(ratV,pres,rect=[0,0,1,1000])
+xtitle('p2 vs V2/V1','V2/V1','p2 psi')