integrales mediante rectángulos, cambiando el numero de ptos tomados

clear all;
close all;
clc;
%------------
NUM=[256, 512, 1024, 2048, 4096]; %diferentes numeros de ptos tomados
numPUN=length(NUM);
tic; % indica el tiempo que tarda en realizar el ejercicio
for m=1:numPUN
ndata=NUM(m);%número de puntos de xx
xx=linspace(0,10,ndata);
fun=xx.^3 + 0.25*xx +3;
integral=0;
paso=xx(2)-xx(1);
for k=1:ndata-1
integral=integral + fun(k)*paso;
end%fin del bucle con contador 'k'
resul(m)=integral;
end%fin del bucle con contador 'm'
timeE=toc;