sympy dsolve TypeError: Invalid Input

代码：
from sympy import *
f=Function('f')
t,p,m=symbols('t p M')
eq=tdiff(f(t),(t,2))+pdiff(f(t),t)+tm**2f(t)
dsolve(eq,f(t))

dsolve不能给出该方程的解，它的解应该为贝塞尔函数乘t的(1-p)/2次方。

dsolve(eq,t)
这样吧