Contents

Vänster rektangelregel i Matlab

f = @(x)x.*cos(x);
a = 0; b = 1; n = 3;
x = linspace(a,b,n+1);
h = (b-a)/n;
qv = sum(h*f(x(1:n)))
qv =
   0.279636829763181

Exakt svar

sin(1)+cos(1)-1
ans =
   0.381773290676036

Funktionen min_integral (uppgift 2 i lab6)

qv = min_integral(@(x)x.*cos(x),[0,1],3,1);
exakt = sin(1)+cos(1)-1;

Uppgift 4 i lab6

f=@(x)x.*cos(x);
q50=min_integral(f,[0,1],50,1);
q100=min_integral(f,[0,1],100,1);
exakt = sin(1)+cos(1)-1;
(q50-exakt)/(q100-exakt)
ans =
   2.007995479710493