Question: Bug using Heaviside and PDF in statistics package

Integrating a positive definite function (normal distribution) and a >= 0 function (Heaviside) should not return a negative value.

 

 

with(Statistics);
X := RandomVariable(Normal(1, sqrt(2.25)));
int(PDF(X, x)*Heaviside(x^7-5*x^4-3*x+1), x = -infinity .. infinity);
                         -0.08507120131

 

 

Bug.mw

Please Wait...