More info on the wrong evaluation from Maxima:
This seems to be happening when the second argument of spherical_bessel_J
is an RR, but not if it is in ZZ o QQ.
On Friday, July 19, 2013 12:18:51 PM UTC+2, Jesús Torrado wrote:
>
> Hi there,
>
> Have a look at the second line:
>
> sage: maxima('spherical_bessel_j(50,9.5)')
> 18.003620332195756l-33
> sage: spherical_bessel_J(50,9.5, algorithm="maxima") # default
> algorithm
> 18.0036203322*l - 33
> sage: spherical_bessel_J(50, 9.5, algorithm="scipy")
> 1.8003620332195507e-32
>
> It seems Sage is not making sense out of that Maxima expression, is it?
>
> Also, notice that:
>
> sage: spherical_bessel_J(50,10000, algorithm="maxima").n()
> 0.0000424207788467889
> sage: spherical_bessel_J(50,10000, algorithm="scipy")
> nan
>
> which is an open issue of scipy as of today:
> https://github.com/scipy/scipy/issues/1641
>
> Cheers,
> Jesús Torrado
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-support.
For more options, visit https://groups.google.com/groups/opt_out.
