function xn=simul_1(n)
    xn=1;
    for k=2:n
        xn=grand(1,1,"uin",1,3)
    end
endfunction

function pk=loi_1(m)
    sk=[0 0 0];
    for i=1:m
        test=simul_1(50)
        sk(test)=sk(test)+1
    end
    pk=sk/m
endfunction

function xn=simul_2(n)
    xn=1;
    for k=2:n
        if xn==1 then
            if rand()<.5 then 
                xn=2 
                else xn=3
            end
        elseif xn==2 then 
            if rand()<.5 then
                xn=1 
                else xn=3
            end
        elseif rand()<.5 then 
                xn=1 
                else xn=2
             end
    end
endfunction

function pk=loi_2(m)
    sk=[0 0 0];
    for i=1:m
        test=simul_2(50)
        sk(test)=sk(test)+1
    end
    pk=sk/m
endfunction