quad_qag and quad_qags fail with imagpart( bessel_i(1, %i*x)) arg
Subject: quad_qag and quad_qags fail with imagpart( bessel_i(1, %i*x)) arg
From: Stavros Macrakis
Date: Sat, 22 Oct 2011 18:06:22 -0400
On Sat, Oct 22, 2011 at 17:45, Dieter Kaiser <drdieterkaiser at web.de> wrote:
> ...The reason is, that Maxima knows symbolically that bessel_i(1,%i*x) is a
> pure imaginary. Therefore the imagpart automatically simplifies:
>
> (%i1) imagpart(bessel_i(1,%i*x));
> (%o1) bessel_i(1, %i x)
> ...
Yes, b: bessel_i(1,%i*x) is pure imaginary. But the *imagpart* is real,
namely b/%i = -%i*b.
-s