1 files changed, 15 insertions, 0 deletions
diff --git a/548/CH4/EX4.20/4_20.sce b/548/CH4/EX4.20/4_20.sce
new file mode 100755
index 000000000..2d375c7eb
--- /dev/null
+++ b/548/CH4/EX4.20/4_20.sce
@@ -0,0 +1,15 @@
+pathname=get_absolute_file_path('4_20.sce')
+filename=pathname+filesep()+'4_20data.sci'
+exec(filename)
+To=Te*A
+disp(To,"To=","To=Te*(1+(y-1)*Me^2/2)","let reservoir temperature required is To ")
+Po=Pe*A^(y/(y-1));
+disp(Po,"Po=","Po=Pe*((1+(y-1)*Me^2/2))^y/y-1","let reservoir pressure required is Po ")
+r=sqrt((2*A/(y+1))^((y+1)/(y-1))/Me^2)
+disp(r,"Ae/At=","Ae/At=sqrt((2*(1+(y-1)*Me^2/2)/(y+1))^((y+1)/(y-1))/Me^2","Area ratio required is equal to")
+printf("\Answer:\n")
+printf("\n\ required reservoir temperature: %f K\n\n",To)
+printf("\n\ required reservoir pressure: %f N/m^2\n\n",Po)
+printf("\n\ required Area Ratio: %f \n\n",r)
+
+