function y = simpson(F,a,b,n)

//
// Integracion usando la regla de Simpson
//

y = 0;
x = a;
h = (b-a)/n;


impar = n - int(n/2)*2;   // Verificando si se tiene un numero impar de intervalos


for i = 1:(n-impar)/2
   y = y + (h/3)*(feval(x,F)+4*feval(x+h,F)+feval(x+2*h,F));
   x = x + 2*h;
end
if impar ~= 0
   y = y + (h/2)*(feval(x,F)+feval(x+h,F));
end

endfunction
