% Khaliq03 model driver for Matlab

global Ih_m K3_m K2_m K1_h K1_m CaP_m CaBK_h CaBK_m comp95_cac comp92_e_Leak comp92_gbar_Leak comp88_e_Ih comp88_gbar_Ih comp84_e comp84_gbar comp84_Na_z comp78_e_K3 comp78_gbar_K3 comp72_e_K2 comp72_gbar_K2 comp57_e_K1 comp57_gbar_K1 comp22_gbar_CaP comp0_e_CaBK comp0_gbar_CaBK comp0_CaBK_z ik ina i ek ena e ica cai cao v comp84_Na_gamma comp84_Na_epsilon comp72_e_K2 comp84_Na_Oon comp84_Na_zeta comp84_Na_x6 comp84_Na_x5 comp84_Na_x4 comp84_Na_x3 comp84_Na_x2 comp84_Na_x1 comp78_e_K3 comp84_Na_delta comp84_Na_alfac comp84_Na_Con comp84_gbar comp88_gbar_Ih comp57_e_K1 ca0 comp84_Na_alpha comp78_gbar_K3 comp72_gbar_K2 comp57_gbar_K1 comp0_gbar_CaBK comp95_F comp84_Na_Ooff comp22_cao comp0_e_CaBK comp0_CaBK_ztau comp95_ca_depth comp84_Na_btfac comp92_e_Leak comp84_Na_beta comp22_gmax_CaP comp84_Na_Coff comp84_e comp92_gbar_Leak comp95_ca_beta comp88_e_Ih comp95_cac comp84_b01 comp84_b02 comp84_b03 comp84_b04 comp72_K2_v comp84_bi1 comp84_bi2 comp84_bi3 comp84_bi4 comp84_bi5 comp84_b11 comp84_b12 comp84_b13 comp84_b14 comp84_b0O comp84_bin comp84_bip comp84_b1n comp22_CaP_tau comp88_Ih_tau comp84_f01 comp84_f02 comp84_f03 comp84_f04 comp84_fi1 comp84_fi2 comp84_fi3 comp84_fi4 comp84_fi5 comp84_f11 comp84_f12 comp84_f13 comp84_f14 comp0_CaBK_zinf comp84_f0O comp84_fin comp84_fip comp22_CaP_inf comp84_f1n comp88_Ih_inf comp78_K3_v comp57_K1_v comp0_CaBK_v comp22_gbar_CaP CaP_m_tau comp72_K2_minf comp57_K1_htau comp57_K1_mtau Ih_m_tau comp0_CaBK_htau comp0_CaBK_mtau CaP_m_inf comp57_K1_hinf comp78_K3_mtau comp57_K1_minf Ih_m_inf comp0_CaBK_hinf comp0_CaBK_minf comp78_K3_minf comp72_K2_mtau K3_m_inf CaBK_m_inf CaBK_h_tau K1_m_tau CaBK_h_inf K2_m_tau K1_m_inf K1_h_tau K3_m_tau K2_m_inf CaBK_m_tau K1_h_inf i_Leak i_Ih i_Narsg i_K3 i_K2 i_K1 i_CaP i_CaBK

y0 = Khaliq03_init(-65)

t0 = 0.0;
t1 = 2000.0;
dt = 1e-3;

reltol=1e-2;
abstol=1e-3;

P = odeset ('RelTol', reltol, 'AbsTol', abstol, 'InitialStep', dt, 'BDF', 'on');
[t,y] = ode15s (@Khaliq03_dy, [t0 t1], y0, P);