clear

f=@(z)integral((250*z)/(6+z)*exp(-z/10),0,30);
vastusSimpson=quadv(f,0,30)%vaikimisi täpsus e-6

%trapsetsvalemiga
%teeme kaks vektorit
x=0:1e-3:30;
m=f(x);
vastusTrapets=trapz(x,m)

%õige vastus: 1,15e3