#3426: [with new patch, needs review] bessel_K function is broken
-----------------------------------------------+----------------------------
 Reporter:  bober                              |        Owner:  gfurnish
     Type:  defect                             |       Status:  new     
 Priority:  major                              |    Milestone:  sage-3.2
Component:  calculus                           |   Resolution:          
 Keywords:  bessel, bessel_K, editor_gfurnish  |  
-----------------------------------------------+----------------------------
Comment (by AlexGhitza):

 Thanks for catching this.  It actually comes from a bug in Pari:

 {{{
 ? besselj(0, 0)
 %1 = 1.000000000000000000000000000
 ? besselj(0.E-19, 0)
   *** besselj: gpow: 0 to a non positive exponent.
 }}}

 I've reported it upstream, but I will post  a patch with a workaround
 while we wait.

-- 
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/3426#comment:29>
Sage <http://sagemath.org/>
Sage - Open Source Mathematical Software: Building the Car Instead of 
Reinventing the Wheel
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups 
"sage-trac" group.
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to [EMAIL PROTECTED]
For more options, visit this group at 
http://groups.google.com/group/sage-trac?hl=en
-~----------~----~----~----~------~----~------~--~---

Reply via email to