fun=@（x）exp（x）;
a=0;
b=1;
epsilon=10^-6;

n=1;
h=（b-a）/2;
y0=h*（feval（funa）+feval（funb））;
yiter=y0;
while?1
????????f=sum（feval（funa+（1:2:2*n-1）*h））;
????????y=y0/2+h*f;
???