Matlab/Octave/freemat Implementation of Simpsons 1/3 rule for numerical integration