//Problem 1.3
h1=2.5
h2=1.5
s1=1
s2=0.8
w=9810
p2=s2*w*h2 //Pressure intensity at interface
p1=s1*w*h1
p=p1+p2
disp(p2, "pressure intensity at interface(N/m2)")
disp(p, "pressure intensity at bottom(N/m2)")