I'm trying to integrate using Simpson's rule

I'm trying to integrate sin(x) from 0 to pi using the simpson's rule. The only thing is I need it to output several iterations each with a different number of intervals. For example, n= 4, 8 , 16, 32, 64. Here is my code so far clear; close all; clc;

 ezplot('sin(x)', [0, pi]), hold on

 f=@(x) sin(x); 
a=0;b=pi; 
n=4; 
h=(b-a)/n;
p=0;

 for i=a:b
p=p+1;
x(p)=i;
y(p)= sin(i); 
plot(x, sin(x), 'r')
end

 l=length(x);
x
y

any help appreciated!

Answer 1

[Walter Roberson]

So instead of "n=4" use

for n = [4, 8, 16, 32, 64]

and put an "end" after.

[ray sanchez]

Where would I put the "end"? I tried it in several places and it didn't work

[Walter Roberson]

for n = [4, 8, 16, 32, 64]
  h=(b-a)/n;
  p=0;

    for i=a:b
      p=p+1;
      x(p)=i;
      y(p)= sin(i); 
      plot(x, sin(x), 'r')
      drawnow();
    end

     l=length(x);
    x
    y
  end

Right at the end.

Note that I also added a drawnow() call so you have a hope of seeing the plots you are generating. I have no idea why you would want to plot in the manner you are plotting, but I guess you must have a reason.