Status: New
Owner: ----
Labels: Type-Defect Priority-Medium

New issue 4083 by [email protected]: mellin_transform(besselj(3, x)/x**3, x, s) wring result
http://code.google.com/p/sympy/issues/detail?id=4083

mellin_transform(besselj(3, x)/x**3, x, s)

Wrong value returned:
-2**(s - 4)*gamma(s/2 - 3)/gamma(-s/2 + 1)

Correct result is:
-2**(s - 4)*gamma(s/2)/gamma(-s/2 + 4)


In [1]: from sympy import *

In [2]: x, s = symbols('x s')
   ...: F, _, _  = mellin_transform(besselj(3, x)/x**3, x, s)
   ...: assert F == -2**(s - 4)*gamma(s/2)/gamma(-s/2 + 4)
   ...:
---------------------------------------------------------------------------
AssertionError                            Traceback (most recent call last)
<ipython-input-2-13c7c0886cab> in <module>()
      1 x, s = symbols('x s')
      2 F, _, _  = mellin_transform(besselj(3, x)/x**3, x, s)
----> 3 assert F == -2**(s - 4)*gamma(s/2)/gamma(-s/2 + 4)
      4

AssertionError:

In [3]: F
Out[3]: -2**(s - 4)*gamma(s/2 - 3)/gamma(-s/2 + 1)

--
You received this message because this project is configured to send all issue notifications to this address.
You may adjust your notification preferences at:
https://code.google.com/hosting/settings

--
You received this message because you are subscribed to the Google Groups 
"sympy-issues" 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/sympy-issues.
For more options, visit https://groups.google.com/groups/opt_out.

Reply via email to