...
Feilen for Simpsons metode er gitt som en funksjon av den fjerdederiverte til funksjonen. Ofte er det vanskelig å finne gode skranker for denne, men vi ønsker likevel å ha en viss kontroll på hvor stor feil vi gjør. Vi skal nå se på en mulig praktisk løsning på problemet.
...
Skriv funksjonen simpsons_error(start, stop, error, fn) som tar inn start og stop som integrasjonsgrenser, feiltoleransen error og funksjonen fn som skal integreres.
Regn ut integralet med en feiltoleranse på 10−8.Test funksjonen med følgende kall:
| Code Block | ||
|---|---|---|
| ||
f = @(t) sin(t) simpsons_error(0, 1, 10^-4, f) % skal returnere 0.459707744927311 % Skriv format long før du kjører koden. |
...