// FORMULA DI SIMPSON COMPOSTA PER IL CALCOLO APPROSSIMATO DI INTEGRALI
//
// function I = csimp(funz,a,b,n)
//
// DATI DI INPUT:
//  funz = stringa contenente il nome della funzione integranda
//         (definita in un omonimo file .m)
//     a = estremo inferiore dell'intervallo di integrazione
//     b = estremo superiore dell'intervallo di integrazione
//     n = numero di sottointervalli a cui applicare la formula di Simpson
//
// DATI DI OUTPUT:
//     I = approssimazione dell'integrale
//

function I = csimp(funz,a,b,n)

h = (b-a)/(n);
x = a + [0:n]*h;
I = (feval(x(1),funz) + feval(x(n+1),funz) + 2*sum(feval(x(2:n),funz)) + 4*sum(feval((x(1:n)+x(2:n+1))/2,funz))) * h/6;

//inf = feval(x(1),funz);
//sup = feval(x(n+1),funz);
//inte = 2*sum(feval(x(2:n),funz),'m');
//mid = 4*sum(feval((x(1:n)+x(2:n+1))/2,funz),'m');
//I = (inf + sup + inte + mid) * h/6;