Simpson's rule is a method for numerical integration.
intsim(f, delta)
f
Function values to be integrated.
delta
Spacing size.
intsim returns the value of the intergral.
intsim