summaryrefslogtreecommitdiff
path: root/51/CH1/EX1.3/1_3.sce
blob: e8f5b712f979ea846937661f803083635249802a (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
clc;
clear;
V=0.84;//ft^3
p=50;//psi
T=70;//degree farenheit
atmp=14.7;//psi
//the air density d=P/(RT)
//1ft^2=144 inches^2
d=((p+atmp)*144)/((1716)*(T+460));
disp(d)
//slugs/ft^3
//weight of air
W=d*32.2*V;
//1lb=1 slug.ft/sq sec
disp("lb",W,"W=")
//taking various values of p a graph is plotted between W and p
x= -20:60;
for p= -20: 60
    i=p+21;
    y(1,i)=((p+atmp)*144/((1716)*(T+460)))*32.2*V;
    
end
plot(x,y)
xtitle('W vs p','p.psi','W,lb')