Write M-function function S myquad (fun, a, b, n) that calculates a numerical value of the integral of a function using Simpson quadrature (Simpson 1/3 rule), where a and b are lower and upper limits and n is the number of subintervals. Since the number of subintervals has to be multiple of 2, you may generate a warning and then force it to be multiple of 2 by adding one more interval.