clacula la integral mediante rectángulos

%integral
clear all;clc; close all;
num=1024;
xx=linspace(0,10,num);
fun=xx.^3+0.25*xx+3;
plot(xx,fun);
deltax=xx(2)-xx(1); %almacena el incremento
integral=0;
for k=1:(num-1);
integral=integral+fun(k)*deltax; % realiza la integral por la aproximación mediante barras
end;