Do the same thing as was done in problem 1 with Simpson's rule. sol) function I=simp(n) x0=0;, xn=pi; ... *(f(1)+4*sum(f([2:2:2*n]))+2*sum(f([3:2:2*n-1]))+f(2*n+1)); >> simp(4) ans = 1.*************9 >> simp ... (8) ans = 1.29278827445628 >> simp(16) ans = 1.3************* >> simp(32) ans = 1.30235988785545 >> simp