1 files changed, 17 insertions, 0 deletions
diff --git a/548/CH9/EX9.1/9_01.sce b/548/CH9/EX9.1/9_01.sce
new file mode 100755
index 000000000..9bdd8718c
--- /dev/null
+++ b/548/CH9/EX9.1/9_01.sce
@@ -0,0 +1,17 @@
+pathname=get_absolute_file_path('9_01.sce')
+filename=pathname+filesep()+'9_01data.sci'
+exec(filename)
+x=poly(0,'x');
+P=x-10*x+9.5;
+t=roots(P);
+V2=%pi*b^2*(Stroke+t)*(10^-6)/4;disp(V2,"V2=%pi*b^2*(Stroke+t)/4");
+V3=V2/r;disp(V3,"V3=V2/r");
+V5=V2;V4=V3;
+Wcomp=(P2*V2-P3*V3)/(1-y);
+disp(Wcomp,"Wcomp=","Wcomp=P2*V2-P3*V3/(1-y);","work done in compression cycle Wcomp:")
+Wpower=(P5*V5-P4*V4)/(1-y);
+disp(Wpower,"Wpower=","Wpower=P5*V5-P4*V4/(1-y);","work done in power stroke Wpower:")
+Pa=6*n*nm*(rpm)*(Wpower-Wcomp)/120;
+disp(Pa,"Pa=n*nm*(rpm)*(Wpower-Wcomp)/120","power available Pa:")
+printf("\Answer:\n")
+printf("\n\Power available from the engine propeller combination: %f J/s \n\n",Pa)