function res=estisortie(x,a,b,m,T,h)
    X=procOU(m,h,T,x);
    sorta=0;
    sortb=0;
    for k=1:1:m
        z=X(k,:);
        testa=bool2s(z<=a);
        testb=bool2s(z>=b);
        test=sum(testa) + sum(testb);
        if (test>0) // on regarde si on sort de l'intervalle
            auxa=cumsum(testa);
            auxb=cumsum(testb);
            auxa=sum(bool2s((auxa<1)))+1; // le temps de sortie en a
            auxb=sum(bool2s((auxb<1)))+1; // le temps de sortie en b
            if (auxa<auxb) then
                sorta=sorta+1;
            else 
                sortb=sortb+1;
            end     
        end   
    end
    total=sorta+sortb;
    proba=primigauss(x,b)/primigauss(a,b);
    probb=primigauss(a,x)/primigauss(a,b);
    if (total>0)
        res=[sorta/total sortb/total ; proba probb]; // premiere ligne = valeur empirique;  deuxieme ligne = valeur empirique
    end
endfunction