Matlab function Simpson given below implements Simpson's 1/3 Rule to estimate

Question

Matlab function Simpson given below implements Simpson's 1/3 Rule to estimate

Image transcription text

Problem 5. Matlab function Simpson given below implements Simpson's 1/3 Rule to estimate I = So f(x)dx using n odd points. function I = Simpson (fun, a, b, n) % simpson Composite 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 integral n = number of nodes, odd and > 2. % Output : I = approximate value of the integral from a to b of f(x)*dx if (mod (n, 2) == 0) || (n < 3) disp(' Warning: n must be > 1 and odd, program terminates . .. ') ; return; end dx = (b-a) / (n-1) ; % stepsize x = a: dx: b; divide the interval f = feval (fun, x) ; % evaluate integrand I = (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 I 1 = 37 sin(x/3) da 2 - cos(2/2) (6) correct to at least 4 significant digits, by the use of Simpson's 1/3 rule. For this purpose start with n = 3 and increase n by 2 at a time until your result meets the required accuracy. Use Matlab function Simpson given above for the solution of this problem. Compare your answer with Matlab's quad () function. Comment on your computation.

Answer 1

llentesque dapibus efficitur laoreet. Nam risus ante, dapi

entesque dapibus efficitur laoreet. Nam risus ante, dapibus a molestie consequat, ultrices ac magna. Fusce dui lectus, congue vel laoreet ac,

l

sus ante, dapibus a molestie consequat, ultrices ac magna. Fusce dui lectus, congue vel laoreet ac, dictum

Answer 2

amet, consectetur adipiscing elit. Nam lacinia pulvinar tortor nec facilisis. P

Fusce dui lectus

e vel laoreet ac, dictum vitae odio. Donec aliquet. Lorem ipsum dolor sit amet, consectetur adipiscing elit. Nam lacinia pulvinar tortor nec facilisis. Pellentesque dapibus efficitur laoreet. Nam risus ante,