function int = simpsonc(a,b,n,fv)

% A program az összetett Simpson-formulával közelíti az fv függvény
% integrálját az [a,b] intervallumon úgy, hogy az intervallumot 
% n egyenlő részre osztjuk.

h=(b-a)/n; x=[a:h/2:b]; 
y=eval(fv);
int=(h/6)*(y(1)+2*sum(y(3:2:2*n-1))+4*sum(y(2:2:2*n))+y(2*n+1));