A&S 13.6.13 tells us that
1F1(1,2,x) = 2 * exp(x/2) * sinh(x/2) / x
but
(%i72) hgfred([1],[2],x); (%o72) (sqrt(%pi)*bessel_i(1/2,x/2)*%e^(x/2))/sqrt(x)
I think (%o72) isn't wrong, but I don't think it's what we want.
Logged In: YES user_id=28849 Originator: NO
ev(%o72,besselexpand:true) will produce what you want.
We could also add this as a special case, if that's what people really want.
Logged In: YES user_id=895922 Originator: YES
This isn't a bug and there is a workaround. I'll close the bug.
Log in to post a comment.
Logged In: YES
user_id=28849
Originator: NO
ev(%o72,besselexpand:true) will produce what you want.
We could also add this as a special case, if that's what people really want.
Logged In: YES
user_id=895922
Originator: YES
This isn't a bug and there is a workaround. I'll close the bug.