In deze opdracht bepalen we een primitieve functie langs numerieke weg, gebruik makend van de standaard trapeziumregel. Schrijf een functie numerieke_primitieve(), met als argumenten:

• f: de functie waarvan de primitieve $$F$$ bepaald moet worden
• c: de keuze van $$F$$ waarvoor $$F(c) = 0$$
• N: het aantal subintervallen

Het resultaat van deze functie is een nieuwe, gevectoriseerde functie, die voor het argument $$x$$ de waarde $$\int_c^x f(u)du$$ oplevert. Je mag aannemen dat het argument een NumPy-rij is. Probeer dit te programmeren zonder gebruik te maken van for-lussen.

Voorbeeld

p = numerieke_primitieve(lambda x:3*x**2 + 2*x**1 + 1, 1.0, 10)
r = p(np.array([-2.0, 0.0, 2.0])) #[-9.1350, -3.0050, 11.0050]
p = numerieke_primitieve(lambda x:3*np.sin(3*x) + 2*np.sin(2*x) + 1*np.sin(1*x), 1.0, 10)
r = p(np.array([-2.0, 0.0, 2.0])) #[-0.6371, -3.8458, -0.7432]