membraneODE_ca.m

function dy = membraneODE_ca(t, y) %%%HAS EVERY CHANNEL EXCEPT FOR THE CALCIUM-DEP POTASSIUM
[g, e, elecparam, calciump] = loadglobal; %load params. g(i) is conductance, e(j) is reversal potential for i,jth ion
dy = zeros(12,1);
if t>=elecparam(5) && t<elecparam(6)
dy(1) = (elecparam(1) ... %external current
        - g(1)*(y(2)^3).* y(3)*(y(1)-e(1)) ... %SODIUM
        - g(2)*(y(4)^4).*(y(1)-e(2)) ... %FAST POTASSIUM
        - g(3)*(y(5)^4).*(y(1)-e(2))...%SLOW POTASSIUM
        - g(4)*y(6)*y(7).*(y(1)-e(2)) ...%A-TYPE POTASSIUM
        - g(5)*(y(8)^2)*y(9).*(y(1)-100) ...%T-TYPE CALCIUM
        - g(6)*(y(10)^2)*y(11).*(y(1)-100) ...%N-TYPE CALCIUM 
        - g(7)*(y(12)^2).*(y(1)-100))/elecparam(2); %L-TYPE CALCIUM
else
    dy(1) = (0 ... %external current
        - g(1)*(y(2)^3).* y(3)*(y(1)-e(1)) ... %SODIUM
        - g(2)*(y(4)^4).*(y(1)-e(2)) ... %FAST POTASSIUM
        - g(3)*(y(5)^4).*(y(1)-e(2))...%SLOW POTASSIUM
        - g(4)*y(6)*y(7).*(y(1)-e(2)) ...%A-TYPE POTASSIUM
        - g(5)*(y(8)^2)*y(9).*(y(1)-100) ...%T-TYPE CALCIUM
        - g(6)*(y(10)^2)*y(11).*(y(1)-100) ...%N-TYPE CALCIUM 
        - g(7)*(y(12)^2).*(y(1)-100))/elecparam(2); %L-TYPE CALCIUM
end %end if
        

dy(2) = a_m(y(1)).*(1-y(2)) - b_m(y(1)).* y(2); %dm/dt SODIUM
dy(3) = a_h(y(1)).*(1-y(3)) - b_h(y(1)).* y(3); %dh/dt SODIUM
dy(4) = a_nf(y(1)).*(1-y(4)) - b_nf(y(1)).* y(4); %dnf/dt FAST POTASSIUM
dy(5) = a_ns(y(1)).*(1-y(5)) - b_ns(y(1)).* y(5); %dns/dt SLOW POTASSIUM 
dy(6) = ((1/(1+a_k(y(1))))-y(6)) / ((b_k(y(1)/(0.015*(1+a_k(y(1))))))); %dk/dt A-TYPE POTASSIUM
dy(7) = ((1/(1+a_l(y(1))))-y(7)) / ((b_l(y(1))/(0.015*(1+a_l(y(1)))))); %dl/dt A-TYPE POTASSIUM
dy(8) = a_a(y(1)).*(1-y(8)) - b_a(y(1)).* y(8); %da/dt T-TYPE CALCIUM 
dy(9) = a_b(y(1)).*(1-y(9)) - b_b(y(1)).* y(9); %db/dt T-TYPE CALCIUM [WORKS--adding this increases current tolerance]
dy(10) = a_c(y(1)).*(1-y(10)) - b_c(y(1)).* y(10); %dc/dt N-TYPE CALCIUM 
dy(11) = a_d(y(1)).*(1-y(11)) - b_d(y(1)).* y(11); %dd/dt N-TYPE CALCIUM [WORKS--adding this introduces fast bursting]
dy(12) = a_e(y(1)).*(1-y(12)) - b_e(y(1)).* y(12);%de/dt L-TYPE CALCIUM [WORKS--adding this doesn't change much, more stability at higher currents?]

end % end function