x = (p*q)/(2*p_0) * sqrt(1+(4*p*g*h)/(p_0))+1