function logerreur(A,B,fun,Nmin,Nmax,N,Iexa)

% trace le diagramme log log et affiche les 4 pentes.
% on intègre la fonction f sur l'intervalle [A,B] en prenant N points
% dans l'intervalle fix(logspace(Nmin,Nmax,N)
% Iexa est la valeur exacte de l'intégrale.
% cf Exemple 3.27 du poly de cours
% Attention, utilise la fonction int_fcn, à remplacer éventuellement par
% d'autres fonctions d'intégration.
% Exemple :
% logerreur(0,pi/2,'sin',1,4,100,1)

er=zeros(4,N);
Ntot=fix(logspace(Nmin,Nmax,N));
H=(B-A)./Ntot;
for i=1:N
    for j=1:4
        er(j,i)=int_fcn(j,Ntot(i),A,B,fun);
    end
end
er=abs(er-Iexa);
loglog(H,er);
legend('rectangles','trapèzes','point milieu','Simpson','location','best');
for j=1:4
    ind=er(j,:)~=0&er(j,:)>=1e-13;
    p=polyfit(log10(H(ind)),log10(er(j,ind)),1);
    disp(['pente méthode',int2str(j)]);
    disp(p(1));
end