blob: e92a73568c2f6011954b9f0bb450730485b1b82f (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
clc
ro=16
gama=1.4
T1=310.15 //temperature in kelvin
T2=T1*(ro^(gama-1))
mprintf("T2=%fK\n",T2)//ans vary due to roundoff error
deltas=1.2
CP=1.0047
T3=(%e^(deltas/CP))*T2
mprintf("T3=%fK\n",T3)//ans vary due to roundoff error
q1=CP*(T3-T2)
rc=T3/T2
mprintf("rc=%f\n",rc)//ans vary due to roundoff error
mprintf("q1=%fkJ/kg\n",q1)//ans vary due to roundoff error
Eta=1-(((rc^gama)-1)/((gama*ro^(gama-1))*(rc-1)))
mprintf("Eta=%f",Eta)//ans vary due to roundoff error
|
