Message Boards Message Boards


How to obtain the right solution for modified Bessel ODE?

Posted 2 months ago
1 Reply
0 Total Likes

Hi, I am trying to solve the modified Bessel ODE of real variables:

DSolve[u''[r] + u'[r] / r - a * u[r] ⩵ 0, u[r], r, Assumptions → {a ≥ 0}]

but I obtain the result in terms of the BesselJ and BesselY functions of imaginary arguments:

u[r] → BesselJ[0, ⅈ a^1/2 r] C[1] + BesselY[0, -ⅈ a^1/2 r] C[2]. 

However, I would like to verify that the solution is, in fact, a linear combination of the BesselI and BesselK functions of real arguments. Calling FunctionExpand[BesselJ[0, ⅈ a^(1/2) r] gives indeed BesselI[0, a^(1/2) r], but I cannot find a way to prove that the second component, BesselY[0, -ⅈ a^1/2 r] , would be equivalent to BesselK[0, a^1/2 r], or would contribute such a term. FunctionExpand does not give such a result. Is there any working method to obtain the desired result? Lesław.


 FullSimplify[BesselY[0, -I Sqrt[a] r] // FunctionExpand // Re, 
  Assumptions -> {a > 0, r > 0}]

 (*-((2 BesselK[0, Sqrt[a] r])/\[Pi])*)
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
or Discard

Group Abstract Group Abstract