clear all
clc

x=linspace(-1,1);
a=x(1);
b=x(end);

f=@(x) cos(x.^3);

integral(f,a,b)

w=ones(size(x))*(x(2)-x(1))/3;
for i=2:length(x)-1
    if rem(i,2)==0
        w(i)=4*w(i);
    else
        w(i)=2*w(i);
    end
end

F=f(x);

I=w*F'