Solution of Bessel Differential Equation

How can I get the general solution in this way? My command's result appears like the integer order of the Bessel differential equation.

ODE2 := diff(T[1](t), t, t) + 2*diff(T[1](t), t)/t + (Omega^2/(a^2*t^(4/3)) - 20/(9*t^2))*T[1](t)

generalsol := dsolve(ODE2, T[1](t))

Thank you.

Answer 1

[Walter Roberson]

syms T_1(t)
syms a Omega positive
ODE2 = diff(T_1(t), t, t) + 2*diff(T_1(t), t)/t + (Omega^2/(a^2*t^(4/3)) - 20/(9*t^2))*T_1(t)

ODE2 = 

generalsol = dsolve(ODE2)

generalsol = 

Note that is a constant, and so can be folded into and

https://www.mathworks.com/help/matlab/ref/bessely.html