Do the same thing as was done in problem 1 with Simpson's rule.sol)function I=simp(n)x0=0;, xn=pi;h=( ... *(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.08301683149489>> simp(8)ans ... =1.29278827445628>> simp(16)ans =1.30184166534723>> simp(32)ans =1.30235988785545>> simp(64)ans =1.302391582830424