How to plot graph of Simpson's rule and Trapezoidal rule between error and n?

Hi there,

I have written 2 programs of Simpson's rule and Trapezoidal rule to compute the integral of a function f. Now I would like to plot graphs of both rules between the errors and n but I don't know how to do it. Could anyone help me with this problem please?

Here is the code of my program:

a=0.1;
b=0.2;
n=30;
x=linspace(a,b,n);
h=(b-a)/n;
f=inline('x.^4-x.*5+1');
y=x.^4-x.*5+1;

H=h*ones(1,n); %trapezoid
H([1,n])=h/2;
tr=sum (H.*f(x));
tr

H(1:2:n)=2*h/3; %simpson
H(2:2:n)=4*h/3;
H([1,n])=h/3;
simp=sum(H.*f(x));
simp

Thank you in advance!

Answer 1

[Alan Stevens]

If you use cumsum instead of just sum then you can plot tr and simp against x.

Answer 2

[skerdi hymeraj]

where is ploting