...
Code Block |
---|
language | none |
---|
linenumbers | true |
---|
|
% 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 . |
d)
Test funksjonen med simpsons(0,pi, 14,
@(x)sin(x)
på funksjonen ved følgende kall:
Code Block |
---|
f = @(t) sin(x)
simpsons(0,pi, 14, f) |
...
Code Block |
---|
|
n =200;
c = 1;
f = @(t) sin(x)
for i = [2:2:n]
e(c) = simpsons(0,pi,i,@sinf);
c = c+1;
end
loglog([1:n/2], e-2) |
...
Nyttig:
Du kan få bruk for atfølgende som ikke er pensum:
Definerer en vanlig (mattematisk) funksjon:
...