// Ma fonction integration
function [y,n]=integration(f,minx,maxx,dy,n1)
n2=n1+2;
dx=(maxx-minx)/n1;
x=minx:dx:maxx;

fx=eval(f);
y1=(sum(fx)-0.5*(maxx-minx))*dx;
d=2*dy;
plot(x,fx,'+-r');
while (d>dy);
// Calcul de la seconde intégrale  
    dx=(maxx-minx)/n2;
    x=minx:dx:maxx ;
    fx=eval(f);
    y=(sum(fx)-0.5*(fx(1)+fx(length(fx))))*dx;
// comparaison des résultats    
    d=abs(y1-y);
    y1=y;
    n2=n2*2;
// fin de la boucle    
end
n=n2;
plot(x,fx,'--g');
xlabel('Longueur (%lambda)');
ylabel('Energie %epsilon (  10^6 ergs) ');
ST1 =['Fonction sur '+mtlb_num2str(n1)+' points'] // ici on utilise la fonction matlab num2str
ST2 =['Fonction sur '+string(n2)+' points']       // ici on utilie la fonction scilab string (meme resultat)
legend([ST1;ST2]);


