Thanks for your response. I tried this but still, the result is too long. How can I further improve it?
If sol is the formula for X[Y], for example, I would try this as a start:
sol
X[Y],
Collect[TrigToExp@Numerator[Together[sol]] /. E^u_ :> E^Collect[u, Y], Y^_.*E^(_. + _*Y), Style[Simplify[#], RandomColor[]] &]
It needs more work, but I have to go now.
