function I=integrale(x,N)
% Calcola l'integrale tra 1 e x (x>0) della funzione
%   f(x) = cos(log(x))
% usando la formula dei punti sinistri, con una partizione di N
% intervalli.

% Esempio:
% >> integrale(3,2000)
% 1.5184
    
if x<=0 | N<=0
    'dati iniziali non accettabili'
    return
end
h=(x-1)/N;      % Lunghezza di ogni itervallo della partizione.
r=1:h:(x-h);    % Punti della partizione.
                % Nota: il secondo estremo e' "x-h" non "x".
I=h*sum(cos(log(r)));