clc
G=15.7
a=0
b=15
c=30
H=3
Kp=4.977 // from table 13.9
Pp=Kp*G*H^2/2
printf('a)the passive force = %f kN/m\n',Pp)
// for part b
Kp=4.53
Pp=Kp*G*H^2/2
printf(' b)the passive force = %f kN/m\n',Pp)
// for part c
Kp=4.13
Pp=Kp*G*H^2/2
printf(' c)the passive force = %f kN/m\n',Pp)
//for part d
Kp=4.56
Pp=Kp*G*H^2/2
printf(' d)the passive force = %f kN/m\n',Pp)