- #1
Hepth
Gold Member
- 464
- 40
I have found that while these two should be the same, mathematica does not evaluate them equally.
I believe the second one is correct.
Also please try :
Notice its just a simple sum of terms, but the last two are zero while the first is not. As soon as you include the first variable in the integrand, it doesn't seem to work.
If you change the region of the u integration to -1..1, then its fine.
Does Mathematica have a problem with evaluating delta's on when the function sits on the boundary of integration? I assumed it would just use the definition of:
Any ideas? This is Mathematica 7 for Ubuntu.
Code:
Integrate[(a + b el + c el^2) DiracDelta[u], {u, 0, 1}, {el, e1, e2}]
Integrate[ Integrate[(a + b el + c el^2) DiracDelta[u], {el, e1, e2}], {u, 0, 1}]
I believe the second one is correct.
Also please try :
Code:
Integrate[(1 + u) DiracDelta[u], {u, 0, 1}, {el, e1, e2}]
Integrate[(el) DiracDelta[u], {u, 0, 1}, {el, e1, e2}]
Integrate[(1 + u + el) DiracDelta[u], {u, 0, 1}, {el, e1, e2}]
Notice its just a simple sum of terms, but the last two are zero while the first is not. As soon as you include the first variable in the integrand, it doesn't seem to work.
If you change the region of the u integration to -1..1, then its fine.
Does Mathematica have a problem with evaluating delta's on when the function sits on the boundary of integration? I assumed it would just use the definition of:
Code:
Integrate[ DiracDelta[u], {u, -1, 1}]
Integrate[ DiracDelta[u], {u, 0, 1}]
Any ideas? This is Mathematica 7 for Ubuntu.