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: Dieter Kaiser
Date: Sun, 23 Oct 2011 00:01:52 +0200
Am Samstag, den 22.10.2011, 17:36 -0400 schrieb Stavros Macrakis:
> There is a bug in imagpart. imagpart(bessel_i(1,%i*x)) should be -%
> i*bessel_i(1,%i*x), but imagpart is returning bessel_i(1,%i*x).
>
>
> Could you please report this -- see bug_report() for instructions.
Yes, this is a bug.
Dieter Kaiser