please use matlab to do this problem. Thank you.

Write a program to evaluate I = integral^b_a f(x)dx using the composite Simpson's rule with n subdivisions, calling the result I_n. Use the program to calculate the following integrals with n = 2^k, k = 1, ...9. (a) integral^1_0 x^5/2 dx (b) integral^pi_0 sin(sin(x))dx Analyze empirically the rate of convergence of I_n to I by calculating R_n = I_2n - I_n/I_4n - I_2n, Discuss your findings