NumericalIntegration.jl icon indicating copy to clipboard operation
NumericalIntegration.jl copied to clipboard

initial value for cumulative integration

Open jmcclena opened this issue 4 years ago • 1 comments

When I use culum_integrate, it appears that it integrates 1 index beyond the array.

using NumericalIntegration
a = LinRange(0,1,100)
print(integrate(a,a), ' ')
print(last(cumul_integrate(a,a)))

returns:

0.5 0.5000510152025299

Is there a way to define the initial value of the integration so that it returns an integration for 0 to 1 with the first index being the initial value?

jmcclena avatar May 25 '21 22:05 jmcclena

Sorry, my description of the problem was wrong. It is not an indexing thing. Its just that the initial value is chosen defaulting to 5.1e-5 for cumul_integrate. For my purposes, I can just subtract off the index 1.

jmcclena avatar May 26 '21 21:05 jmcclena