Din uppgift är att skriva ett program i Python som implementerar trapetsregeln. Om du vill ha hjälp så finns pseudokod för programmet i länken nedan.
Inled filen med import math för att kunna använda matematikfunktionerna i Python.
Lägg sedan till programkoden nedan.
def minfunktion(x):
return math.exp(-x**2)
Du har nu definierat funktionen minfunktion som \(f(x)=e^{-x^2}\). Vill du beräkna \(e^{-3^2}\) kan du alltså skriva minfunktion(3)
i programmet. Du slipper då skriva in hela ditt krångliga funktionsuttryck varje gång. När du ska byta funktion byter du math.exp(-x**2)
till lämpligt uttryck.
I programmet ska du kunna variera funktionen, integrationsgränserna och antalet delintervall.
Beräkna integralen \(\int_0^1e^{-x^2}dx\) med ditt program. Testa först att du får samma resultat som i uppgift 1 med fyra delintervall. Bestäm sedan integralen med minst 5 korrekta decimaler och beräkna det procentuella felet när endast 4 delintervall används.