解微分方程，求matlab代码

利用龙格库塔法求贝塞尔方程的数值解，这个方程的n为常数，应该怎么表达呢

能再说的具体一些吗？