Simpson's 1/3 rule on the 50g

[EngineerX]

Hello,

I would like to implement the following on the 50g in RPL but don't know how. Seems simple but I am new to the language. I have a feeling the Prime would be a better choice for me who has some experience with Matlab, but the thing is that I prefer the 50g for nostalgic reasons, although I never really learned more than the trivial.

I would like to be able to calculate the area under a curve using Simpson's 1/3 rule for a given set of data.

For exemple:

For x = {x0 x1 x2 x3 ... xN} where N is an even number.

Compute
A = h/3*(y0 + 4*y1 + 2*y2 + 4*y3 + 2*y4 + 4*y5 +...+ yN)

Very much appreciated.