Displaying 1 result from an estimated 1 matches for "shiftedgauss".
2010 Sep 21
1
puzzle with integrate over infinite range
...#39;s usually better to
specify Inf explicitly as a limit rather than an arbitrary large
number, as in this case integrate() performs a trick to do the
integration better.
However, I do not understand the following, if I shift the Gauss
function by some amount the integral should not be affected,
shiftedGauss <- function(x0=500){
integrate(function(x) exp(-(x-x0)^2/100^2), 0, Inf)$value
}
shift <- seq(500, 800, by=10)
plot(shift, sapply(shift, shiftedGauss))
Suddenly, just after 700, the value of the integral drops to nearly 0
when it should be constant all the way. Any clue as to what's go...