// calculate angle to which tube is incliend to horizontal
clc;
P=9.81;
g=9.81;
dm=0.864*10^3;
R=4*10^-3;
d=2;
D=20;
th=asind(P/(g*dm*R*(1+(d/D)^2)));
disp(th,'angle to which tube is incliend to horizontal(degree)')