% x = simpsonQuadrature(fun,tLow,tUpp,n)
% This function uses simpson quadrature over each of n uniform segments to
% approximate the integral of fun(t) on the interval tLow <= t <= tUpp
% fun = function handle
% f = fun(t)
% t = [1, nt] = time query points on [a,b]
% f = [nx, nt] = vector function at query points
% tLow = scalar lower bound on time
% tUpp = scalar upper bound on time
% n = number of segments to divide interval into
% x = [nx, 1] = integral along each dimension
Matthew Kelly (2019). simpsonQuadrature (https://www.mathworks.com/matlabcentral/fileexchange/54457-simpsonquadrature), MATLAB Central File Exchange. Retrieved .
Thanks for the heads-up. I've added that file (and one other) to my acknowledgements, since they were posted first.
added reference to existing simpson quadrature routine on file exchange.