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