diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /1457/CH9/EX9.6 | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '1457/CH9/EX9.6')
| -rwxr-xr-x | 1457/CH9/EX9.6/9_6.sce | 23 |
1 files changed, 23 insertions, 0 deletions
diff --git a/1457/CH9/EX9.6/9_6.sce b/1457/CH9/EX9.6/9_6.sce new file mode 100755 index 000000000..eda956d11 --- /dev/null +++ b/1457/CH9/EX9.6/9_6.sce @@ -0,0 +1,23 @@ +clc
+//Initialization of variables
+P1=10 //psia
+T1=460+40 //R
+R=1715
+k=1.4
+V1=1400 //fps
+//calculations
+rho1=P1/(R*T1)
+c1=sqrt(k*R*T1)
+M1=V1/c1
+P2= P1 * (2*k*M1^2 - (k-1))/(k+1)
+V2 =V1*((k-1)*M1^2 +2)/((k+1)*M1^2)
+rho2=rho1*V1/V2
+T2=P2/rho2/R
+P22=122.5
+V22=286
+T22=328
+//results
+printf("Pressure at point 2 = %.1f psia and %.1f N/m^2",P2,P22)
+printf("\n Velocity at point 2 = %d fps and %d m/s",V2,V22)
+printf("\n Temperature at point 2 = %d R and %d K",T2,T22)
+disp("Similarly it can be done for SI units")
|