clc;
clear all;
close all;

a=0; b=10; 
x=a:b;
f= @(x)x.^5; %función cuadrática
I_simp=quad(f,a,b);

I_a= (1/6)*(b^6-a^6); % método analítico

% comparación entre el método analítico y el de Simpson
fprintf('Valor analítico= %f \n Simpson: %f \n',I_a,I_simp)