>
>
> sum(x^(3*k)/factorial(2*k),k,0,oo)
>
> I understand that Sage has limited exploitation of Maxima's hypergeometric
> functionality, and I suspect this is the main issue. Are there any
> conceivable workarounds?
>
>
This actually works out of the box in the latest development release,
proba
>
>
> > which, to me, is a very useful answer. But other sums are simply wrong.
> >
> > k = var('k')
> > sum(x^(2*k)/factorial(2*k),k,0,oo)
>
> I'm working with Maxima 5.33.0. I get
>
> simplify_sum ('sum(x^(2*k)/factorial(2*k),k,0,inf));
> => sqrt(%pi)*bessel_i(-1/2,x)*sqrt(x)/sqrt(2)
On 2014-04-23, Karl S wrote:
> #taylor coefficient for erf(3x)
> a_erf(m) = (3)^(2*m+1)*(-1)^m*2/sqrt(pi)/(factorial(m)*(2*m+1))
>
> #coefficient of chebyshev polynomial
> c_erf_cheb(p) = sum(a_erf(m)*binomial(2*m+1,m-p)*4^-m,m,p,oo).simplify_full
> ()
>
> Here the function c_erf_cheb(p) ends up
On Wednesday, April 23, 2014 2:40:46 PM UTC-7, Karl S wrote:
>
>
> I understand that Sage has limited exploitation of Maxima's hypergeometric
> functionality, and I suspect this is the main issue. Are there any
> conceivable workarounds?
>
>
http://trac.sagemath.org/ticket/2516 should basically
>
> But other sums are simply wrong.
>
> k = var('k')
> sum(x^(2*k)/factorial(2*k),k,0,oo)
>
> gives
>
> -1/4*sqrt(2)*sqrt(pi)*x^(3/2)
>
> but the answer should be sinh(x).
>
>
Hmm. That shouldn't be happening, though I wouldn't be surprised if it
didn't turn out as easy as that.
(%i1) load(si