Comment #9 on issue 2624 by asmeurer: Sympy 0.7.1 can't integrate Gaussians
http://code.google.com/p/sympy/issues/detail?id=2624
By the way, that would indeed be very cool if we could do the integral from
that paper. The author wants to know how Mathematica does the integral.
But with Tom's code, you just have to set SYMPY_DEBUG to True and it will
give you detailed output of how it is working:
In [7]: integrate(exp(-(x - 1)**2), (x, -oo, oo))
Integrating exp(-(x - 1)**2) wrt x from -oo to oo.
Integrating -oo to +oo.
Sensible splitting points: set([1])
Trying to split at 1
Trying originial integrand
Integrating exp(-_x**2)
Could rewrite as single G function: 1 1 [(1, 0, meijerg(((), ()), ((0,),
()), _x**2))]
Checking antecedents for 1 function:
delta=1/2, eta=exp_polar(0), m=1, n=0, p=0, q=1
ap = [], []
bq = [-1/2], []
cond_3=True, cond_3*=True, cond_4=True
case 1: [False, False, False]
case 2: [True]
case 3: [False, False]
extra case: [False]
second extra case: [True]
Trying originial integrand
Integrating exp(-_x**2)
Could rewrite as single G function: 1 1 [(1, 0, meijerg(((), ()), ((0,),
()), _x**2))]
Checking antecedents for 1 function:
delta=1/2, eta=exp_polar(0), m=1, n=0, p=0, q=1
ap = [], []
bq = [-1/2], []
cond_3=True, cond_3*=True, cond_4=True
case 1: [False, False, False]
case 2: [True]
case 3: [False, False]
extra case: [False]
second extra case: [True]
Out[7]:
⎽⎽⎽
╲╱ π
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" 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-issues?hl=en.
