clc
c=30
b=15
a=10
Ka=0.3872 // from table 13.8
H=4
G=15
Pa=G*H^2*Ka/2
printf('The active force per unit length Pa = %f kN/m\n',Pa)
printf(' The resultant will act at a vertical distance equal to H/3 = 4/3 = 1.33 m above \n the bottom of the wall and will be inclined at an angle of 15to the back face of the wall.')