global KV_m comp223_vcbdur comp223_vchdur comp223_vcsteps comp223_vcinc comp223_vcbase comp223_vchold comp63_e comp63_gbar ik ek ica v cai comp63_V0alpha_n comp63_e fix_celsius comp63_Kalpha_n comp223_vcbdur Vrest comp63_Kbeta_n comp223_vchold comp63_Abeta_n comp63_gbar comp63_Aalpha_n comp63_Q10 comp223_vchdur comp223_vcinc comp223_vcsteps comp63_V0beta_n comp223_vcbase i_KV


function dy = CGC_KV (y, t)
  global KV_m comp223_vcbdur comp223_vchdur comp223_vcsteps comp223_vcinc comp223_vcbase comp223_vchold comp63_e comp63_gbar ik ek ica v cai comp63_V0alpha_n comp63_e fix_celsius comp63_Kalpha_n comp223_vcbdur Vrest comp63_Kbeta_n comp223_vchold comp63_Abeta_n comp63_gbar comp63_Aalpha_n comp63_Q10 comp223_vchdur comp223_vcinc comp223_vcsteps comp63_V0beta_n comp223_vcbase i_KV
  v  =  y(1);
  % rate eq KV_mO
  KV_mO  =  y(2);
  % equation for KV_m
  KV_m  =  KV_mO;
  dy  =  zeros(length(y), 1);
  % state KV_mO
  v354  =  KV_mO;; 
dy(2)  =  -(KV_mO * comp63_beta_n(v)) + (1 - v354) * (comp63_alpha_n(v));
  v351  =  KV_m;; 
i_KV  =  (comp63_gbar * v351 * v351 * v351 * v351) * (v - comp63_e);
  ik  =  i_KV;
  dy(1)  =  0.0;


end


function y0 = CGC_KV_init (Vinit)
  global KV_m comp223_vcbdur comp223_vchdur comp223_vcsteps comp223_vcinc comp223_vcbase comp223_vchold comp63_e comp63_gbar ik ek ica v cai comp63_V0alpha_n comp63_e fix_celsius comp63_Kalpha_n comp223_vcbdur Vrest comp63_Kbeta_n comp223_vchold comp63_Abeta_n comp63_gbar comp63_Aalpha_n comp63_Q10 comp223_vchdur comp223_vcinc comp223_vcsteps comp63_V0beta_n comp223_vcbase i_KV
  y0  =  zeros(2, 1);
  v  =  Vinit;
  y0(1)  =  v;
  comp63_V0alpha_n  =  -25.0;
  comp63_e  =  -84.69;
  fix_celsius  =  30.0;
  comp63_Kalpha_n  =  -10.0;
  comp223_vcbdur  =  100.0;
  Vrest  =  -68.0;
  comp63_Kbeta_n  =  -80.0;
  comp223_vchold  =  -71.0;
  comp63_Abeta_n  =  0.125;
  comp63_gbar  =  0.003;
  comp63_Aalpha_n  =  -0.01;
  comp63_Q10  =  13.5137964673603;
  comp223_vchdur  =  30.0;
  comp223_vcinc  =  10.0;
  comp223_vcsteps  =  8.0;
  comp63_V0beta_n  =  -35.0;
  comp223_vcbase  =  -69.0;
  KV_m  =  (comp63_alpha_n(v)) / (comp63_alpha_n(v) + comp63_beta_n(v));
  KV_mO  =  KV_m;
  y0(2)  =  KV_mO;
  v351  =  KV_m;; 
i_KV  =  (comp63_gbar * v351 * v351 * v351 * v351) * (v - comp63_e);
  ik  =  i_KV;


end




function retval355 = comp63_alpha_n (v)
  global KV_m comp223_vcbdur comp223_vchdur comp223_vcsteps comp223_vcinc comp223_vcbase comp223_vchold comp63_e comp63_gbar ik ek ica cai comp63_V0alpha_n comp63_e fix_celsius comp63_Kalpha_n comp223_vcbdur Vrest comp63_Kbeta_n comp223_vchold comp63_Abeta_n comp63_gbar comp63_Aalpha_n comp63_Q10 comp223_vchdur comp223_vcinc comp223_vcsteps comp63_V0beta_n comp223_vcbase i_KV
  retval355  =  comp63_Q10 * comp63_Aalpha_n * linoid(v + -(comp63_V0alpha_n), comp63_Kalpha_n);
end




function retval356 = comp63_beta_n (v)
  global KV_m comp223_vcbdur comp223_vchdur comp223_vcsteps comp223_vcinc comp223_vcbase comp223_vchold comp63_e comp63_gbar ik ek ica cai comp63_V0alpha_n comp63_e fix_celsius comp63_Kalpha_n comp223_vcbdur Vrest comp63_Kbeta_n comp223_vchold comp63_Abeta_n comp63_gbar comp63_Aalpha_n comp63_Q10 comp223_vchdur comp223_vcinc comp223_vcsteps comp63_V0beta_n comp223_vcbase i_KV
  retval356  =  comp63_Q10 * comp63_Abeta_n * exp((v + -(comp63_V0beta_n)) / comp63_Kbeta_n);
end




function retval357 = linoid (x, y)
  global KV_m comp223_vcbdur comp223_vchdur comp223_vcsteps comp223_vcinc comp223_vcbase comp223_vchold comp63_e comp63_gbar ik ek ica v cai comp63_V0alpha_n comp63_e fix_celsius comp63_Kalpha_n comp223_vcbdur Vrest comp63_Kbeta_n comp223_vchold comp63_Abeta_n comp63_gbar comp63_Aalpha_n comp63_Q10 comp223_vchdur comp223_vcinc comp223_vcsteps comp63_V0beta_n comp223_vcbase i_KV
  if (abs(x / y) < 1e-06) 
      v359  =  y * (1.0 + -(x / y / 2.0)); 
      else 
        v359  =  x / (exp(x / y) + -1.0); 
      end; 
retval357  =  v359;
end