function ApprocheMorris(Fonction,Parametres,Entrees,EspacePhases,Nd,Nr)

    for i=1:size(EspacePhases,1)
    
        Grille(i,:)=[EspacePhases(i,1)+(EspacePhases(i,2)-EspacePhases(i,1))/(2*Nd):(EspacePhases(i,2)-EspacePhases(i,1))/(Nd):EspacePhases(i,2)-(EspacePhases(i,2)-EspacePhases(i,1))/(2*Nd)];
    
    end
    
    for t=1:Nr
        
       IndicesEntrees=randi([1 size(Grille,2)],length(Entrees),1);
       
       EntreesTestees=[];
       
       for i=1:length(Entrees)
           
           EntreesTestees=[EntreesTestees,Grille(i,IndicesEntrees(i))];
           
       end
       
       ValInf=feval(Fonction,Parametres,EntreesTestees);
       
       EntreesDispos=[1:1:length(Entrees)];
       
       for e=1:length(Entrees)
          
           Indice=randi([1 length(Entrees)+1-e],1,1);
           
           EntreeInteret=EntreesDispos(Indice);
           
           EntreesDispos(Indice)=[];
           
           Variation=randi([1 2],1,1);
           
           IndicesEntreesNew=IndicesEntrees;
           
           NewIndice=IndicesEntreesNew(EntreeInteret)+(2*Variation-3);
           
           if (NewIndice>Nd) 
               
               NewIndice=Nd-1;
               
           elseif (NewIndice==0)
           
               NewIndice=2;
               
           end
           
           IndicesEntreesNew(EntreeInteret)=NewIndice;
           
           
            EntreesTesteesNew=[];
       
            for i=1:length(Entrees)
                
                EntreesTesteesNew=[EntreesTesteesNew,Grille(i,IndicesEntreesNew(i))];
                
            end
       
            ValSup=feval(Fonction,Parametres,EntreesTesteesNew);
       
            E(t,EntreeInteret)=abs(ValSup-ValInf)*Nd/(EspacePhases(EntreeInteret,2)-EspacePhases(EntreeInteret,1))*EntreesTestees(EntreeInteret)/ValInf;
       
            IndicesEntrees=IndicesEntreesNew;
            
            EntreesTestees=EntreesTesteesNew;
           
       end
        
    end
    
    MoyE=mean(E);
    
    EType=std(E);
    
    MAXI=max(max(MoyE),max(EType));
    
    figure; plot(MoyE,EType,'+',[0,MAXI],[0,MAXI]); grid on; xlabel('\mu'); ylabel('\sigma'); title('Approche de Morris')
    
    hold on ;
    
    for i=1:length(Entrees)
        
        text(MoyE(i)+MAXI/40,EType(i),num2str(i));
        
    end

end