sqrt(5 E - BesselI[0, 2] - BesselI[0, 1] - (2 E - BesselI[0, 2])^2)