d=0;n=0
for k=1:2500
    tirage=rand()
    if tirage<1/3 then d=d+1;n=n+1
    else 
        if tirage<2/3 then
            if rand()<1/2 then n=n+1
            end
        end
    end
end
afficher(d/n)