clc;
clear all;
close all;
disp('Método de Simpson')
% los valores del intervalo [a, b]  son proporcionados
% por el usuario desde el teclado
a=input('Introduzca el valor inicial del intervalo de integración: ');
b=input('Introduzca el valor final del intervalo de integración: ' );
if a>=0 && b>=0 
    if a>b
        disp('Error: el intervalo positivo [a,b] debe cumplir a<b.')
        return
    end
% regla de Simpson
% la función raíz cuadrada (sqrt) es indicada como
% cadena de caracteres en la función quad
I_simp=quad('sqrt',a,b);
% método analítico
I_a=2/3*(b^(3/2)-a^(3/2));
fprintf('Valor analítico= %f \n  Simpson: %f \n',I_a,I_simp)
else % Para el caso de valores negativos
    disp('Error en los valores del intervalo de integración.')
    disp('El intervalo [a,b] con a<b debe contener únicamente valores positivos.')
end