//Example 3.20
clc;
clear;
close;
format('v',7);
//Given data :
a=2;//meter
b=1;//meter
d=2;//meter
w=9.81;//kN/m^2
xbar=2+a/2;//meter
A=a*b;//m^2
P=w*A*xbar;//kN
disp(P,"Total Pressure in kN : ");
IG=(b*d^3)/12;//in m^4
h_bar=IG/A/xbar+xbar;//in meter
disp(h_bar,"Position of centre of pressure in meter : ");