Bug [1954846] bessel_i(1/2,0) -> divide by zero error

On Sat, Jul 5, 2008 at 4:07 AM, Dieter Kaiser <drdieterkaiser at web.de> wrote:
> Might it be possible to introduce to Maxima the correct handling of
> infinities, so we can define the singularities of functions more completly?

That would be nice, but several attempts over the years have not
resulted in usable code....  Also, it is not clear what the correct
handling is.

> (%i25) tan(%pi/2);
> (%o25) tan(%pi/2)             ; we get a noun

This seems to be new behavior.  I am not sure it is the best behavior,
but it is better than returning inf/minf/infinity.

> (%i29) tan(%pi/2),numer;
> (%o29) 1.6331778728383844E+16 ; a result!

Since float(%pi/2) is not exactly %pi/2, tan(float(%pi/2)) will not be
exactly infinite.  This is correct behavior.  Anyone using approximate
numbers needs to be aware of cases like this.

         -s