bessel_y sign error with integrate?

• Subject: bessel_y sign error with integrate?
• From: Edwin Woollett
• Date: Fri, 21 Oct 2011 15:51:41 -0700

On Oct. 21, 2011, I wrote about bessel_i integrate
sign error. Using the same session, here is an 
apparent sign error when using integrate
with bessel_y(2,x) . (The reasoning is the same.)

first the 'height' (the 'base' is unity):

(%i23) f(bessel_y(2,4));
(%o23) 0.21590359460361

(%i24) integrate (bessel_y(2,y),y,3.5,4.5);
(%o24) (46650016*%pi-292890325)/445824550-(6893566*%pi-5764575)/63526475
(%i25) f(%);
(%o25) -0.57840103823513

So both sign and magnitude error apparently.

Ted Woollett

• Prev by Date: integrate of bessel_i sign error?
• Next by Date: computing generating elements of ZZ_p^* ?
• Previous by thread: bessel.lisp error codes and noisy quad_qagi
• Next by thread: bessel_y sign error with integrate?
• Index(es):
  • Date
  • Thread