function trapecioc=trapecioc(a,b,n)
s=f(a);
h=(b-a)/n;
for i=1:n-1
    s=2*f(a+i*h)+s;
end
s=s+f(b);
trapecioc=s*h/2;