On Oct. 26, 2011, I wrote:
------------------
>(%i26) f(integrate(bessel_k(2,%i*y),y,1,100));
>
>(%o26) 1.736293756153568-2.1367387806763253E+42*%i
>----------------------------------------------------------------
>Where does the 10^42 come from??
----------------------------
I should have included some quad_qag output:
(%i27) quad_qag(imagpart(bessel_k(2.0,%i*y)),y,1,100,3,limit=700);
(%o27) [1.62942476406395,1.1838977094676129E-9,217,0]
Ted