mmSaeule = @(dP,rho) dP*100/(g*rho)*1000; % rho in kg/m^3
p_mmW_2 = mmSaeule(p_2-p_L,rho_W*1000)
Undefined function or variable 'rho_W'.
p_mmW_3 = mmSaeule(p_3-p_L,rho_W*1000)
p_mmHg_2 = mmSaeule(p_2-p_L,rho_Hg*1000)
p_mmHg_3 = mmSaeule(p_3-p_L,rho_Hg*1000)
hk = @(gamma,alpha,rho,r) (2*gamma*cosd(alpha))/(rho*g*r);
hk_W = hk(gamma_W,alpha_W,rho_W*1000,d/2*0.001)
hk_Hg = hk(gamma_Hg,alpha_Hg,rho_Hg*1000,d/2*0.001)