Could you open issues for these (or better yet, pull requests)? Otherwise they will be forgotten.
Aaron Meurer On Thu, Aug 16, 2012 at 6:09 AM, someone <[email protected]> wrote: >> Similar issues arise for cos, sinh and cosh. > > To be clear, there are also two bugs here. Sinh > works like sin: > > In [27]: integrate(sinh(x**n)/x**n, x, meijerg=True) > Out[27]: x*gamma(1/(2*n))*hyper((1/(2*n),), (3/2, 1 + 1/(2*n)), > x**(2*n)/4)/(2*n*gamma(1 + 1/(2*n))) > > In [28]: combsimp(_) > Out[28]: x*hyper((1/(2*n),), (3/2, 1 + 1/(2*n)), x**(2*n)/4) > > In [29]: hyperexpand(_) > Out[29]: x*hyper((1/(2*n),), (3/2, 1 + 1/(2*n)), x**(2*n)/4) > > In [30]: _.subs(n,1) > Out[30]: x*hyper((1/2,), (3/2, 3/2), x**2/4) > > In [31]: hyperexpand(_) > Out[31]: Shi(x) > > But the other two fail with a zoo comming (probably) > from wrong gamma prefactors: > > In [23]: integrate(cos(x**n)/x**n, x, meijerg=True) > Out[23]: x*x**(-n)*gamma(-1/2 + 1/(2*n))*hyper((-1/2 + 1/(2*n),), (1/2, 1/2 + > 1/(2*n)), -x**(2*n)/4)/(2*n*gamma(1/2 + 1/(2*n))) > > In [24]: combsimp(_) > Out[24]: x*x**(-n)*gamma(-1/2 + 1/(2*n))*hyper((-1/2 + 1/(2*n),), (1/2, 1/2 + > 1/(2*n)), -x**(2*n)/4)/(2*n*gamma(1/2 + 1/(2*n))) > > In [25]: _.subs(n,1) > Out[25]: zoo*hyper((0,), (1/2, 1), -x**2/4) > > In [26]: hyperexpand(_) > Out[26]: zoo > > > and > > > In [32]: integrate(cosh(x**n)/x**n, x, meijerg=True) > Out[32]: x*x**(-n)*gamma(-1/2 + 1/(2*n))*hyper((-1/2 + 1/(2*n),), (1/2, 1/2 + > 1/(2*n)), x**(2*n)/4)/(2*n*gamma(1/2 + 1/(2*n))) > > In [33]: combsimp(_) > Out[33]: x*x**(-n)*gamma(-1/2 + 1/(2*n))*hyper((-1/2 + 1/(2*n),), (1/2, 1/2 + > 1/(2*n)), x**(2*n)/4)/(2*n*gamma(1/2 + 1/(2*n))) > > In [34]: hyperexpand(_) > Out[34]: x*x**(-n)*gamma(-1/2 + 1/(2*n))*hyper((-1/2 + 1/(2*n),), (1/2, 1/2 + > 1/(2*n)), x**(2*n)/4)/(2*n*gamma(1/2 + 1/(2*n))) > > In [35]: _.subs(n,1) > Out[35]: zoo*hyper((0,), (1/2, 1), x**2/4) > > In [36]: hyperexpand(_) > Out[36]: zoo > > -- > You received this message because you are subscribed to the Google Groups > "sympy" 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/sympy?hl=en. > -- You received this message because you are subscribed to the Google Groups "sympy" 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/sympy?hl=en.
