summaryrefslogtreecommitdiff
path: root/3785/CH2/EX2.1/Ex2_1.sce
blob: 1f3772062621b8921bc31dde377360773ac33205 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
// Example 2_1
clc;funcprot(0);
// Given data
p_0=1.0133*10^5;// The sea level pressure in Pa
alpha=1.2*10^-4;// m^-1

// Calculation
z=0;// km
minusdelp=alpha*p_0*exp(-alpha*z);// The pressure force per unit volume in N/m^3
printf("\nAt z=0,The pressure force per unit volume -delp=(%0.2f N/m^3)i_z",minusdelp);
z=5;// km
minusdelp=alpha*p_0*exp(-alpha*z*10^3);// The pressure force per unit volume in N/m^3
printf("\nAt z=5,The pressure force per unit volume -delp=(%0.3f N/m^3)i_z",minusdelp);