NumericalIntegration.jl
NumericalIntegration.jl copied to clipboard
initial value for cumulative integration
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?
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.