Question

Need help with MATLAB, what code should be added to print the graph of the trapezoid rule code below?

%%Matlab code for Question 5. syms X y intlx exp(x),0,2); %integration of x*exp(x) fprintf("htlntegration of x"2*exp(x) for [O 3] is %2.5f.\n,y) %Integration using Trapizoidal rule nn [100 1000]; for ii 1:2 a:0; b-3; %Integration limit ns-1:nn(i) integral Values-zeros(size(ns)); ourFunction-@(x) x.*2.*exp(x); for n-ns val(n)-trapizoidal (ourFunction,a,b,nn(i); end fprintf nlntegration of x 2*exp(x) using Trapizoidal rule for [O 3]\n') fprintf('1tTrapizoidal integration value for n-%d is 96f .\n.nn(ii),val(end)) end %%Matlab function for Trapizoidal integration function val-trapizoidal(func,a,b,N) % func is the function for integration 96 a is the lower limit of integration % b is the upper limit of integration % N number of rectangles to be used val-0 %splits interval a to b into N+1 subintervals xx linspace(a,b,N+1); dx-xx(2)-xx(1); %x interval %loop for Riemann integration for i-2:length(xx)-1 val-val+dx*double(func(xxl)); end val-val+dx*(0.5*double(func(xx(1))+0.5*double(func(xx(end))); end

Answer 1

`Hey,

Note: Brother in case of any queries, just comment in box I would be very happy to assist all your queries

clear all
clc
syms x
y=int(x^2*exp(x),0,3);
fprintf('\tIntegration of x^2*exp(x) is %2.5f\n',y);
nn=[100,1000];
for i=1:2
a=0;
b=3;
f=@(x) x^2*exp(x);
val=[];
for j=1:nn(i)
val(j)=trapezoid(f,a,b,j);
end
fprintf('Integration using trapezoidal for n=%d is %2.5f\n',nn(i),val(end));
figure;
plot(val);
xlabel('Number of intervals');
ylabel('Integral');
end
function I = trapezoid(func,a,b,n)

% I = trap(func,a,b,n):

% multiple-application trapezoidal rule.

% input:

% func = name of function to be integrated

% a, b = integration limits

% n = number of segments

% output:

% I = integral estimate

x = a;

h = (b - a)/n;

s = feval(func,a);

for i = 1 : n-1

x = x + h;

s = s + 2*feval(func,x);

end

s = s + feval(func,b);

I = (b - a) * s/(2*n);
end

MATLAB R2018a E C rch Documentation Log In Figure 1 ー Figure 2 File Edit View Insert Tools Desktop Window Help ile Edit View Insert Iools Desktop Window Help New Open 280 260 240 220 200 180 160 280 260 240 220 200 180 160 140 120 100 Current Folder mathlab in a1.m A1.dat A1.mat A2.dat A3.dat 120 100 80 10 20 30 40 50 6070 80 90 100 Number of intervals Name ▲ 0 100 200 300 400 500 600 700 800 900 1000 Number of intervals Integration using trapezoidal for n- Integration using trapezoidal tor n-1000 1s 98.42791 2 expx) 1000 100,1000] x1000 double 1x1 sym x1 sym trapezoid Ln 20 Col 4 ENG 913 AM yA ^ 4) IN 04/19/19 O Type here to search

Kindly revert for any queries

Thanks.