diff options
Diffstat (limited to '797/CH3/EX3.4.s')
| -rw-r--r-- | 797/CH3/EX3.4.s/3_04_solution.sce | 10 |
1 files changed, 10 insertions, 0 deletions
diff --git a/797/CH3/EX3.4.s/3_04_solution.sce b/797/CH3/EX3.4.s/3_04_solution.sce new file mode 100644 index 000000000..3aefd7c9f --- /dev/null +++ b/797/CH3/EX3.4.s/3_04_solution.sce @@ -0,0 +1,10 @@ +//Soultion 3-04 +WD=get_absolute_file_path('3_04_solution.sce'); +datafile=WD+filesep()+'3_04_example.sci'; +clc; +exec(datafile) +P_1 = rho_0 * g * h_1; //Gauge pressure at the bottom of surface zone [Pa] +P_2 = P_1 + rho_0 * g * 4 *H / %pi * asinh(tan( %pi * H / (4 * H))); //After integrating w.r.t depth s +P_2 = P_2 / 1000; //conversion from [Pa] to [kPa] +//result +printf("Pressure at the bottom of Gradient layer is %1.1f kPa",P_2); |