%
%
%       Biofisica 2005 - FCEyN - UBA
%       Random Walks canales
%       Jacobo Sitt
%

function w=rw_canal(k1,k1p,k2,k2p,dt,n)

% 0 = c1 , 1 = c2 , 2 = c3  (codificacion de estados)

w(1)=0;  % Lugar de comienzo del random walk


for i=1:n
    r=rand ;
    
    if w(i)==0  
        if r<k1*dt
            w(i+1)=1;  % paso de c1 a c2 con prob k1*dt
        else
            w(i+1)=0;  % el canal se queda en c1
        end
    end
    
    if w(i)==1  
        if r<k1p*dt
            w(i+1)=0;  % paso de c2 a c1 con prob k1p*dt
        end
        if r>k1p*dt & r<(k1p+k2)*dt
            w(i+1)=2;  % paso de c2 a o con prob k2*dt
        end
        if r>(k1p+k2)*dt
            w(i+1)=1;  % el canal se queda en c2
        end        
    end
    
    if w(i)==2  
        if r<k2p*dt
            w(i+1)=1;  % paso de abierto a c2 con prob k1*dt
        else
            w(i+1)=2;  % el canal se queda abierto
        end
    end
    
end