Trapèzes et Simpson

hello tout le monde!! J'espère que vous allez tous bien!
g vraiment du mal à faire les methodes de simpson et des trapezes sans faire intervenir de fonctions....
Quelqu'un aurait-il des programmes qui fonctionnent???
merci bcp!

Code:

%% TRAPEZE/SIMPSON

clc;clear;
format long;
a=1;
b=4;
N=2000;
H=(b-a)/N;
x=a:H:b;
f=inline('x./(1+x.^4)');
y=f(x);

solTrapeze=(((f(1)+y(end))/2)+(sum(y(2:end-1))))*H
solSimpson=(y(1)+y(end)+4*sum(y(2:2:end-1))+2*sum(y(3:2:end-2)))*H/3

% Vérification :

quad(f,a,b)

Pour simpson tu dois juste remplacer la derniere ligne par:

I=h/3*((y(1)+y(end))+(4*sum(y(1:2:end-1)))+(2*sum(y(2:2:end-2))))

c'est juste une autre maniere de calculer