powered by
Internal function to calculate integration weights for Simpson's rule
simpson(dim)
a double vector for the integration weights
integer for the number of nodes