%numericka integrace

 f=inline('1./((x-.3).^2+.01)+1./((x-.9).^2+.04)-6');   %inline funkce
 fplot(f,[0 1]);    %graf
 hold on
 stem(0:.02:1,f(0:.02:1));  %stopkový graf
 I1=quad(f,0,1) %urèitý integrál numericky pomocí adaptivního Simsponova pravidla
 I2=quadl(f,0,1)    %pomocí adaptivní Gaussovy-Lobattovy kvadratury
 I3=quadv(f,0,1)    %vektorizované rekurzivní Simsponovo pravidlo
 
 f2=@(x,y) y*sin(x)+x*cos(y);   %anonymní funkce
  D1=dblquad(f2,pi,2*pi,0,pi)    %dvojný integrál, kde jsou meze pi <= x <= 2*pi a 0 <= y <= pi
 D2=dblquad(@integrand,pi,2*pi,0,pi, 'tol',@quadv)    %stejný integrál zadaný pomocí handlu funkce a metody @quadv
 
 f3=@(x,y,z) y*sin(x)+z*cos(y);
 T1=triplequad(f3,0,pi,0,1,-1,1)
 T2=triplequad(f3,0,pi,0,1,-1,1, 'tol',@quadl)
 