ASKSAGE: Sage Q&A Forum - Latest question feedhttp://ask.sagemath.org/questions/Q&A Forum for SageenCopyright Sage, 2010. Some rights reserved under creative commons license.Wed, 24 Apr 2019 09:19:34 -0500How to increase maxterms for hypergeometric?http://ask.sagemath.org/question/46337/how-to-increase-maxterms-for-hypergeometric/I'm trying to compute a (bunch of) hypergeometric functions for which I get a NoConvergence error
e.g.
hypergeometric([4.14 + 15*I, -3.14 + 15*I],[1. - 1.12e7*I], -500000)
The suggestion of the error message is to try increasing 'maxterms'. However, I don't seem to find a way to do this as e.g.` hypergeometric([4.14 + 15*I, -3.14 + 15*I],[1. - 1.12e7*I], -500000,maxterms=10^6)` does not work.
(I could be missing something very basic, I'm rather new to sage)virtual_neutrinoWed, 24 Apr 2019 09:19:34 -0500http://ask.sagemath.org/question/46337/A simple hypergeometric function fails.http://ask.sagemath.org/question/24677/a-simple-hypergeometric-function-fails/There is a nice method to compute the Narayana polynomials.
With Maple we can write
P := n -> simplify(hypergeom([-n,-n+1], [2], 1/x));
seq(expand(x^k*P(k)), k=0..5);
and get the answer
1, x, x^2+x, x^3+3*x^2+x, x^4+6*x^3+6*x^2+x, x^5+10*x^4+20*x^3+10*x^2+x.
Trying the same with Sage
P = lambda n: simplify(hypergeometric([-n,-n+1],[2], 1/x))
[expand(x^k*P(k)) for k in (0..5)]
gives the answer
[1, x, ..., x^n*hypergeometric((-n, -n-1),(2),1/x)]
This is really disappointing. Is there a workaround?
Peter LuschnyWed, 29 Oct 2014 10:43:29 -0500http://ask.sagemath.org/question/24677/