% a) Input
nn=10^6
aa=1
bb=2
f=@(x) x.^2

% b) Listen 
hh=(bb-aa)/nn;
ii=0:nn;
xi=aa+ii*hh;
yi=f(xi);

% c) Summieren (Trapezmethode)
A=hh*(sum(yi(2:end-1))+(yi(1)+yi(end))/2)
A=hh*(sum(yi(2:(nn)))+(yi(1)+yi(nn+1))/2)

% d) Untersumme/Obersumme (für monoton steigende Funktion)
Au=hh*sum(yi(1:(nn-1)))
Ao=hh*sum(yi(2:nn))