...
Code Block | ||
---|---|---|
| ||
% Ta inn a, b, n og funksjonen fn. % Sjekk at n er partall, 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 . |
#anchor
d)
Test funksjonen med simpsons(0,pi, 14, @sin).
...
Lag funksjonen deviation(a,b,n, fn, int). int er det analytiske svaret
...
svaret.
f)
Mattematisk bonusspørsmål:
Simpsons metode er en annenordens metode. Hva har det med det logaritmiske plottet under å gjøre?
Code Block | ||
---|---|---|
| ||
n =200; c = 1; for i = [2:2:n] |
...
e(c) = simpsons(0,pi,i,@sin); |
...
c = c+1; |
...
end |
...
loglog([1:n/2], e-2) |
...
...
Nyttig:
...
Du får bruk for at:
...
Definerer en vanlig (mattematisk) funksjon:
Code Block | ||
---|---|---|
| ||
function y = f(x) y = x; end |
...