a)
Code Block | ||
---|---|---|
| ||
function result = simpsons(a,b,n,f) if mod(n, 2) ~= 0 error('Parameteren n må være et partall, var %d\n', n); end dx = (b-a)/n; y= f(a:dx:b); yodd = 2*sum(y(3:2:n-1)); yeven = 4*sum(y(2:2:n)); ysum = y(1) + yodd + yeven + y(n+1); result = dx/3 * ysum; end |
b)
Code Block | ||
---|---|---|
| ||
function result = simpsons_error(start, stop, error, fn) N = 2; % N må være et partall Si = simpsons(start, stop, N, fn); S2i = simpsons(start, stop, 2*N, fn); while abs(Si-S2i) >= error N = 2 * N; Si = S2i; S2i = simpsons(start, stop, 2*N, fn); end fprintf('Antal ledd: %d\n', N); result = Si; end |