> 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.
