function [Prob,Temps]=Ivrogne(a,b)
    S=0;
    T2=0;
    for k=1:2000
    X=0;
    T=0;
    while((X>(-a))& (X<b))
        X=X+(2*(rand(1)<0.5)-1);
        T=T+1;
    end;
    T2=T2+T;
    S=S+(X==(-a));
    Prob=S/2000;
    Temps=T2/2000;
    end
endfunction

Ivrogne(3,4)

a=3; 
Pr=zeros(1,10);
Temps=zeros(1,10);
for b=1:10
    [Pr(b),Temps(b)]=Ivrogne(a,b);
end
b=[1:10];
plot2d(b,Temps)
figure
plot2d(b,Pr);
plot2d(b,b./(a+b),2)