Write a Matlab program/function myInt(a, b, n). It will calculate the integral of a user defined function f(x) (Stored in a .m file called f) over the interval [a, b] using both Trapezoidal and Simpsons Rule with n partitions. Round output answers to 8 decimals.

Example output

>>myInt(0,1,2)

Trapezoid: 0.70833333

Simpson: 0.69444444

Do not make any assumptions on the input n. If n is odd, the Simpsons output should be Error: n is not even