// Monografia de Analisis Numerico
// Proyecto FONDEF D05I-10211
// Autor: Axel Osses, 2008  
// Herramientas para la formacin de profesores
// Prohibida la reproduccin parcial o total de este software
//
// Formulas de cuadratura (Capitulo 3)
// Descripcion: Comparacion de las formulas
// de rectangulos, trapecios y Simpson

function func=f(x)
 func=sin(x);
endfunction;

//Intervalo
a=0
b=%pi
// Paso
N=10
h=(b-a)/N
//Puntos de integracion
x=[a:h:b];

// Cuadraturas
suma0=0;
suma1=0;
suma2=0;
suma3=0;

for i=1:N
 suma0=suma1+h*f(x(i));
 suma1=suma1+h*f(x(i)+h/2);
 suma2=suma2+h/2*(f(x(i))+f(x(i+1)));
 suma3=suma3+h/6*(f(x(i))+4*f(x(i)+h/2)+f(x(i+1)));
end;

integral_rectangulos1=suma0
integral_rectangulos2=suma1
integral_trapecios=suma2
integral_Simpson=suma3



