function [s]=integrale(a,b,n,f)
delta=(b-a)/n;
s=(f(a)+f(b))/2;
x=a;
for k=1:n-1
    x=x+delta;
    s=s+f(x);
end
s=s*delta;
end
