I guess there is a limit bug here somewhere.  With MultiSeries things are OK...

L:=(x*BesselI(k-1, x)+x*BesselI(k+1, x))*BesselI(-k, x);
L := (x BesselI(k - 1, x) + x BesselI(k + 1, x)) BesselI(-k, x)
simplify(limit(L,x=0)) assuming k>0;
0
with(MultiSeries): simplify(limit(L,x=0)) assuming k>0;
2 sin(Pi k)
-----------
Pi

BesselLimit.mw

---

G A Edgar


Please Wait...