solve for e : c = Sqrt[ Sqrt[ (1/4)*( e^2 - b^2 )^2 + (p^2)*(e^2) ] + (1/2)* (e^2 - b^2) ]