Changeset 27233 in project


Ignore:
Timestamp:
08/09/12 09:20:06 (9 years ago)
Author:
Ivan Raikov
Message:

nemo: modifications to some parts of the CGC example

Location:
release/4/nemo/trunk/examples/CGC/parts
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • release/4/nemo/trunk/examples/CGC/parts/CGC_CaHVA.m

    r27218 r27233  
    1 global CaHVA_h CaHVA_m comp307_vcbdur comp307_vchdur comp307_vcsteps comp307_vcinc comp307_vcbase comp307_vchold comp33_e comp33_gbar ica eca v comp307_vcinc comp33_gbar comp33_Q10 comp307_vchold fix_celsius comp33_Abeta_s comp307_vcbdur comp33_Kalpha_u comp307_vcsteps comp33_Kalpha_s comp33_Kbeta_s comp33_Kbeta_u comp33_Abeta_u comp307_vcbase comp307_vchdur comp33_V0alpha_u comp33_V0alpha_s comp33_e comp33_V0beta_s comp33_V0beta_u comp33_Aalpha_s comp33_Aalpha_u i_CaHVA
     1global CaHVA_h CaHVA_m comp307_vcbdur comp307_vchdur comp307_vcsteps comp307_vcinc comp307_vcbase comp307_vchold comp33_e comp33_gbar ica eca v comp33_V0alpha_s comp307_vcinc comp307_vchdur fix_celsius comp33_gbar comp33_Kalpha_s comp33_Kalpha_u comp307_vcbase comp33_V0beta_s comp307_vcbdur comp307_vchold comp33_Q10 comp33_e comp33_V0alpha_u comp33_Kbeta_u comp33_Aalpha_s comp33_Kbeta_s comp33_Aalpha_u comp33_V0beta_u comp307_vcsteps comp33_Abeta_s comp33_Abeta_u i_CaHVA
    22
    33
    44function dy = CGC_CaHVA (y, t)
    5   global CaHVA_h CaHVA_m comp307_vcbdur comp307_vchdur comp307_vcsteps comp307_vcinc comp307_vcbase comp307_vchold comp33_e comp33_gbar ica eca v comp307_vcinc comp33_gbar comp33_Q10 comp307_vchold fix_celsius comp33_Abeta_s comp307_vcbdur comp33_Kalpha_u comp307_vcsteps comp33_Kalpha_s comp33_Kbeta_s comp33_Kbeta_u comp33_Abeta_u comp307_vcbase comp307_vchdur comp33_V0alpha_u comp33_V0alpha_s comp33_e comp33_V0beta_s comp33_V0beta_u comp33_Aalpha_s comp33_Aalpha_u i_CaHVA
     5  global CaHVA_h CaHVA_m comp307_vcbdur comp307_vchdur comp307_vcsteps comp307_vcinc comp307_vcbase comp307_vchold comp33_e comp33_gbar ica eca v comp33_V0alpha_s comp307_vcinc comp307_vchdur fix_celsius comp33_gbar comp33_Kalpha_s comp33_Kalpha_u comp307_vcbase comp33_V0beta_s comp307_vcbdur comp307_vchold comp33_Q10 comp33_e comp33_V0alpha_u comp33_Kbeta_u comp33_Aalpha_s comp33_Kbeta_s comp33_Aalpha_u comp33_V0beta_u comp307_vcsteps comp33_Abeta_s comp33_Abeta_u i_CaHVA
    66  v  =  y(1);
     7  % rate eq CaHVA_hO
     8  CaHVA_hO  =  y(2);
    79  % rate eq CaHVA_mO
    8   CaHVA_mO  =  y(2);
    9   % rate eq CaHVA_hO
    10   CaHVA_hO  =  y(3);
     10  CaHVA_mO  =  y(3);
     11  % equation for CaHVA_m
     12  CaHVA_m  =  CaHVA_mO;
    1113  % equation for CaHVA_h
    1214  CaHVA_h  =  CaHVA_hO;
    13   % equation for CaHVA_m
    14   CaHVA_m  =  CaHVA_mO;
    1515  dy  =  zeros(length(y), 1);
     16  % state CaHVA_hO
     17  v452  =  CaHVA_hO;;
     18dy(2)  =  -(CaHVA_hO * comp33_beta_u(v)) + (1 - v452) * (comp33_alpha_u(v));
    1619  % state CaHVA_mO
    17   v452  =  CaHVA_mO;;
    18 dy(2)  =  -(CaHVA_mO * comp33_beta_s(v)) + (1 - v452) * (comp33_alpha_s(v));
    19   % state CaHVA_hO
    20   v455  =  CaHVA_hO;;
    21 dy(3)  =  -(CaHVA_hO * comp33_beta_u(v)) + (1 - v455) * (comp33_alpha_u(v));
     20  v455  =  CaHVA_mO;;
     21dy(3)  =  -(CaHVA_mO * comp33_beta_s(v)) + (1 - v455) * (comp33_alpha_s(v));
    2222  v449  =  CaHVA_m;;
    2323i_CaHVA  =  (comp33_gbar * v449 * v449 * CaHVA_h) * (v - comp33_e);
     
    3030
    3131function y0 = CGC_CaHVA_init (Vinit)
    32   global CaHVA_h CaHVA_m comp307_vcbdur comp307_vchdur comp307_vcsteps comp307_vcinc comp307_vcbase comp307_vchold comp33_e comp33_gbar ica eca v comp307_vcinc comp33_gbar comp33_Q10 comp307_vchold fix_celsius comp33_Abeta_s comp307_vcbdur comp33_Kalpha_u comp307_vcsteps comp33_Kalpha_s comp33_Kbeta_s comp33_Kbeta_u comp33_Abeta_u comp307_vcbase comp307_vchdur comp33_V0alpha_u comp33_V0alpha_s comp33_e comp33_V0beta_s comp33_V0beta_u comp33_Aalpha_s comp33_Aalpha_u i_CaHVA
     32  global CaHVA_h CaHVA_m comp307_vcbdur comp307_vchdur comp307_vcsteps comp307_vcinc comp307_vcbase comp307_vchold comp33_e comp33_gbar ica eca v comp33_V0alpha_s comp307_vcinc comp307_vchdur fix_celsius comp33_gbar comp33_Kalpha_s comp33_Kalpha_u comp307_vcbase comp33_V0beta_s comp307_vcbdur comp307_vchold comp33_Q10 comp33_e comp33_V0alpha_u comp33_Kbeta_u comp33_Aalpha_s comp33_Kbeta_s comp33_Aalpha_u comp33_V0beta_u comp307_vcsteps comp33_Abeta_s comp33_Abeta_u i_CaHVA
    3333  y0  =  zeros(3, 1);
    3434  v  =  Vinit;
    3535  y0(1)  =  v;
     36  comp33_V0alpha_s  =  -29.06;
    3637  comp307_vcinc  =  10.0;
     38  comp307_vchdur  =  30.0;
     39  fix_celsius  =  30.0;
    3740  comp33_gbar  =  0.00046;
     41  comp33_Kalpha_s  =  15.87301587302;
     42  comp33_Kalpha_u  =  -18.183;
     43  comp307_vcbase  =  -69.0;
     44  comp33_V0beta_s  =  -18.66;
     45  comp307_vcbdur  =  100.0;
     46  comp307_vchold  =  -71.0;
    3847  comp33_Q10  =  3.0;
    39   comp307_vchold  =  -71.0;
    40   fix_celsius  =  30.0;
     48  comp33_e  =  129.33;
     49  comp33_V0alpha_u  =  -48.0;
     50  comp33_Kbeta_u  =  83.33;
     51  comp33_Aalpha_s  =  0.04944;
     52  comp33_Kbeta_s  =  -25.641;
     53  comp33_Aalpha_u  =  0.0013;
     54  comp33_V0beta_u  =  -48.0;
     55  comp307_vcsteps  =  8.0;
    4156  comp33_Abeta_s  =  0.08298;
    42   comp307_vcbdur  =  100.0;
    43   comp33_Kalpha_u  =  -18.183;
    44   comp307_vcsteps  =  8.0;
    45   comp33_Kalpha_s  =  15.87301587302;
    46   comp33_Kbeta_s  =  -25.641;
    47   comp33_Kbeta_u  =  83.33;
    4857  comp33_Abeta_u  =  0.0013;
    49   comp307_vcbase  =  -69.0;
    50   comp307_vchdur  =  30.0;
    51   comp33_V0alpha_u  =  -48.0;
    52   comp33_V0alpha_s  =  -29.06;
    53   comp33_e  =  129.33;
    54   comp33_V0beta_s  =  -18.66;
    55   comp33_V0beta_u  =  -48.0;
    56   comp33_Aalpha_s  =  0.04944;
    57   comp33_Aalpha_u  =  0.0013;
     58  CaHVA_h  =  (comp33_alpha_u(v)) / (comp33_alpha_u(v) + comp33_beta_u(v));
     59  CaHVA_hO  =  CaHVA_h;
    5860  CaHVA_m  =  (comp33_alpha_s(v)) / (comp33_alpha_s(v) + comp33_beta_s(v));
    5961  CaHVA_mO  =  CaHVA_m;
    60   CaHVA_h  =  (comp33_alpha_u(v)) / (comp33_alpha_u(v) + comp33_beta_u(v));
    61   CaHVA_hO  =  CaHVA_h;
    62   y0(2)  =  CaHVA_mO;
    63   y0(3)  =  CaHVA_hO;
     62  y0(2)  =  CaHVA_hO;
     63  y0(3)  =  CaHVA_mO;
    6464  v449  =  CaHVA_m;;
    6565i_CaHVA  =  (comp33_gbar * v449 * v449 * CaHVA_h) * (v - comp33_e);
     
    7272
    7373
    74 function retval456 = comp33_beta_u (v)
    75   global CaHVA_h CaHVA_m comp307_vcbdur comp307_vchdur comp307_vcsteps comp307_vcinc comp307_vcbase comp307_vchold comp33_e comp33_gbar ica eca comp307_vcinc comp33_gbar comp33_Q10 comp307_vchold fix_celsius comp33_Abeta_s comp307_vcbdur comp33_Kalpha_u comp307_vcsteps comp33_Kalpha_s comp33_Kbeta_s comp33_Kbeta_u comp33_Abeta_u comp307_vcbase comp307_vchdur comp33_V0alpha_u comp33_V0alpha_s comp33_e comp33_V0beta_s comp33_V0beta_u comp33_Aalpha_s comp33_Aalpha_u i_CaHVA
    76   retval456  =  comp33_Q10 * comp33_Abeta_u * exp((v + -(comp33_V0beta_u)) / comp33_Kbeta_u);
     74function retval456 = comp33_alpha_u (v)
     75  global CaHVA_h CaHVA_m comp307_vcbdur comp307_vchdur comp307_vcsteps comp307_vcinc comp307_vcbase comp307_vchold comp33_e comp33_gbar ica eca comp33_V0alpha_s comp307_vcinc comp307_vchdur fix_celsius comp33_gbar comp33_Kalpha_s comp33_Kalpha_u comp307_vcbase comp33_V0beta_s comp307_vcbdur comp307_vchold comp33_Q10 comp33_e comp33_V0alpha_u comp33_Kbeta_u comp33_Aalpha_s comp33_Kbeta_s comp33_Aalpha_u comp33_V0beta_u comp307_vcsteps comp33_Abeta_s comp33_Abeta_u i_CaHVA
     76  retval456  =  comp33_Q10 * comp33_Aalpha_u * exp((v + -(comp33_V0alpha_u)) / comp33_Kalpha_u);
    7777end
    7878
     
    8080
    8181
    82 function retval457 = comp33_beta_s (v)
    83   global CaHVA_h CaHVA_m comp307_vcbdur comp307_vchdur comp307_vcsteps comp307_vcinc comp307_vcbase comp307_vchold comp33_e comp33_gbar ica eca comp307_vcinc comp33_gbar comp33_Q10 comp307_vchold fix_celsius comp33_Abeta_s comp307_vcbdur comp33_Kalpha_u comp307_vcsteps comp33_Kalpha_s comp33_Kbeta_s comp33_Kbeta_u comp33_Abeta_u comp307_vcbase comp307_vchdur comp33_V0alpha_u comp33_V0alpha_s comp33_e comp33_V0beta_s comp33_V0beta_u comp33_Aalpha_s comp33_Aalpha_u i_CaHVA
    84   retval457  =  comp33_Q10 * comp33_Abeta_s * exp((v + -(comp33_V0beta_s)) / comp33_Kbeta_s);
     82function retval457 = comp33_alpha_s (v)
     83  global CaHVA_h CaHVA_m comp307_vcbdur comp307_vchdur comp307_vcsteps comp307_vcinc comp307_vcbase comp307_vchold comp33_e comp33_gbar ica eca comp33_V0alpha_s comp307_vcinc comp307_vchdur fix_celsius comp33_gbar comp33_Kalpha_s comp33_Kalpha_u comp307_vcbase comp33_V0beta_s comp307_vcbdur comp307_vchold comp33_Q10 comp33_e comp33_V0alpha_u comp33_Kbeta_u comp33_Aalpha_s comp33_Kbeta_s comp33_Aalpha_u comp33_V0beta_u comp307_vcsteps comp33_Abeta_s comp33_Abeta_u i_CaHVA
     84  retval457  =  comp33_Q10 * comp33_Aalpha_s * exp((v + -(comp33_V0alpha_s)) / comp33_Kalpha_s);
    8585end
    8686
     
    8888
    8989
    90 function retval458 = comp33_alpha_s (v)
    91   global CaHVA_h CaHVA_m comp307_vcbdur comp307_vchdur comp307_vcsteps comp307_vcinc comp307_vcbase comp307_vchold comp33_e comp33_gbar ica eca comp307_vcinc comp33_gbar comp33_Q10 comp307_vchold fix_celsius comp33_Abeta_s comp307_vcbdur comp33_Kalpha_u comp307_vcsteps comp33_Kalpha_s comp33_Kbeta_s comp33_Kbeta_u comp33_Abeta_u comp307_vcbase comp307_vchdur comp33_V0alpha_u comp33_V0alpha_s comp33_e comp33_V0beta_s comp33_V0beta_u comp33_Aalpha_s comp33_Aalpha_u i_CaHVA
    92   retval458  =  comp33_Q10 * comp33_Aalpha_s * exp((v + -(comp33_V0alpha_s)) / comp33_Kalpha_s);
     90function retval458 = comp33_beta_u (v)
     91  global CaHVA_h CaHVA_m comp307_vcbdur comp307_vchdur comp307_vcsteps comp307_vcinc comp307_vcbase comp307_vchold comp33_e comp33_gbar ica eca comp33_V0alpha_s comp307_vcinc comp307_vchdur fix_celsius comp33_gbar comp33_Kalpha_s comp33_Kalpha_u comp307_vcbase comp33_V0beta_s comp307_vcbdur comp307_vchold comp33_Q10 comp33_e comp33_V0alpha_u comp33_Kbeta_u comp33_Aalpha_s comp33_Kbeta_s comp33_Aalpha_u comp33_V0beta_u comp307_vcsteps comp33_Abeta_s comp33_Abeta_u i_CaHVA
     92  retval458  =  comp33_Q10 * comp33_Abeta_u * exp((v + -(comp33_V0beta_u)) / comp33_Kbeta_u);
    9393end
    9494
     
    9696
    9797
    98 function retval459 = comp33_alpha_u (v)
    99   global CaHVA_h CaHVA_m comp307_vcbdur comp307_vchdur comp307_vcsteps comp307_vcinc comp307_vcbase comp307_vchold comp33_e comp33_gbar ica eca comp307_vcinc comp33_gbar comp33_Q10 comp307_vchold fix_celsius comp33_Abeta_s comp307_vcbdur comp33_Kalpha_u comp307_vcsteps comp33_Kalpha_s comp33_Kbeta_s comp33_Kbeta_u comp33_Abeta_u comp307_vcbase comp307_vchdur comp33_V0alpha_u comp33_V0alpha_s comp33_e comp33_V0beta_s comp33_V0beta_u comp33_Aalpha_s comp33_Aalpha_u i_CaHVA
    100   retval459  =  comp33_Q10 * comp33_Aalpha_u * exp((v + -(comp33_V0alpha_u)) / comp33_Kalpha_u);
     98function retval459 = comp33_beta_s (v)
     99  global CaHVA_h CaHVA_m comp307_vcbdur comp307_vchdur comp307_vcsteps comp307_vcinc comp307_vcbase comp307_vchold comp33_e comp33_gbar ica eca comp33_V0alpha_s comp307_vcinc comp307_vchdur fix_celsius comp33_gbar comp33_Kalpha_s comp33_Kalpha_u comp307_vcbase comp33_V0beta_s comp307_vcbdur comp307_vchold comp33_Q10 comp33_e comp33_V0alpha_u comp33_Kbeta_u comp33_Aalpha_s comp33_Kbeta_s comp33_Aalpha_u comp33_V0beta_u comp307_vcsteps comp33_Abeta_s comp33_Abeta_u i_CaHVA
     100  retval459  =  comp33_Q10 * comp33_Abeta_s * exp((v + -(comp33_V0beta_s)) / comp33_Kbeta_s);
    101101end
    102102
  • release/4/nemo/trunk/examples/CGC/parts/CGC_KCa.m

    r27218 r27233  
    1 global KCa_m comp47_e comp47_gbar ik ek cai v ica comp47_gbar comp47_Kalpha_c comp47_Kbeta_c Vrest comp47_Q10 fix_celsius comp47_Balpha_c comp47_e comp47_Abeta_c comp47_Aalpha_c comp47_Bbeta_c KCa_m_inf KCa_m_tau i_KCa
     1global KCa_m comp209_vcbdur comp209_vchdur comp209_vcsteps comp209_vcinc comp209_vcbase comp209_vchold comp47_e comp47_gbar ik ek ica v cai comp47_Balpha_c fix_celsius Vrest comp209_vcbdur comp47_Kalpha_c comp47_Abeta_c comp47_Bbeta_c comp47_Kbeta_c comp209_vcinc comp47_Q10 comp47_Aalpha_c comp209_vchdur comp209_vcbase comp209_vcsteps comp47_gbar comp47_e comp209_vchold KCa_m_inf KCa_m_tau i_KCa
    22
    33
    44function dy = CGC_KCa (y, t)
    5   global KCa_m comp47_e comp47_gbar ik ek cai v ica comp47_gbar comp47_Kalpha_c comp47_Kbeta_c Vrest comp47_Q10 fix_celsius comp47_Balpha_c comp47_e comp47_Abeta_c comp47_Aalpha_c comp47_Bbeta_c KCa_m_inf KCa_m_tau i_KCa
     5  global KCa_m comp209_vcbdur comp209_vchdur comp209_vcsteps comp209_vcinc comp209_vcbase comp209_vchold comp47_e comp47_gbar ik ek ica v cai comp47_Balpha_c fix_celsius Vrest comp209_vcbdur comp47_Kalpha_c comp47_Abeta_c comp47_Bbeta_c comp47_Kbeta_c comp209_vcinc comp47_Q10 comp47_Aalpha_c comp209_vchdur comp209_vcbase comp209_vcsteps comp47_gbar comp47_e comp209_vchold KCa_m_inf KCa_m_tau i_KCa
    66  v  =  y(1);
    77  % rate eq KCa_m
     
    2323
    2424function y0 = CGC_KCa_init (Vinit)
    25   global KCa_m comp47_e comp47_gbar ik ek cai v ica comp47_gbar comp47_Kalpha_c comp47_Kbeta_c Vrest comp47_Q10 fix_celsius comp47_Balpha_c comp47_e comp47_Abeta_c comp47_Aalpha_c comp47_Bbeta_c KCa_m_inf KCa_m_tau i_KCa
     25  global KCa_m comp209_vcbdur comp209_vchdur comp209_vcsteps comp209_vcinc comp209_vcbase comp209_vchold comp47_e comp47_gbar ik ek ica v cai comp47_Balpha_c fix_celsius Vrest comp209_vcbdur comp47_Kalpha_c comp47_Abeta_c comp47_Bbeta_c comp47_Kbeta_c comp209_vcinc comp47_Q10 comp47_Aalpha_c comp209_vchdur comp209_vcbase comp209_vcsteps comp47_gbar comp47_e comp209_vchold KCa_m_inf KCa_m_tau i_KCa
    2626  y0  =  zeros(2, 1);
    2727  v  =  Vinit;
    2828  y0(1)  =  v;
     29  comp47_Balpha_c  =  0.0015;
     30  fix_celsius  =  30.0;
     31  Vrest  =  -68.0;
     32  comp209_vcbdur  =  100.0;
     33  comp47_Kalpha_c  =  -11.765;
     34  comp47_Abeta_c  =  1.5;
     35  comp47_Bbeta_c  =  0.00015;
     36  comp47_Kbeta_c  =  -11.765;
     37  comp209_vcinc  =  10.0;
     38  comp47_Q10  =  1.0;
     39  comp47_Aalpha_c  =  2.5;
     40  comp209_vchdur  =  30.0;
     41  comp209_vcbase  =  -69.0;
     42  comp209_vcsteps  =  8.0;
    2943  comp47_gbar  =  0.004;
    30   comp47_Kalpha_c  =  -11.765;
    31   comp47_Kbeta_c  =  -11.765;
    32   Vrest  =  -68.0;
    33   comp47_Q10  =  1.0;
    34   fix_celsius  =  30.0;
    35   comp47_Balpha_c  =  0.0015;
    3644  comp47_e  =  -84.69;
    37   comp47_Abeta_c  =  1.5;
    38   comp47_Aalpha_c  =  2.5;
    39   comp47_Bbeta_c  =  0.00015;
     45  comp209_vchold  =  -71.0;
    4046  KCa_m_inf  =  comp47_alpha_c(v, cai);
    4147  KCa_m_tau  =  comp47_beta_c(v, cai);
     
    5157
    5258
    53 function retval251 = comp47_alpha_c (v, cai)
    54   global KCa_m comp47_e comp47_gbar ik ek ica comp47_gbar comp47_Kalpha_c comp47_Kbeta_c Vrest comp47_Q10 fix_celsius comp47_Balpha_c comp47_e comp47_Abeta_c comp47_Aalpha_c comp47_Bbeta_c KCa_m_inf KCa_m_tau i_KCa
    55   retval251  =  (comp47_Q10 * comp47_Aalpha_c) / (1.0 + (comp47_Balpha_c * exp(v / comp47_Kalpha_c)) / cai);
     59function retval336 = comp47_alpha_c (v, cai)
     60  global KCa_m comp209_vcbdur comp209_vchdur comp209_vcsteps comp209_vcinc comp209_vcbase comp209_vchold comp47_e comp47_gbar ik ek ica comp47_Balpha_c fix_celsius Vrest comp209_vcbdur comp47_Kalpha_c comp47_Abeta_c comp47_Bbeta_c comp47_Kbeta_c comp209_vcinc comp47_Q10 comp47_Aalpha_c comp209_vchdur comp209_vcbase comp209_vcsteps comp47_gbar comp47_e comp209_vchold KCa_m_inf KCa_m_tau i_KCa
     61  retval336  =  (comp47_Q10 * comp47_Aalpha_c) / (1.0 + (comp47_Balpha_c * exp(v / comp47_Kalpha_c)) / cai);
    5662end
    5763
     
    5965
    6066
    61 function retval252 = comp47_beta_c (v, cai)
    62   global KCa_m comp47_e comp47_gbar ik ek ica comp47_gbar comp47_Kalpha_c comp47_Kbeta_c Vrest comp47_Q10 fix_celsius comp47_Balpha_c comp47_e comp47_Abeta_c comp47_Aalpha_c comp47_Bbeta_c KCa_m_inf KCa_m_tau i_KCa
    63   retval252  =  (comp47_Q10 * comp47_Abeta_c) / (1.0 + cai / (comp47_Bbeta_c * exp(v / comp47_Kbeta_c)));
     67function retval337 = comp47_beta_c (v, cai)
     68  global KCa_m comp209_vcbdur comp209_vchdur comp209_vcsteps comp209_vcinc comp209_vcbase comp209_vchold comp47_e comp47_gbar ik ek ica comp47_Balpha_c fix_celsius Vrest comp209_vcbdur comp47_Kalpha_c comp47_Abeta_c comp47_Bbeta_c comp47_Kbeta_c comp209_vcinc comp47_Q10 comp47_Aalpha_c comp209_vchdur comp209_vcbase comp209_vcsteps comp47_gbar comp47_e comp209_vchold KCa_m_inf KCa_m_tau i_KCa
     69  retval337  =  (comp47_Q10 * comp47_Abeta_c) / (1.0 + cai / (comp47_Bbeta_c * exp(v / comp47_Kbeta_c)));
    6470end
    6571
  • release/4/nemo/trunk/examples/CGC/parts/CGC_KV.m

    r27218 r27233  
    1 global KV_m comp223_vcbdur comp223_vchdur comp223_vcsteps comp223_vcinc comp223_vcbase comp223_vchold comp63_e comp63_gbar ik ek ica cai v comp223_vcbdur comp63_Abeta_n comp223_vcinc comp223_vchold Vrest fix_celsius comp63_Q10 comp63_Kbeta_n comp223_vchdur comp63_gbar comp63_Kalpha_n comp63_V0beta_n comp223_vcbase comp63_V0alpha_n comp63_Aalpha_n comp223_vcsteps comp63_e i_KV
     1global 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
    22
    33
    44function dy = CGC_KV (y, t)
    5   global KV_m comp223_vcbdur comp223_vchdur comp223_vcsteps comp223_vcinc comp223_vcbase comp223_vchold comp63_e comp63_gbar ik ek ica cai v comp223_vcbdur comp63_Abeta_n comp223_vcinc comp223_vchold Vrest fix_celsius comp63_Q10 comp63_Kbeta_n comp223_vchdur comp63_gbar comp63_Kalpha_n comp63_V0beta_n comp223_vcbase comp63_V0alpha_n comp63_Aalpha_n comp223_vcsteps comp63_e i_KV
     5  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
    66  v  =  y(1);
    77  % rate eq KV_mO
     
    2323
    2424function y0 = CGC_KV_init (Vinit)
    25   global KV_m comp223_vcbdur comp223_vchdur comp223_vcsteps comp223_vcinc comp223_vcbase comp223_vchold comp63_e comp63_gbar ik ek ica cai v comp223_vcbdur comp63_Abeta_n comp223_vcinc comp223_vchold Vrest fix_celsius comp63_Q10 comp63_Kbeta_n comp223_vchdur comp63_gbar comp63_Kalpha_n comp63_V0beta_n comp223_vcbase comp63_V0alpha_n comp63_Aalpha_n comp223_vcsteps comp63_e i_KV
     25  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
    2626  y0  =  zeros(2, 1);
    2727  v  =  Vinit;
    2828  y0(1)  =  v;
     29  comp63_V0alpha_n  =  -25.0;
     30  comp63_e  =  -84.69;
     31  fix_celsius  =  30.0;
     32  comp63_Kalpha_n  =  -10.0;
    2933  comp223_vcbdur  =  100.0;
     34  Vrest  =  -68.0;
     35  comp63_Kbeta_n  =  -80.0;
     36  comp223_vchold  =  -71.0;
    3037  comp63_Abeta_n  =  0.125;
     38  comp63_gbar  =  0.003;
     39  comp63_Aalpha_n  =  -0.01;
     40  comp63_Q10  =  13.5137964673603;
     41  comp223_vchdur  =  30.0;
    3142  comp223_vcinc  =  10.0;
    32   comp223_vchold  =  -71.0;
    33   Vrest  =  -68.0;
    34   fix_celsius  =  30.0;
    35   comp63_Q10  =  13.5137964673603;
    36   comp63_Kbeta_n  =  -80.0;
    37   comp223_vchdur  =  30.0;
    38   comp63_gbar  =  0.003;
    39   comp63_Kalpha_n  =  -10.0;
     43  comp223_vcsteps  =  8.0;
    4044  comp63_V0beta_n  =  -35.0;
    4145  comp223_vcbase  =  -69.0;
    42   comp63_V0alpha_n  =  -25.0;
    43   comp63_Aalpha_n  =  -0.01;
    44   comp223_vcsteps  =  8.0;
    45   comp63_e  =  -84.69;
    4646  KV_m  =  (comp63_alpha_n(v)) / (comp63_alpha_n(v) + comp63_beta_n(v));
    4747  KV_mO  =  KV_m;
     
    5757
    5858
    59 function retval355 = comp63_beta_n (v)
    60      global KV_m comp223_vcbdur comp223_vchdur comp223_vcsteps comp223_vcinc comp223_vcbase comp223_vchold comp63_e comp63_gbar ik ek ica cai v comp223_vcbdur comp63_Abeta_n comp223_vcinc comp223_vchold Vrest fix_celsius comp63_Q10 comp63_Kbeta_n comp223_vchdur comp63_gbar comp63_Kalpha_n comp63_V0beta_n comp223_vcbase comp63_V0alpha_n comp63_Aalpha_n comp223_vcsteps comp63_e i_KV
    61 
    62   retval355  =  comp63_Q10 * comp63_Abeta_n * exp((v + -(comp63_V0beta_n)) / comp63_Kbeta_n);
     59function retval355 = comp63_alpha_n (v)
     60  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
     61  retval355  =  comp63_Q10 * comp63_Aalpha_n * linoid(v + -(comp63_V0alpha_n), comp63_Kalpha_n);
    6362end
    6463
     
    6665
    6766
    68 function retval356 = comp63_alpha_n (v)
    69      global KV_m comp223_vcbdur comp223_vchdur comp223_vcsteps comp223_vcinc comp223_vcbase comp223_vchold comp63_e comp63_gbar ik ek ica cai v comp223_vcbdur comp63_Abeta_n comp223_vcinc comp223_vchold Vrest fix_celsius comp63_Q10 comp63_Kbeta_n comp223_vchdur comp63_gbar comp63_Kalpha_n comp63_V0beta_n comp223_vcbase comp63_V0alpha_n comp63_Aalpha_n comp223_vcsteps comp63_e i_KV
    70 
    71   retval356  =  comp63_Q10 * comp63_Aalpha_n * linoid(v + -(comp63_V0alpha_n), comp63_Kalpha_n);
     67function retval356 = comp63_beta_n (v)
     68  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
     69  retval356  =  comp63_Q10 * comp63_Abeta_n * exp((v + -(comp63_V0beta_n)) / comp63_Kbeta_n);
    7270end
    7371
     
    7674
    7775function retval357 = linoid (x, y)
     76  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
    7877  if (abs(x / y) < 1e-06)
    7978      v359  =  y * (1.0 + -(x / y / 2.0));
  • release/4/nemo/trunk/examples/CGC/parts/CGC_KV_vclamp.m

    r27218 r27233  
    4848  lsode_options ("absolute tolerance" , abstol) ;
    4949  lsode_options ("relative tolerance" , reltol) ;
     50  lsode_options ("integration method" , "bdf") ;
    5051  lsode_options ("initial step size" , dt) ;
    5152  t0 = 0.0 ;
     
    7172  ys = vertcat (ys , [ thold1' y ]) ;
    7273endfunction
    73 
    7474function ilog = CGC_KV_vclamp_fn (N , dt , vchold , vchdur , vcbase , vcdur , vcinc , vcsteps)
    7575  vc = vcbase ;
     
    9191  endfor ;
    9292endfunction
    93 
    9493global reltol abstol
    95 
    9694reltol = 1e-03 ;
    9795abstol = 1e-03 ;
    9896dt = 0.001 ;
    99 
    10097## i_KV plot
    10198i_KV_index = 2 ;
    10299global comp223_vchold comp223_vchdur comp223_vcbase comp223_vcbdur comp223_vcinc comp223_vcsteps
    103100CGC_KV_ilog = CGC_KV_vclamp_fn (N , dt , comp223_vchold , comp223_vchdur , comp223_vcbase , comp223_vcbdur , comp223_vcinc , comp223_vcsteps) ;
    104 
    105101subplot (1 , 3 , 1) ;
    106102plot (CGC_KV_ilog{1}(:,1) , CGC_KV_ilog{1}(:,i_KV_index) , CGC_KV_ilog{2}(:,1) , CGC_KV_ilog{2}(:,i_KV_index) , CGC_KV_ilog{3}(:,1) , CGC_KV_ilog{3}(:,i_KV_index) , CGC_KV_ilog{4}(:,1) , CGC_KV_ilog{4}(:,i_KV_index) , CGC_KV_ilog{5}(:,1) , CGC_KV_ilog{5}(:,i_KV_index) , CGC_KV_ilog{6}(:,1) , CGC_KV_ilog{6}(:,i_KV_index) , CGC_KV_ilog{7}(:,1) , CGC_KV_ilog{7}(:,i_KV_index) , CGC_KV_ilog{8}(:,1) , CGC_KV_ilog{8}(:,i_KV_index)) ;
  • release/4/nemo/trunk/examples/CGC/parts/plot_vclamp.m

    r27218 r27233  
    5959title ("pNa current");
    6060
    61 print  ("NEURON_Vclamp.eps", "-depsc");
     61print  ("Octave_Vclamp.eps", "-depsc");
    6262
Note: See TracChangeset for help on using the changeset viewer.