Updates:
        Status: Started
        Owner: julien.r...@gmail.com
        Labels: NeedsReview

Comment #8 on issue 2630 by julien.r...@gmail.com: DeltaIntegrate provides wrong answer on Multivariable integrals
http://code.google.com/p/sympy/issues/detail?id=2630

This works in https://github.com/sympy/sympy/pull/1620:

In [1]: densityX = exp(-x**2/2)/sqrt(2*pi)

In [2]: densityY = exp(-y**2/2)/sqrt(2*pi)

In [3]: integrate(densityX*densityY, (x, -oo, oo), (y, -oo, oo))
Out[3]: 1

In [4]: integrate(densityX*densityY * DiracDelta(10*x+y-1), (x, -oo, oo), (y, -oo, oo)).evalf()
Out[4]: 0.0395002101563691

In [5]: integrate(densityX*densityY * DiracDelta(x+10*y-1), (x, -oo, oo), (y, -oo, oo)).evalf()
Out[5]: 0.0395002101563691

The example from comment #8 was added as a test.

--
You received this message because you are subscribed to the Google Groups 
"sympy-patches" group.
To post to this group, send email to sympy-patches@googlegroups.com.
To unsubscribe from this group, send email to 
sympy-patches+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/groups/opt_out.


Reply via email to