Численное интегрирование функций

clc
clear
a=xlsread('lab6',1,'B1'); %нижняя граница
b=xlsread('lab6',1,'C1'); %верхняя
n=10;  
h=(b-a)/n;%шаг
x=a:h:b;
y= @(x) (4*x^2+1.9*x).^1/3;
m4 = max(diff(y(x),4));
eps=0.0001; %точность
f=1;
while f==1
    if(h <= (180*eps)/((b-a)*m4)^(1/4))
        f=0;
    else
        n = n-1;
        h=(b-a)/n;
    end
end
disp('Резульат, который мы должны получить:')
proverka = trapz(x,y(x))
temp = y(x);
p=temp(length(temp));%значение последнего элемента функции у
f=1;
s=[0 0]; 
%реализация метода
while f==1
    for i=1:2
        for x=(a+h):(2*h):(b-h)%нечетные числа
            s(i)=s(i)+(4.*y(x));
        end
        for x=(a+2*h):(2*h):(b-(2*h))%четные
            s(i)=s(i)+(2.*y(x));
        end
        s(i)=(s(i)+p)*h/3;
        n=2*n;
        h=(b-a)/n;
    end
%если проходит проверку на точность то выход из цикла
if abs(s(1)-s(2))<(eps/2)
    f=0;
else
    n=n/2;
    h=(b-a)/n;
end
end
disp('Резульат, который мы получили:')
rezultat = s(i)
%проверка с помощью встроенной функции
xlswrite('lab6',rezultat,1,'A1');
xlswrite('lab6',proverka,1,'A2');