...
Noen funksjoner er transcendentale og kan derfor ikke integreres analytisk. Men de kan (av og til) integreres nummerisk, altså tilnærmes.
Nedenfor er Simpsons metode, som dere kanskje har lært om i matematikk 1, forklart kort.
Hvor
Her er for en i = [0,1,2,...,n],og funksjonen du vil integrere, f(x). n må være et oddetall.
a)
Fordi n må være partallodedtall, må funksjonen isEvenisOdd(n)
lages. Returvariabelen er 1 eller 0 om n er partall oddetall eller ikke.
Digresjon: Er det egentlig nødvendig å lage en egen funksjon for dette?
...
Code Block | ||||
---|---|---|---|---|
| ||||
% Ta inn a, b, n og funksjonen fn. % Sjekk at n er partalloddetall, gi feilmelding hvis ikke. % Gå videre: % La variabelen y inneha alle funksjonsverdiene til fn fra a til b med steglende deltaX. % velg odde yer: y1, y3, y5 ... y_(n-1) % regn ut 4*summen av disse % velg partall yer: y2, y4, y6 ... y_(n-2) % regn ut 2*summen av disse % legg sammen de to summene sammen med y_0 og y_n % gang alt med \deltaX/3 . |
...