Spero che quello che dici sia vero .... non si scherza sulle malattie.
Grazie a questa presunzione di innocenza ti fornisco alcune indicazioni.
Se segui questo link ....
http://mathworld.wolfram.com/Newton-CotesFormulas.htmlVedrai che le formule altro non sono che le classiche formule di
Newton - Cotes per cui i coefficienti sono quelli che trovi in letteratura. Una verifica era comunque d'obbligo per non avere sorprese a posteriori.
Per quanto riguarda il codice
Matlab basta guardare la relativa documentazione che trovi qui :
http://www.mathworks.it/it/help/matlab/ref/quad.htmlQuindi per fare le tue prove il codice per una funzione del tipo :

Ti scrivi una funzione e la salvi con lo stesso nome che hai dato all'interno del codice ( in questo caso myfun.m )
Poi da workspace digiti il comando istantaneo
E se tutto funziona bene la risposta dovrebbe essere il valore dell'integrale.
Questo deve essere confrontato con il valore ottenuto tramite un metodo di integrazione che implementi tu.
Se per esempio utilizzi la formula di
Cavalieri - Simpson ( nel tuo caso

e

) i coefficienti sono 1 4 1 per cui :

Devi quindi scrivere un semplice programmino che ti calcola la funzione in 3 punti ed esegue il calcolo impiegando anche lo step di integrazione

.