%
%
%       Biofisica 2005 - FCEyN - UBA
%       Simulacion Canales
%       Jacobo Sitt
%


function canales


    ci=[1 0]; % Condiciones iniciales


[t,y]=ode23('ode_canal',[0 15],ci); % Integrador

c1=y(:,1); % Cerrado 1
c2=y(:,2); % Cerrado 2
o=1-c1-c2; % Abierto


% plot (c1,c2,'g') % Grafico del diagrama de fases
% xlabel('c_1')
% ylabel('c_2')


figure
plot(t,c1,t,c2,t,o) % Grafico de la evolucion temporal
hold on
xlabel('t')


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

l=150;   %Longitud del RW
n=100;   %cantidad de canales

for i=1:n
    w(i,:)=rw_canal(1,.3,1.2,.4,.1,l); % pone en la fila i de w el RW del i-esimo canal
                                       % con parametros k1, k1p, k2, k2p,
                                       % delta t y longitud l
                                       % respectivamente
end


% El siguiente loop cuenta para cada columna z (representa el tiempo)
% cuantos canales estan en cada estado

for i=1:l
    c1(i)=0;
    c2(i)=0;
    o(i)=0;
    for z=1:n
        if w(z,i)==0
            c1(i)=c1(i)+1;
        end
        if w(z,i)==1
            c2(i)=c2(i)+1;
        end
        if w(z,i)==2
            o(i)=o(i)+1;
        end

    end
end

x=0.1:.1:15;  % escalea el tiempo con delta t en este caso .1


plot(x,c1/n,x,c2/n,x,o/n) % plotea la probabilidad de cada estado 
                          % (como la cantidad de canales en cada estado
                          % sobre el total de canales)

hold off