(a) Write a Matlab function that compute the numerical integral using the composite Boole's rule with the definition function num1 = compBoole(f, a, b, N) where f function to be integrated left and right end points of the interval, a, b : N : number of intervals, 1 integral. (b) Use your function compBoole.m to compute the integral of f(x) = e-x2/2 on the interval [-3,3] with N uniform intervals. Take number of intervals N to be 5 * 2, with k = 1, 2, 3, 4, 5, The exact integral is V2? erf(3V2/2), where erf is a matlab built-in function. Print out the number of divided intervals along with the absolute numerical error. Plot the error against the divided interval length (b - a)/N 6/N using loglog scale. Numerically compute the rate of convergence for this rule. You shall observe the rate is approximately 6