x = -(P_0/(p*g))/2 + sqrt( ((P_0/(p*g))/2)^2 + (P_0*h/(p*g))