blob: f795ecf3708b5f86bb8eb47facd2b2c36459f193 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
clc
//initialization of variables
P1=30 //lb/in^2
T1=300+460 //R
T2=60 +460 //R
cp=0.25 //B/lbm F
R=53.3 //ft-lb/lbm R
//calculations
Q=cp*(T2-T1)
du=(cp-R/778)*(T2-T1)
W=778*(Q-du)
function [ds]=c(T)
ds=cp/T
endfunction
S=intg(T1,T2,c)
//results
printf("Change in entropy = %.3f B/lbm R",S)
|