Schrijf een functie simpson() met als argumenten:

• f: de te integreren functie
• a: ondergrens van het integratie-interval
• b: bovengrens van het integratie-interval
• N: aantal subintervallen

De functie f wordt verondersteld 1 positioneel argument te hebben, dat zowel scalair als een NumPy-rij kan zijn. Je mag aannemen dat de functie gevectoriseerd is en dat de parameter N even is.

Het resultaat van de functie is de gezochte integraalwaarde, gebruik makende van de regel van Simpson. Beperk ook hier het aantal oproepen naar de functie f. Probeer je functie dus efficiënt te implementeren (vermijd m.a.w. onnodige herhalingen van berekingen van functiewaarden).

NOOT: het verbeterscript maakt gebruik van de methode veelterm(). Deze hoef je NIET mee in te dienen.

Voorbeeld

v = veelterm([1, 2, 3])
simpson(v, 1, 5, 100) #152.0000