Problem 5. Matlab function Simpson given below implements Simpson’s 1/3 Rule to estimateI = So f(x)dx using n odd points.function I = Simpson (fun, a, b, n)% simpsonComposite Simpson’s rule%%% Synopsis: I = simpson (fun, a, b, npanel)% Input :fun= (string) name of the function that evaluates f (x)%a, b= lower and upper limits of the integraln= number of nodes, odd and > 2.% Output :I = approximate value of the integral from a to b of f(x)*dxif (mod (n, 2) == 0) || (n < 3)disp(' Warning: n must be > 1 and odd, program terminates . .. ‘) ;return;enddx = (b-a) / (n-1) ;% stepsizex = a: dx: b;divide the intervalf = feval (fun, x) ;% evaluate integrandI = (dx/3) *(f(1) + 4*sum(f (2:2:n-1) ) + 2*sum(f (3:2:n-2) ) + f(n));Now, it is desired to estimate the definite integral I1 =37 sin(x/3) da2 – cos(2/2)(6)correct to at least 4 significant digits, by the use of Simpson’s 1/3 rule. For this purpose start withn = 3 and increase n by 2 at a time until your result meets the required accuracy. Use Matlab functionSimpson given above for the solution of this problem. Compare your answer with Matlab’s quad ()function. Comment on your computation.

Don't use plagiarized sources. Get Your Custom Essay on
Matlab function Simpson given below implements Simpson’s 1/3 Rule to estimate
Just from $10/Page
Order Essay

Order your essay today and save 30% with the discount code ESSAYSHELP