
%Nummerisk integrasjon med trapes metoden. 
%Intervallet funksjonen f integreres over er [a,b]
%Intervallet deles i n like brede delintervaller, med
%bredde lik d. 

a=0;
b=1;
f=@(x) x*2 ;

for k=0:4
n=10^k; 

d=(b-a)/n;
S=0;
for i=1:n
 S=S+ d*(f(a+ (i-1)*d) + f(a+i*d)) ;
end
S/2
end