High frequency stimulation of the Subthalamic Nucleus (Rubin and Terman 2004)

 Download zip file 
Help downloading and running models
Accession:116867
" ... Using a computational model, this paper considers the hypothesis that DBS works by replacing pathologically rhythmic basal ganglia output with tonic, high frequency firing. In our simulations of parkinsonian conditions, rhythmic inhibition from GPi to the thalamus compromises the ability of thalamocortical relay (TC) cells to respond to depolarizing inputs, such as sensorimotor signals. High frequency stimulation of STN regularizes GPi firing, and this restores TC responsiveness, despite the increased frequency and amplitude of GPi inhibition to thalamus that result. We provide a mathematical phase plane analysis of the mechanisms that determine TC relay capabilities in normal, parkinsonian, and DBS states in a reduced model. This analysis highlights the differences in deinactivation of the low-threshold calcium T -current that we observe in TC cells in these different conditions. ..."
Reference:
1 . Rubin JE, Terman D (2004) High frequency stimulation of the subthalamic nucleus eliminates pathological thalamic rhythmicity in a computational model. J Comput Neurosci 16:211-35 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network; Neuron or other electrically excitable cell;
Brain Region(s)/Organism: Thalamus; Basal ganglia; Subthalamic Nucleus;
Cell Type(s): Thalamus geniculate nucleus/lateral principal neuron; Subthalamus nucleus projection neuron; Globus pallidus neuron;
Channel(s): I Na,t; I T low threshold; I K; I Calcium;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: XPP;
Model Concept(s): Activity Patterns; Parkinson's; Deep brain stimulation;
Implementer(s):
Search NeuronDB for information about:  Thalamus geniculate nucleus/lateral principal neuron; I Na,t; I T low threshold; I K; I Calcium;
## Set file for test2dbsfcn.ode on Sat Aug  2 15:13:26 2003
336   Number of equations and auxiliaries
133   Number of parameters
# Numerical stuff
1    nout
40    nullcline mesh
8   Qual RK
1250  total
0.1  DeltaT
0  T0
0  Transient
50000  Bound
1e-12  DtMin
1  DtMax
0.001  Tolerance
0.001  Abs. Tolerance
0  Max Delay
100   Eigenvector iterates
0.001  Eigenvector tolerance
0.001  Newton tolerance
0  Poincare plane
1e-05  Boundary value tolerance
1e-05  Boundary value epsilon
20   Boundary value iterates
0   Poincare None
1   Poincare variable
1   Poincare sign
0   Stop on Section
0   Delay flag
0  Current time
0  Last Time
0   MyStart
0   INFLAG
# Delays
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
# Bndry conds
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
# Old ICs
-69.12811098270906  V1
-69.12717763396287  V2
-56.97710097612938  V3
-60.40240307615291  V4
-70.35685487158057  V5
-68.04097200574483  V6
-62.71441018987391  V7
-62.35273706279698  V8
-61.2315138249805  V9
-62.49025842554087  V10
-69.62193980245195  V11
-68.96718916171575  V12
-61.27761451007183  V13
-63.15653567004609  V14
-68.84549539802367  V15
-69.62849202858642  V16
0.281676673132536  H1
0.281390198092845  H2
0.1147224796380693  H3
0.1013289788582662  H4
0.2691768417469897  H5
0.2669279539353307  H6
0.09389971242310781  H7
0.09473859460937538  H8
0.0981880124132382  H9
0.09507382778249453  H10
0.3034863026177871  H11
0.2867011121924212  H12
0.09873570268068801  H13
0.0920559289897373  H14
0.3712568369348424  H15
0.2914551076847784  H16
0.05216349524028844  N1
0.0522201660137626  N2
0.3457226452235606  N3
0.3981199126105596  N4
0.05439525881862976  N5
0.07509678542136203  N6
0.4358811063874363  N7
0.4309843836844034  N8
0.4156227363654603  N9
0.4328154047590343  N10
0.04860954618978539  N11
0.05362014676979376  N12
0.4132918529573688  N13
0.4502171925294871  N14
0.03695776167753602  N15
0.04779848200791642  N16
0.36575531780011  R1
0.3656741638734325  R2
0.1047063811532236  R3
0.05999002684181447  R4
0.4194505094909791  R5
0.3550002990863161  R6
0.06559762689099767  R7
0.06597829319364709  R8
0.05695024808284845  R9
0.05674833587777332  R10
0.2711077844193491  R11
0.2850342711476045  R12
0.05378337783814302  R13
0.06013298744883041  R14
0.2746483360857667  R15
0.2624220317855703  R16
0.2228356374439819  CA1
0.2215569285793844  CA2
0.2962883874138816  CA3
0.3276944759695554  CA4
0.2939980512025778  CA5
0.3208947959028472  CA6
0.2232225254924457  CA7
0.2095096020136324  CA8
0.2289913049589802  CA9
0.2639352805340913  CA10
0.2469446827702531  CA11
0.2833206038455554  CA12
0.3050966085138562  CA13
0.3082744801459322  CA14
0.2039989465586884  CA15
0.246911851514305  CA16
0.002644073084260944  S1
0.002644391634392721  S2
0.01150797503546001  S3
0.008168033091150667  S4
0.002409921890862083  S5
0.00291084685070629  S6
0.04911381127204784  S7
0.0300452593169284  S8
0.01043022380858413  S9
0.03456253311546535  S10
0.002558367871739318  S11
0.002638159168738268  S12
0.01008366987484816  S13
0.2388454131193506  S14
0.002678519783062835  S15
0.00254962123117104  S16
0.01800015438740976  SI1
0.01800694380308231  SI2
0.4047074514140117  SI3
0.5493493021160042  SI4
0.01768949008318783  SI5
0.02713611723761502  SI6
0.6716028216189945  SI7
0.6523636537411861  SI8
0.5942330304232434  SI9
0.6574015263277148  SI10
0.02650859805415025  SI11
0.02661376487730824  SI12
0.592579255053132  SI13
0.7297887374392803  SI14
0.02567391682567561  SI15
0.02561118806721522  SI16
-74.80387067121256  VG1
-79.5933546345955  VG2
-75.10830345522299  VG3
-73.14431518147096  VG4
-74.78093122364646  VG5
-79.5936246493865  VG6
-75.0465776200093  VG7
-71.94802990173955  VG8
-73.12505041116347  VG9
-73.33942678168052  VG10
-78.58263497396057  VG11
-77.3436804599869  VG12
-75.38872451219477  VG13
-77.14231198118823  VG14
-77.46020045962824  VG15
-69.9908680793304  VG16
0.265637111102081  NG1
0.5965927493794086  NG2
0.1518792436218246  NG3
0.1593398284607819  NG4
0.265610123198105  NG5
0.5966397550443858  NG6
0.1527359779576022  NG7
0.1472996663779491  NG8
0.1626166192073982  NG9
0.1565779802759985  NG10
0.2934156566790282  NG11
0.5251653039713522  NG12
0.1550376609079995  NG13
0.1334418495247361  NG14
0.3264726655263877  NG15
0.3588338974278034  NG16
0.5646497233238033  HG1
0.3336732679621591  HG2
0.8004278076154859  HG3
0.795996694161231  HG4
0.5646750341027785  HG5
0.3336520523566163  HG6
0.799171619106732  HG7
0.8112550161172151  HG8
0.7913316688787225  HG9
0.7999982196752538  HG10
0.5408706956806638  HG11
0.380775784776796  HG12
0.7909635752807689  HG13
0.8134626100927261  HG14
0.5121332279008678  HG15
0.4853485513240676  HG16
0.8101947256348955  RG1
0.776948971508445  RG2
0.9278270734170959  RG3
0.9174616344319636  RG4
0.8092597031797867  RG5
0.7765293857140189  RG6
0.9264232411850234  RG7
0.9382953837301133  RG8
0.9099203451743211  RG9
0.9292055440668339  RG10
0.8528017665089627  RG11
0.8338190781348963  RG12
0.8994807128785964  RG13
0.9286127195280354  RG14
0.8407639891053467  RG15
0.8163555753368956  RG16
0.1447958067443786  CAG1
0.1325758888634619  CAG2
0.1389009147956837  CAG3
0.128588039855197  CAG4
0.1431961092510587  CAG5
0.1334101964027299  CAG6
0.134814165664149  CAG7
0.1396682654966533  CAG8
0.1331644797565681  CAG9
0.1617392143914989  CAG10
0.1410568935524672  CAG11
0.1300281140399894  CAG12
0.1326887867423998  CAG13
0.1528950320228091  CAG14
0.1440662185684193  CAG15
0.1238347298630465  CAG16
0.7549064846063573  SG1
0.9047755660974316  SG2
0.007622877289237556  SG3
0.0176641239537329  SG4
0.7547877372519034  SG5
0.9047887200864839  SG6
0.007632081028952417  SG7
0.01769144991225911  SG8
0.006547384460456033  SG9
0.004151071485089863  SG10
0.7895202726223848  SG11
0.8735162410906138  SG12
0.006545585340253204  SG13
0.004148078191404642  SG14
0.8026147527277904  SG15
0.8054989662144936  SG16
0.4203091199521011  SGGI1
0.6693490508275813  SGGI2
4.393110121391565e-06  SGGI3
3.519507195504246e-05  SGGI4
0.4201537832177801  SGGI5
0.6693774810896055  SGGI6
4.409822083517933e-06  SGGI7
3.54308033205531e-05  SGGI8
3.103020259753987e-06  SGGI9
1.025602346199523e-06  SGGI10
0.4712291110699801  SGGI11
0.6061977420195918  SGGI12
3.112810699328764e-06  SGGI13
1.001706067049635e-06  SGGI14
0.4907986013509174  SGGI15
0.4947031090826632  SGGI16
-85.71294578085116  VGI1
-89.25365728340867  VGI2
-77.88412953379547  VGI3
-66.89831172646826  VGI4
-85.67924606255713  VGI5
-88.98683288877176  VGI6
-53.36099261592358  VGI7
-72.59192691357568  VGI8
-72.89501867509094  VGI9
10.73282746213917  VGI10
-86.50085659033785  VGI11
-88.02586763518649  VGI12
-62.92906387287897  VGI13
-78.18122232323546  VGI14
-84.95285981031955  VGI15
-86.93966607201804  VGI16
0.08191257161198202  NGI1
0.07036633372246581  NGI2
0.6897159868014461  NGI3
0.3970849111215521  NGI4
0.08186454731080141  NGI5
0.07186996861009079  NGI6
0.3745988460453169  NGI7
0.5087029210697127  NGI8
0.5046083496991685  NGI9
0.8317284424837751  NGI10
0.1058531332620342  NGI11
0.1443602610809459  NGI12
0.371325009134022  NGI13
0.8097917218004893  NGI14
0.2836531345386089  NGI15
0.1066665386787778  NGI16
0.8661095197428768  HGI1
0.9028922934462152  HGI2
0.2663069119390469  HGI3
0.4350118708639981  HGI4
0.8663760245198685  HGI5
0.9005726131990338  HGI6
0.4589532548896784  HGI7
0.3652089633224476  HGI8
0.3648014209381595  HGI9
0.2356346457271996  HGI10
0.8045209509867527  HGI11
0.764617259699407  HGI12
0.4574388928021531  HGI13
0.2132693765353696  HGI14
0.5503758289856476  HGI15
0.8081448282807453  HGI16
0.8299335246503925  RGI1
0.9166774842568947  RGI2
0.4223187708498386  RGI3
0.4775717006628248  RGI4
0.8280939361636874  RGI5
0.9011828271734866  RGI6
0.4267660096419663  RGI7
0.4895324814531421  RGI8
0.4526750260796286  RGI9
0.4127822484137393  RGI10
0.6631376511929418  RGI11
0.635608087264053  RGI12
0.4590995708366686  RGI13
0.4273802721929745  RGI14
0.6319027616874602  RGI15
0.6760987714014427  RGI16
0.2362025622457039  CAGI1
0.2399045565679337  CAGI2
0.2486303704241155  CAGI3
0.2506019026595908  CAGI4
0.2423033492147151  CAGI5
0.2540310663151549  CAGI6
0.2291180611422414  CAGI7
0.2252025022951157  CAGI8
0.2464722687126124  CAGI9
0.2453491212783939  CAGI10
0.2300063596636656  CAGI11
0.2519495117364803  CAGI12
0.2631584767380465  CAGI13
0.2526884444384774  CAGI14
0.2339812027396742  CAGI15
0.2286575089272404  CAGI16
0.001346085366589612  SGI1
0.0001737250905970536  SGI2
0.8608756060841254  SGI3
0.7031025010198553  SGI4
0.001076354353086035  SGI5
0.001660042860242276  SGI6
0.6292586362228153  SGI7
0.7820907629776708  SGI8
0.7795982848387322  SGI9
0.8804031544461797  SGI10
0.04665689491279193  SGI11
0.06382034599086839  SGI12
0.6672442595451157  SGI13
0.9034916153990793  SGI14
0.6352237744660592  SGI15
0.05637476682614843  SGI16
-79.31646988450696  VT1
0.1095674472253051  HTT1
0.9999226323835648  HT1
-80.52420295721285  VT2
0.1125995252257795  HTT2
0.9999432568069376  HT2
0  COUN1
0  COUN2
50  Z1
50  Z2
0  GOOD1
0  GOOD2
# Ending  ICs
-63.02561821174685  V1
-63.22236109295277  V2
-69.3184130317683  V3
-70.33210019492124  V4
-55.91818580878974  V5
-61.67644045880005  V6
-69.69915029226914  V7
-69.80110092376265  V8
-68.80685941352412  V9
-68.84124870182687  V10
-61.33588957987941  V11
-61.28655575271043  V12
-69.10437169607052  V13
-68.76186715164563  V14
-62.44109554692372  V15
-61.07742120056528  V16
0.09236395408604231  H1
0.09177457008118874  H2
0.1654216901761081  H3
0.3237705216448455  H4
0.1191232744190006  H5
0.09686024018052675  H6
0.2108382594045349  H7
0.2059622904058425  H8
0.4283525312346126  H9
0.4181906302429899  H10
0.09861124185086752  H11
0.09900367512148547  H12
0.3779711530280475  H13
0.3823005045296959  H14
0.09544353075790456  H15
0.09915854638619065  H16
0.4142218564359074  N1
0.4171691121932414  N2
0.1909194596523165  N3
0.1042220295986705  N4
0.2887656564053019  N5
0.4021960490588918  N6
0.1244604205755661  N7
0.1266607341988308  N8
0.02888501723854752  N9
0.02868160409895836  N10
0.4162993273397929  N11
0.4142882799429862  N12
0.03186144817981399  N13
0.03177378528826483  N14
0.4335577835045417  N15
0.4135063759481261  N16
0.1901692726225  R1
0.1888743804096129  R2
0.2637046352610132  R3
0.3272666991395335  R4
0.2471564284812515  R5
0.1607541695606363  R6
0.2550152079843789  R7
0.2550820330075756  R8
0.2802677952417053  R9
0.2825448306733358  R10
0.03472611962384477  R11
0.03271595239384381  R12
0.2871535206165031  R13
0.2954482791041  R14
0.03042743183349372  R15
0.03272080679861695  R16
0.250561473184899  CA1
0.2529896534173813  CA2
0.3177093558393729  CA3
0.3443524175364035  CA4
0.3746948289235348  CA5
0.2990774013111289  CA6
0.2438102886247215  CA7
0.2460717483975722  CA8
0.2103199206008577  CA9
0.2314593437137446  CA10
0.2796841286551143  CA11
0.3094900949344483  CA12
0.2739823949242481  CA13
0.3073557498648748  CA14
0.2581647090341286  CA15
0.2571281364735475  CA16
0.03293769451784053  S1
0.04500903292300461  S2
0.002540857116947778  S3
0.002346556459179683  S4
0.01330851907858371  S5
0.009846043047751989  S6
0.002443484320157853  S7
0.00243852248011373  S8
0.002693601339092571  S9
0.002682064895341803  S10
0.01063498569578925  S11
0.01011533539768199  S12
0.002624565869195628  S13
0.002693683787101085  S14
0.0327280032092451  S15
0.009573932660346078  S16
0.6709320819573326  SI1
0.6837879906674713  SI2
0.09788535782164799  SI3
0.03455121985368396  SI4
0.3398241067217014  SI5
0.6043275675026154  SI6
0.04444004569608638  SI7
0.04542800318538005  SI8
0.02732970142853463  SI9
0.02704157282898927  SI10
0.5960235007610735  SI11
0.5920547610720707  SI12
0.02730690675786756  SI13
0.02642388311691874  SI14
0.6540414827399319  SI15
0.5845827878954124  SI16
-74.54702730745265  VG1
-73.59557641806853  VG2
-79.44423764445187  VG3
-76.67074738616449  VG4
-74.25223450068874  VG5
-75.03538106878365  VG6
-79.5056111569025  VG7
-77.35238024410889  VG8
-78.6073267108424  VG9
-77.80202027509482  VG10
-71.85810859651606  VG11
-71.4169026595293  VG12
-78.59413221895993  VG13
-77.79217916015861  VG14
-76.49180460287899  VG15
-73.43379135657571  VG16
0.1575105146103429  NG1
0.1572012940945493  NG2
0.5071646385443419  NG3
0.3227570562405988  NG4
0.1595342470744336  NG5
0.1546444358110607  NG6
0.55725097091834  NG7
0.3189235561251694  NG8
0.2863672804991405  NG9
0.3368128115686143  NG10
0.1738539212842649  NG11
0.1737961380496443  NG12
0.2861396785228235  NG13
0.3363983486558816  NG14
0.1378304102809951  NG15
0.171956805104213  NG16
0.7921368550581206  HG1
0.7977902663888276  HG2
0.3877771429259619  HG3
0.5143814480193799  HG4
0.7903321378026048  HG5
0.795986644349623  HG6
0.3596178439172846  HG7
0.517989325655664  HG8
0.5471646358451251  HG9
0.5043612997429548  HG10
0.7784964716704993  HG11
0.7791705853189447  HG12
0.5473446639440214  HG13
0.5047286958036268  HG14
0.8068647361306944  HG15
0.7738945763189974  HG16
0.8897926521337118  RG1
0.8969760333513215  RG2
0.8416474751717132  RG3
0.8697736876296126  RG4
0.8871794596494079  RG5
0.895260812534469  RG6
0.8407425008587452  RG7
0.870620790627048  RG8
0.85716245021237  RG9
0.8365605396778097  RG10
0.8935164814557947  RG11
0.8936849852578129  RG12
0.8489707800735063  RG13
0.8163601562290789  RG14
0.9091199051041607  RG15
0.8722779968147727  RG16
0.112304892087555  CAG1
0.1469910282682155  CAG2
0.168701646562669  CAG3
0.151530300625628  CAG4
0.1079284895041885  CAG5
0.1478871353385584  CAG6
0.1755691964175919  CAG7
0.1524049979108612  CAG8
0.1614408687352443  CAG9
0.1444128036325096  CAG10
0.1091508522155669  CAG11
0.1038126352870794  CAG12
0.1553718904475566  CAG13
0.1444708861168866  CAG14
0.1100330669452639  CAG15
0.1043546691288789  CAG16
0.07051666146677144  SG1
0.1643830033489514  SG2
0.8803221208438273  SG3
0.7955386642295023  SG4
0.06967855066929424  SG5
0.1644183560214746  SG6
0.8930917456172841  SG7
0.7954654805890912  SG8
0.78604384823147  SG9
0.8087389615316  SG10
0.0003692782123582206  SG11
0.0004393280139499427  SG12
0.7858782899745692  SG13
0.8084312252337441  SG14
0.0003692723720557424  SG15
0.0004357412023903713  SG16
0.001133430376739429  SGGI1
0.009301623003177257  SGGI2
0.621949219063838  SGGI3
0.4821007232939579  SGGI4
0.001101140700462458  SGGI5
0.009308971347523903  SGGI6
0.6443887957724811  SGGI7
0.4820189052209531  SGGI8
0.4660891087632231  SGGI9
0.4999663651803494  SGGI10
1.979029863956799e-07  SGGI11
1.950098341373152e-07  SGGI12
0.4658226414884322  SGGI13
0.4994455251094082  SGGI14
8.009272481440118e-08  SGGI15
2.37253551787507e-07  SGGI16
-63.20547694479664  VGI1
-71.32158066167312  VGI2
-82.94000214978436  VGI3
-82.41091471335355  VGI4
-78.42431956444496  VGI5
-58.36265461032063  VGI6
-84.48665605730065  VGI7
-82.25118197576366  VGI8
-82.55386745489491  VGI9
-83.04094868773448  VGI10
-58.6434371997331  VGI11
-77.5093356285725  VGI12
-82.48969822723069  VGI13
-82.95105569664337  VGI14
-73.54042911473645  VGI15
37.84579803391364  VGI16
0.3834250469500803  NGI1
0.4837569692573944  NGI2
0.1180913561909271  NGI3
0.1641191419801904  NGI4
0.7303410579118836  NGI5
0.3616724516108063  NGI6
0.1086267125180665  NGI7
0.1283846366418645  NGI8
0.1334135851388631  NGI9
0.1605964567566105  NGI10
0.3595384557116194  NGI11
0.7019514115329483  NGI12
0.1405607297318956  NGI13
0.1708741103923462  NGI14
0.5307948562130356  NGI15
0.6761812978974358  NGI16
0.4508280665365476  HGI1
0.3848822801090405  HGI2
0.8055130918839301  HGI3
0.6699094757646802  HGI4
0.2480972394737395  HGI5
0.4697495275242991  HGI6
0.8310432825199181  HGI7
0.7293570067562897  HGI8
0.7198177692291279  HGI9
0.677256454901299  HGI10
0.4711366523582273  HGI11
0.2581006365542166  HGI12
0.7080291577068197  HGI13
0.6636356077572623  HGI14
0.3531203813085761  HGI15
0.3334295670829293  HGI16
0.5085981681083616  RGI1
0.5199678498551277  RGI2
0.6491322229622383  RGI3
0.5694029464477712  RGI4
0.4968317754415598  RGI5
0.4810170348326543  RGI6
0.7126667506150736  RGI7
0.584209888620024  RGI8
0.6919512889194444  RGI9
0.6629761636149646  RGI10
0.4181926076123568  RGI11
0.4042933235186801  RGI12
0.6780576640361073  RGI13
0.645852615156695  RGI14
0.4514203534524424  RGI15
0.4045983209299936  RGI16
0.2663231609724284  CAGI1
0.2425279337879178  CAGI2
0.2468544010893826  CAGI3
0.2426173101793528  CAGI4
0.278297980056281  CAGI5
0.2434269395459023  CAGI6
0.2174011370886436  CAGI7
0.2337154634615383  CAGI8
0.2120256925149332  CAGI9
0.224275850269073  CAGI10
0.282319017956355  CAGI11
0.2850179863452723  CAGI12
0.2205622207534466  CAGI13
0.2330374792932023  CAGI14
0.2832901867102043  CAGI15
0.2693448458187996  CAGI16
0.6753282436805498  SGI1
0.7647280127562034  SGI2
0.2038867027691741  SGI3
0.5080688897613265  SGI4
0.8752221493584865  SGI5
0.6284923868052706  SGI6
0.1569032149703416  SGI7
0.4095013487769359  SGI8
0.4411786954563598  SGI9
0.4970111031955169  SGI10
0.6249552915500777  SGI11
0.8661356493252917  SGI12
0.4583159660314014  SGI13
0.5132758083165239  SGI14
0.7935737337258718  SGI15
0.7888677384671987  SGI16
-81.61515370148244  VT1
0.08656244907145687  HTT1
0.9999435566528442  HT1
-81.21429679281508  VT2
0.1367490425467387  HTT2
0.9999504848673988  HT2
2  COUN1
0  COUN2
50.00937131807729  Z1
50.01263145081623  Z2
1  GOOD1
0  GOOD2
# Parameters
-60  vl
55  vna
-80  vk
30  thetam
15  sm
2.25  gl
37.5  gna
45  gk
1  tn
0.05  th
9  gahp
0.5  gca
140  vca
15  k1
5e-05  eps
22.5  kca
39  thetas
8  ss
1  xp
0  i
-39  thetah
3.1  sh
-32  thetan
-8  sn
1  taun0
100  taun1
80  thn
26  sigman
1  tauh0
500  tauh1
57  thh
3  sigmah
0.75  phi
-63  thetat
-7.8  kt
0.5  gt
0.5  phir
-67  thetar
2  kr
7.1  taur0
17.5  taur1
-68  thr
2.2  sigmar
5  alpha
1  beta
-30  ab
0.9  gsyn
-100  vsyn
0.25  rth
-0.07000000000000001  rsig
0.5  rho1
0.9  a1
0  hets
25  inp
1  alphai
0.05  betai
120  gnag
30  gkg
30  gahpg
0.5  gtg
0.1  gcag
0.1  glg
55  vnag
-80  vkg
120  vcag
-55  vlg
-57  thetasg
2  ksg
-35  thetas1g
2  ks1g
-70  thetarg
-2  krg
30  taurg
-37  thetamg
10  sigmamg
-50  thetang
14  sigmang
0.05  taun0g
0.27  taun1g
-40  thng
-12  sng
-58  thetahg
-12  sigmahg
0.05  tauh0g
0.27  tauh1g
-40  thhg
-12  shg
30  k1g
15  kcag
0.0001  epsg
1  phig
0.1  deltang
0.05  deltahg
2  iapp
0.3  gsyngg
-80  vsyngg
0.3  gsyng
0  vsyng
2  alphag
0.04  betag
-20  abg
5  iappi
1  igl
0  thresg
0  hetg
0.08  betagi
15  kcagi
0.5  gsyngi
1  gsynggi
1  alphaggi
0.1  betaggi
-85  vsyntc
0.06  gsyntc
200  asg
0.4  bsg
5  itc
-80  shi
-95  shi2
5  dur
20  dur2
50  period
3  gnabar
5  gkbar
0.05  glbar
50  ena
-90  ek
-70  eleak
5  gtbar
5.5  qht
1  tadj
4  apr
0.3  apt
40  thr0
# Graphics
0.7071097251469876  rm
-0.4999999999826661  rm
0.4999958366025516  rm
0.7071038372138505  rm
0.5000041633974484  rm
-0.4999999999826661  rm
0  rm
0.7071038372138505  rm
0.7071097251469876  rm
0
321
2
1
0
0
1
2
1
0
0
1
2
1
0
0
1
2
1
0
0
1
2
1
0
0
1
2
1
0
0
1
2
1
0
0
1
2
1
0
0
1
2
1
0
0
1
2
1
0
-1000   
1000   
0    
0   3DFlag
1   Timeflag
0   Colorflag
0   Type
1  color scale
0   minscale
1250   xmax
0   xmin
20   ymax
-80   ymin
12   zmax
-12   zmin
500  
0.002
-30
0.02   
0   
0.08333333333333333   
45   Theta
45   Phi
0    xshft
0    yshft
0    zshft
0   xlo
-80   ylo
0   
-80   
1250   xhi
20   yhi
1250   
20   
# Transpose variables etc
V1
2   n columns
1   n rows
1   row skip
1   col skip
1   row 0
# Coupling stuff for H funs
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
# Array plot stuff

1   NCols
0   Row 1
50   NRows
8   RowSkip
0  Zmin
1  Zmax
# Torus information 
0    Torus flag 1=ON
6.283185307179586  Torus period
# Range information
vl
-1   eq-range stab col
0   shoot flag 1=on
10   eq-range steps
0  eq_range low
1  eq_range high
V1
V1
20   Range steps
0   Cycle color 1=on
1   Reset data 1=on
1   Use old I.C.s 1=yes
0  Par1 low
0  Par2 low
1  Par1 high
1  Par2 high
vl
0   BVP side
0   color cycle flag 1=on
10   BVP range steps
0  BVP range low
1  BVP range high
RHS etc ...
dV1/dT=-(IL(V1)+INA(V1,H1)+IK(V1,N1)+IAHP(V1,CA1)+ICA(V1)+IT(V1,R1))-ISYN1+HETS*1+INP+DBS(T)
dV2/dT=-(IL(V2)+INA(V2,H2)+IK(V2,N2)+IAHP(V2,CA2)+ICA(V2)+IT(V2,R2))-ISYN2+HETS*2+INP+DBS(T)
dV3/dT=-(IL(V3)+INA(V3,H3)+IK(V3,N3)+IAHP(V3,CA3)+ICA(V3)+IT(V3,R3))-ISYN3+HETS*3+INP+DBS(T)
dV4/dT=-(IL(V4)+INA(V4,H4)+IK(V4,N4)+IAHP(V4,CA4)+ICA(V4)+IT(V4,R4))-ISYN4+HETS*4+INP+DBS(T)
dV5/dT=-(IL(V5)+INA(V5,H5)+IK(V5,N5)+IAHP(V5,CA5)+ICA(V5)+IT(V5,R5))-ISYN5+HETS*5+INP+DBS(T)
dV6/dT=-(IL(V6)+INA(V6,H6)+IK(V6,N6)+IAHP(V6,CA6)+ICA(V6)+IT(V6,R6))-ISYN6+HETS*6+INP+DBS(T)
dV7/dT=-(IL(V7)+INA(V7,H7)+IK(V7,N7)+IAHP(V7,CA7)+ICA(V7)+IT(V7,R7))-ISYN7+HETS*7+INP+DBS(T)
dV8/dT=-(IL(V8)+INA(V8,H8)+IK(V8,N8)+IAHP(V8,CA8)+ICA(V8)+IT(V8,R8))-ISYN8+HETS*8+INP+DBS(T)
dV9/dT=-(IL(V9)+INA(V9,H9)+IK(V9,N9)+IAHP(V9,CA9)+ICA(V9)+IT(V9,R9))-ISYN9+HETS*9+INP+DBS(T)
dV10/dT=-(IL(V10)+INA(V10,H10)+IK(V10,N10)+IAHP(V10,CA10)+ICA(V10)+IT(V10,R10))-ISYN10+HETS*10+INP+DBS(T)
dV11/dT=-(IL(V11)+INA(V11,H11)+IK(V11,N11)+IAHP(V11,CA11)+ICA(V11)+IT(V11,R11))-ISYN11+HETS*11+INP+DBS(T)
dV12/dT=-(IL(V12)+INA(V12,H12)+IK(V12,N12)+IAHP(V12,CA12)+ICA(V12)+IT(V12,R12))-ISYN12+HETS*12+INP+DBS(T)
dV13/dT=-(IL(V13)+INA(V13,H13)+IK(V13,N13)+IAHP(V13,CA13)+ICA(V13)+IT(V13,R13))-ISYN13+HETS*13+INP+DBS(T)
dV14/dT=-(IL(V14)+INA(V14,H14)+IK(V14,N14)+IAHP(V14,CA14)+ICA(V14)+IT(V14,R14))-ISYN14+HETS*14+INP+DBS(T)
dV15/dT=-(IL(V15)+INA(V15,H15)+IK(V15,N15)+IAHP(V15,CA15)+ICA(V15)+IT(V15,R15))-ISYN15+HETS*15+INP+DBS(T)
dV16/dT=-(IL(V16)+INA(V16,H16)+IK(V16,N16)+IAHP(V16,CA16)+ICA(V16)+IT(V16,R16))-ISYN16+HETS*16+INP+DBS(T)
dH1/dT=PHI*(HINF(V1)-H1)/TAUH(V1)
dH2/dT=PHI*(HINF(V2)-H2)/TAUH(V2)
dH3/dT=PHI*(HINF(V3)-H3)/TAUH(V3)
dH4/dT=PHI*(HINF(V4)-H4)/TAUH(V4)
dH5/dT=PHI*(HINF(V5)-H5)/TAUH(V5)
dH6/dT=PHI*(HINF(V6)-H6)/TAUH(V6)
dH7/dT=PHI*(HINF(V7)-H7)/TAUH(V7)
dH8/dT=PHI*(HINF(V8)-H8)/TAUH(V8)
dH9/dT=PHI*(HINF(V9)-H9)/TAUH(V9)
dH10/dT=PHI*(HINF(V10)-H10)/TAUH(V10)
dH11/dT=PHI*(HINF(V11)-H11)/TAUH(V11)
dH12/dT=PHI*(HINF(V12)-H12)/TAUH(V12)
dH13/dT=PHI*(HINF(V13)-H13)/TAUH(V13)
dH14/dT=PHI*(HINF(V14)-H14)/TAUH(V14)
dH15/dT=PHI*(HINF(V15)-H15)/TAUH(V15)
dH16/dT=PHI*(HINF(V16)-H16)/TAUH(V16)
dN1/dT=PHI*(NINF(V1)-N1)/TAUN(V1)
dN2/dT=PHI*(NINF(V2)-N2)/TAUN(V2)
dN3/dT=PHI*(NINF(V3)-N3)/TAUN(V3)
dN4/dT=PHI*(NINF(V4)-N4)/TAUN(V4)
dN5/dT=PHI*(NINF(V5)-N5)/TAUN(V5)
dN6/dT=PHI*(NINF(V6)-N6)/TAUN(V6)
dN7/dT=PHI*(NINF(V7)-N7)/TAUN(V7)
dN8/dT=PHI*(NINF(V8)-N8)/TAUN(V8)
dN9/dT=PHI*(NINF(V9)-N9)/TAUN(V9)
dN10/dT=PHI*(NINF(V10)-N10)/TAUN(V10)
dN11/dT=PHI*(NINF(V11)-N11)/TAUN(V11)
dN12/dT=PHI*(NINF(V12)-N12)/TAUN(V12)
dN13/dT=PHI*(NINF(V13)-N13)/TAUN(V13)
dN14/dT=PHI*(NINF(V14)-N14)/TAUN(V14)
dN15/dT=PHI*(NINF(V15)-N15)/TAUN(V15)
dN16/dT=PHI*(NINF(V16)-N16)/TAUN(V16)
dR1/dT=PHIR*(RINF(V1)-R1)/TAUR(V1)
dR2/dT=PHIR*(RINF(V2)-R2)/TAUR(V2)
dR3/dT=PHIR*(RINF(V3)-R3)/TAUR(V3)
dR4/dT=PHIR*(RINF(V4)-R4)/TAUR(V4)
dR5/dT=PHIR*(RINF(V5)-R5)/TAUR(V5)
dR6/dT=PHIR*(RINF(V6)-R6)/TAUR(V6)
dR7/dT=PHIR*(RINF(V7)-R7)/TAUR(V7)
dR8/dT=PHIR*(RINF(V8)-R8)/TAUR(V8)
dR9/dT=PHIR*(RINF(V9)-R9)/TAUR(V9)
dR10/dT=PHIR*(RINF(V10)-R10)/TAUR(V10)
dR11/dT=PHIR*(RINF(V11)-R11)/TAUR(V11)
dR12/dT=PHIR*(RINF(V12)-R12)/TAUR(V12)
dR13/dT=PHIR*(RINF(V13)-R13)/TAUR(V13)
dR14/dT=PHIR*(RINF(V14)-R14)/TAUR(V14)
dR15/dT=PHIR*(RINF(V15)-R15)/TAUR(V15)
dR16/dT=PHIR*(RINF(V16)-R16)/TAUR(V16)
dCA1/dT=PHI*EPS*(-ICA(V1)-IT(V1,R1)-KCA*CA1)
dCA2/dT=PHI*EPS*(-ICA(V2)-IT(V2,R2)-KCA*CA2)
dCA3/dT=PHI*EPS*(-ICA(V3)-IT(V3,R3)-KCA*CA3)
dCA4/dT=PHI*EPS*(-ICA(V4)-IT(V4,R4)-KCA*CA4)
dCA5/dT=PHI*EPS*(-ICA(V5)-IT(V5,R5)-KCA*CA5)
dCA6/dT=PHI*EPS*(-ICA(V6)-IT(V6,R6)-KCA*CA6)
dCA7/dT=PHI*EPS*(-ICA(V7)-IT(V7,R7)-KCA*CA7)
dCA8/dT=PHI*EPS*(-ICA(V8)-IT(V8,R8)-KCA*CA8)
dCA9/dT=PHI*EPS*(-ICA(V9)-IT(V9,R9)-KCA*CA9)
dCA10/dT=PHI*EPS*(-ICA(V10)-IT(V10,R10)-KCA*CA10)
dCA11/dT=PHI*EPS*(-ICA(V11)-IT(V11,R11)-KCA*CA11)
dCA12/dT=PHI*EPS*(-ICA(V12)-IT(V12,R12)-KCA*CA12)
dCA13/dT=PHI*EPS*(-ICA(V13)-IT(V13,R13)-KCA*CA13)
dCA14/dT=PHI*EPS*(-ICA(V14)-IT(V14,R14)-KCA*CA14)
dCA15/dT=PHI*EPS*(-ICA(V15)-IT(V15,R15)-KCA*CA15)
dCA16/dT=PHI*EPS*(-ICA(V16)-IT(V16,R16)-KCA*CA16)
dS1/dT=ALPHA*(1-S1)*SINF(V1+AB)-BETA*S1
dS2/dT=ALPHA*(1-S2)*SINF(V2+AB)-BETA*S2
dS3/dT=ALPHA*(1-S3)*SINF(V3+AB)-BETA*S3
dS4/dT=ALPHA*(1-S4)*SINF(V4+AB)-BETA*S4
dS5/dT=ALPHA*(1-S5)*SINF(V5+AB)-BETA*S5
dS6/dT=ALPHA*(1-S6)*SINF(V6+AB)-BETA*S6
dS7/dT=ALPHA*(1-S7)*SINF(V7+AB)-BETA*S7
dS8/dT=ALPHA*(1-S8)*SINF(V8+AB)-BETA*S8
dS9/dT=ALPHA*(1-S9)*SINF(V9+AB)-BETA*S9
dS10/dT=ALPHA*(1-S10)*SINF(V10+AB)-BETA*S10
dS11/dT=ALPHA*(1-S11)*SINF(V11+AB)-BETA*S11
dS12/dT=ALPHA*(1-S12)*SINF(V12+AB)-BETA*S12
dS13/dT=ALPHA*(1-S13)*SINF(V13+AB)-BETA*S13
dS14/dT=ALPHA*(1-S14)*SINF(V14+AB)-BETA*S14
dS15/dT=ALPHA*(1-S15)*SINF(V15+AB)-BETA*S15
dS16/dT=ALPHA*(1-S16)*SINF(V16+AB)-BETA*S16
dSI1/dT=ALPHAI*(1-SI1)*SINF(V1+AB)-BETAI*SI1
dSI2/dT=ALPHAI*(1-SI2)*SINF(V2+AB)-BETAI*SI2
dSI3/dT=ALPHAI*(1-SI3)*SINF(V3+AB)-BETAI*SI3
dSI4/dT=ALPHAI*(1-SI4)*SINF(V4+AB)-BETAI*SI4
dSI5/dT=ALPHAI*(1-SI5)*SINF(V5+AB)-BETAI*SI5
dSI6/dT=ALPHAI*(1-SI6)*SINF(V6+AB)-BETAI*SI6
dSI7/dT=ALPHAI*(1-SI7)*SINF(V7+AB)-BETAI*SI7
dSI8/dT=ALPHAI*(1-SI8)*SINF(V8+AB)-BETAI*SI8
dSI9/dT=ALPHAI*(1-SI9)*SINF(V9+AB)-BETAI*SI9
dSI10/dT=ALPHAI*(1-SI10)*SINF(V10+AB)-BETAI*SI10
dSI11/dT=ALPHAI*(1-SI11)*SINF(V11+AB)-BETAI*SI11
dSI12/dT=ALPHAI*(1-SI12)*SINF(V12+AB)-BETAI*SI12
dSI13/dT=ALPHAI*(1-SI13)*SINF(V13+AB)-BETAI*SI13
dSI14/dT=ALPHAI*(1-SI14)*SINF(V14+AB)-BETAI*SI14
dSI15/dT=ALPHAI*(1-SI15)*SINF(V15+AB)-BETAI*SI15
dSI16/dT=ALPHAI*(1-SI16)*SINF(V16+AB)-BETAI*SI16
dVG1/dT=-(ITG(VG1,RG1)+INAG(VG1,HG1)+IKG(VG1,NG1)+IAHPG(VG1,CAG1)+ICAG(VG1)+ILG(VG1))+IAPP-ISYNGG1-ISYNG1+HETG*1
dVG2/dT=-(ITG(VG2,RG2)+INAG(VG2,HG2)+IKG(VG2,NG2)+IAHPG(VG2,CAG2)+ICAG(VG2)+ILG(VG2))+IAPP-ISYNGG2-ISYNG2+HETG*2
dVG3/dT=-(ITG(VG3,RG3)+INAG(VG3,HG3)+IKG(VG3,NG3)+IAHPG(VG3,CAG3)+ICAG(VG3)+ILG(VG3))+IAPP-ISYNGG3-ISYNG3+HETG*3
dVG4/dT=-(ITG(VG4,RG4)+INAG(VG4,HG4)+IKG(VG4,NG4)+IAHPG(VG4,CAG4)+ICAG(VG4)+ILG(VG4))+IAPP-ISYNGG4-ISYNG4+HETG*4
dVG5/dT=-(ITG(VG5,RG5)+INAG(VG5,HG5)+IKG(VG5,NG5)+IAHPG(VG5,CAG5)+ICAG(VG5)+ILG(VG5))+IAPP-ISYNGG5-ISYNG5+HETG*5
dVG6/dT=-(ITG(VG6,RG6)+INAG(VG6,HG6)+IKG(VG6,NG6)+IAHPG(VG6,CAG6)+ICAG(VG6)+ILG(VG6))+IAPP-ISYNGG6-ISYNG6+HETG*6
dVG7/dT=-(ITG(VG7,RG7)+INAG(VG7,HG7)+IKG(VG7,NG7)+IAHPG(VG7,CAG7)+ICAG(VG7)+ILG(VG7))+IAPP-ISYNGG7-ISYNG7+HETG*7
dVG8/dT=-(ITG(VG8,RG8)+INAG(VG8,HG8)+IKG(VG8,NG8)+IAHPG(VG8,CAG8)+ICAG(VG8)+ILG(VG8))+IAPP-ISYNGG8-ISYNG8+HETG*8
dVG9/dT=-(ITG(VG9,RG9)+INAG(VG9,HG9)+IKG(VG9,NG9)+IAHPG(VG9,CAG9)+ICAG(VG9)+ILG(VG9))+IAPP-ISYNGG9-ISYNG9+HETG*9
dVG10/dT=-(ITG(VG10,RG10)+INAG(VG10,HG10)+IKG(VG10,NG10)+IAHPG(VG10,CAG10)+ICAG(VG10)+ILG(VG10))+IAPP-ISYNGG10-ISYNG10+HETG*10
dVG11/dT=-(ITG(VG11,RG11)+INAG(VG11,HG11)+IKG(VG11,NG11)+IAHPG(VG11,CAG11)+ICAG(VG11)+ILG(VG11))+IAPP-ISYNGG11-ISYNG11+HETG*11
dVG12/dT=-(ITG(VG12,RG12)+INAG(VG12,HG12)+IKG(VG12,NG12)+IAHPG(VG12,CAG12)+ICAG(VG12)+ILG(VG12))+IAPP-ISYNGG12-ISYNG12+HETG*12
dVG13/dT=-(ITG(VG13,RG13)+INAG(VG13,HG13)+IKG(VG13,NG13)+IAHPG(VG13,CAG13)+ICAG(VG13)+ILG(VG13))+IAPP-ISYNGG13-ISYNG13+HETG*13
dVG14/dT=-(ITG(VG14,RG14)+INAG(VG14,HG14)+IKG(VG14,NG14)+IAHPG(VG14,CAG14)+ICAG(VG14)+ILG(VG14))+IAPP-ISYNGG14-ISYNG14+HETG*14
dVG15/dT=-(ITG(VG15,RG15)+INAG(VG15,HG15)+IKG(VG15,NG15)+IAHPG(VG15,CAG15)+ICAG(VG15)+ILG(VG15))+IAPP-ISYNGG15-ISYNG15+HETG*15
dVG16/dT=-(ITG(VG16,RG16)+INAG(VG16,HG16)+IKG(VG16,NG16)+IAHPG(VG16,CAG16)+ICAG(VG16)+ILG(VG16))+IAPP-ISYNGG16-ISYNG16+HETG*16
dNG1/dT=DELTANG*(NINFG(VG1)-NG1)/TAUNG(VG1)
dNG2/dT=DELTANG*(NINFG(VG2)-NG2)/TAUNG(VG2)
dNG3/dT=DELTANG*(NINFG(VG3)-NG3)/TAUNG(VG3)
dNG4/dT=DELTANG*(NINFG(VG4)-NG4)/TAUNG(VG4)
dNG5/dT=DELTANG*(NINFG(VG5)-NG5)/TAUNG(VG5)
dNG6/dT=DELTANG*(NINFG(VG6)-NG6)/TAUNG(VG6)
dNG7/dT=DELTANG*(NINFG(VG7)-NG7)/TAUNG(VG7)
dNG8/dT=DELTANG*(NINFG(VG8)-NG8)/TAUNG(VG8)
dNG9/dT=DELTANG*(NINFG(VG9)-NG9)/TAUNG(VG9)
dNG10/dT=DELTANG*(NINFG(VG10)-NG10)/TAUNG(VG10)
dNG11/dT=DELTANG*(NINFG(VG11)-NG11)/TAUNG(VG11)
dNG12/dT=DELTANG*(NINFG(VG12)-NG12)/TAUNG(VG12)
dNG13/dT=DELTANG*(NINFG(VG13)-NG13)/TAUNG(VG13)
dNG14/dT=DELTANG*(NINFG(VG14)-NG14)/TAUNG(VG14)
dNG15/dT=DELTANG*(NINFG(VG15)-NG15)/TAUNG(VG15)
dNG16/dT=DELTANG*(NINFG(VG16)-NG16)/TAUNG(VG16)
dHG1/dT=DELTAHG*(HINFG(VG1)-HG1)/TAUHG(VG1)
dHG2/dT=DELTAHG*(HINFG(VG2)-HG2)/TAUHG(VG2)
dHG3/dT=DELTAHG*(HINFG(VG3)-HG3)/TAUHG(VG3)
dHG4/dT=DELTAHG*(HINFG(VG4)-HG4)/TAUHG(VG4)
dHG5/dT=DELTAHG*(HINFG(VG5)-HG5)/TAUHG(VG5)
dHG6/dT=DELTAHG*(HINFG(VG6)-HG6)/TAUHG(VG6)
dHG7/dT=DELTAHG*(HINFG(VG7)-HG7)/TAUHG(VG7)
dHG8/dT=DELTAHG*(HINFG(VG8)-HG8)/TAUHG(VG8)
dHG9/dT=DELTAHG*(HINFG(VG9)-HG9)/TAUHG(VG9)
dHG10/dT=DELTAHG*(HINFG(VG10)-HG10)/TAUHG(VG10)
dHG11/dT=DELTAHG*(HINFG(VG11)-HG11)/TAUHG(VG11)
dHG12/dT=DELTAHG*(HINFG(VG12)-HG12)/TAUHG(VG12)
dHG13/dT=DELTAHG*(HINFG(VG13)-HG13)/TAUHG(VG13)
dHG14/dT=DELTAHG*(HINFG(VG14)-HG14)/TAUHG(VG14)
dHG15/dT=DELTAHG*(HINFG(VG15)-HG15)/TAUHG(VG15)
dHG16/dT=DELTAHG*(HINFG(VG16)-HG16)/TAUHG(VG16)
dRG1/dT=PHIG*(RINFG(VG1)-RG1)/TAURG
dRG2/dT=PHIG*(RINFG(VG2)-RG2)/TAURG
dRG3/dT=PHIG*(RINFG(VG3)-RG3)/TAURG
dRG4/dT=PHIG*(RINFG(VG4)-RG4)/TAURG
dRG5/dT=PHIG*(RINFG(VG5)-RG5)/TAURG
dRG6/dT=PHIG*(RINFG(VG6)-RG6)/TAURG
dRG7/dT=PHIG*(RINFG(VG7)-RG7)/TAURG
dRG8/dT=PHIG*(RINFG(VG8)-RG8)/TAURG
dRG9/dT=PHIG*(RINFG(VG9)-RG9)/TAURG
dRG10/dT=PHIG*(RINFG(VG10)-RG10)/TAURG
dRG11/dT=PHIG*(RINFG(VG11)-RG11)/TAURG
dRG12/dT=PHIG*(RINFG(VG12)-RG12)/TAURG
dRG13/dT=PHIG*(RINFG(VG13)-RG13)/TAURG
dRG14/dT=PHIG*(RINFG(VG14)-RG14)/TAURG
dRG15/dT=PHIG*(RINFG(VG15)-RG15)/TAURG
dRG16/dT=PHIG*(RINFG(VG16)-RG16)/TAURG
dCAG1/dT=EPSG*(-ICAG(VG1)-ITG(VG1,RG1)-KCAG*CAG1)
dCAG2/dT=EPSG*(-ICAG(VG2)-ITG(VG2,RG2)-KCAG*CAG2)
dCAG3/dT=EPSG*(-ICAG(VG3)-ITG(VG3,RG3)-KCAG*CAG3)
dCAG4/dT=EPSG*(-ICAG(VG4)-ITG(VG4,RG4)-KCAG*CAG4)
dCAG5/dT=EPSG*(-ICAG(VG5)-ITG(VG5,RG5)-KCAG*CAG5)
dCAG6/dT=EPSG*(-ICAG(VG6)-ITG(VG6,RG6)-KCAG*CAG6)
dCAG7/dT=EPSG*(-ICAG(VG7)-ITG(VG7,RG7)-KCAG*CAG7)
dCAG8/dT=EPSG*(-ICAG(VG8)-ITG(VG8,RG8)-KCAG*CAG8)
dCAG9/dT=EPSG*(-ICAG(VG9)-ITG(VG9,RG9)-KCAG*CAG9)
dCAG10/dT=EPSG*(-ICAG(VG10)-ITG(VG10,RG10)-KCAG*CAG10)
dCAG11/dT=EPSG*(-ICAG(VG11)-ITG(VG11,RG11)-KCAG*CAG11)
dCAG12/dT=EPSG*(-ICAG(VG12)-ITG(VG12,RG12)-KCAG*CAG12)
dCAG13/dT=EPSG*(-ICAG(VG13)-ITG(VG13,RG13)-KCAG*CAG13)
dCAG14/dT=EPSG*(-ICAG(VG14)-ITG(VG14,RG14)-KCAG*CAG14)
dCAG15/dT=EPSG*(-ICAG(VG15)-ITG(VG15,RG15)-KCAG*CAG15)
dCAG16/dT=EPSG*(-ICAG(VG16)-ITG(VG16,RG16)-KCAG*CAG16)
dSG1/dT=ALPHAG*(1-SG1)*SINFG(VG1+ABG)-BETAG*SG1
dSG2/dT=ALPHAG*(1-SG2)*SINFG(VG2+ABG)-BETAG*SG2
dSG3/dT=ALPHAG*(1-SG3)*SINFG(VG3+ABG)-BETAG*SG3
dSG4/dT=ALPHAG*(1-SG4)*SINFG(VG4+ABG)-BETAG*SG4
dSG5/dT=ALPHAG*(1-SG5)*SINFG(VG5+ABG)-BETAG*SG5
dSG6/dT=ALPHAG*(1-SG6)*SINFG(VG6+ABG)-BETAG*SG6
dSG7/dT=ALPHAG*(1-SG7)*SINFG(VG7+ABG)-BETAG*SG7
dSG8/dT=ALPHAG*(1-SG8)*SINFG(VG8+ABG)-BETAG*SG8
dSG9/dT=ALPHAG*(1-SG9)*SINFG(VG9+ABG)-BETAG*SG9
dSG10/dT=ALPHAG*(1-SG10)*SINFG(VG10+ABG)-BETAG*SG10
dSG11/dT=ALPHAG*(1-SG11)*SINFG(VG11+ABG)-BETAG*SG11
dSG12/dT=ALPHAG*(1-SG12)*SINFG(VG12+ABG)-BETAG*SG12
dSG13/dT=ALPHAG*(1-SG13)*SINFG(VG13+ABG)-BETAG*SG13
dSG14/dT=ALPHAG*(1-SG14)*SINFG(VG14+ABG)-BETAG*SG14
dSG15/dT=ALPHAG*(1-SG15)*SINFG(VG15+ABG)-BETAG*SG15
dSG16/dT=ALPHAG*(1-SG16)*SINFG(VG16+ABG)-BETAG*SG16
dSGGI1/dT=ALPHAGGI*(1-SGGI1)*SINFG(VG1+ABG)-BETAGGI*SGGI1
dSGGI2/dT=ALPHAGGI*(1-SGGI2)*SINFG(VG2+ABG)-BETAGGI*SGGI2
dSGGI3/dT=ALPHAGGI*(1-SGGI3)*SINFG(VG3+ABG)-BETAGGI*SGGI3
dSGGI4/dT=ALPHAGGI*(1-SGGI4)*SINFG(VG4+ABG)-BETAGGI*SGGI4
dSGGI5/dT=ALPHAGGI*(1-SGGI5)*SINFG(VG5+ABG)-BETAGGI*SGGI5
dSGGI6/dT=ALPHAGGI*(1-SGGI6)*SINFG(VG6+ABG)-BETAGGI*SGGI6
dSGGI7/dT=ALPHAGGI*(1-SGGI7)*SINFG(VG7+ABG)-BETAGGI*SGGI7
dSGGI8/dT=ALPHAGGI*(1-SGGI8)*SINFG(VG8+ABG)-BETAGGI*SGGI8
dSGGI9/dT=ALPHAGGI*(1-SGGI9)*SINFG(VG9+ABG)-BETAGGI*SGGI9
dSGGI10/dT=ALPHAGGI*(1-SGGI10)*SINFG(VG10+ABG)-BETAGGI*SGGI10
dSGGI11/dT=ALPHAGGI*(1-SGGI11)*SINFG(VG11+ABG)-BETAGGI*SGGI11
dSGGI12/dT=ALPHAGGI*(1-SGGI12)*SINFG(VG12+ABG)-BETAGGI*SGGI12
dSGGI13/dT=ALPHAGGI*(1-SGGI13)*SINFG(VG13+ABG)-BETAGGI*SGGI13
dSGGI14/dT=ALPHAGGI*(1-SGGI14)*SINFG(VG14+ABG)-BETAGGI*SGGI14
dSGGI15/dT=ALPHAGGI*(1-SGGI15)*SINFG(VG15+ABG)-BETAGGI*SGGI15
dSGGI16/dT=ALPHAGGI*(1-SGGI16)*SINFG(VG16+ABG)-BETAGGI*SGGI16
dVGI1/dT=-(ITG(VGI1,RGI1)+INAG(VGI1,HGI1)+IKG(VGI1,NGI1)+IAHPG(VGI1,CAGI1)+ICAG(VGI1)+ILG(VGI1))+IAPPI-ISYNGI1-ISYNGGI1
dVGI2/dT=-(ITG(VGI2,RGI2)+INAG(VGI2,HGI2)+IKG(VGI2,NGI2)+IAHPG(VGI2,CAGI2)+ICAG(VGI2)+ILG(VGI2))+IAPPI-ISYNGI2-ISYNGGI2
dVGI3/dT=-(ITG(VGI3,RGI3)+INAG(VGI3,HGI3)+IKG(VGI3,NGI3)+IAHPG(VGI3,CAGI3)+ICAG(VGI3)+ILG(VGI3))+IAPPI-ISYNGI3-ISYNGGI3
dVGI4/dT=-(ITG(VGI4,RGI4)+INAG(VGI4,HGI4)+IKG(VGI4,NGI4)+IAHPG(VGI4,CAGI4)+ICAG(VGI4)+ILG(VGI4))+IAPPI-ISYNGI4-ISYNGGI4
dVGI5/dT=-(ITG(VGI5,RGI5)+INAG(VGI5,HGI5)+IKG(VGI5,NGI5)+IAHPG(VGI5,CAGI5)+ICAG(VGI5)+ILG(VGI5))+IAPPI-ISYNGI5-ISYNGGI5
dVGI6/dT=-(ITG(VGI6,RGI6)+INAG(VGI6,HGI6)+IKG(VGI6,NGI6)+IAHPG(VGI6,CAGI6)+ICAG(VGI6)+ILG(VGI6))+IAPPI-ISYNGI6-ISYNGGI6
dVGI7/dT=-(ITG(VGI7,RGI7)+INAG(VGI7,HGI7)+IKG(VGI7,NGI7)+IAHPG(VGI7,CAGI7)+ICAG(VGI7)+ILG(VGI7))+IAPPI-ISYNGI7-ISYNGGI7
dVGI8/dT=-(ITG(VGI8,RGI8)+INAG(VGI8,HGI8)+IKG(VGI8,NGI8)+IAHPG(VGI8,CAGI8)+ICAG(VGI8)+ILG(VGI8))+IAPPI-ISYNGI8-ISYNGGI8
dVGI9/dT=-(ITG(VGI9,RGI9)+INAG(VGI9,HGI9)+IKG(VGI9,NGI9)+IAHPG(VGI9,CAGI9)+ICAG(VGI9)+ILG(VGI9))+IAPPI-ISYNGI9-ISYNGGI9
dVGI10/dT=-(ITG(VGI10,RGI10)+INAG(VGI10,HGI10)+IKG(VGI10,NGI10)+IAHPG(VGI10,CAGI10)+ICAG(VGI10)+ILG(VGI10))+IAPPI-ISYNGI10-ISYNGGI10
dVGI11/dT=-(ITG(VGI11,RGI11)+INAG(VGI11,HGI11)+IKG(VGI11,NGI11)+IAHPG(VGI11,CAGI11)+ICAG(VGI11)+ILG(VGI11))+IAPPI-ISYNGI11-ISYNGGI11
dVGI12/dT=-(ITG(VGI12,RGI12)+INAG(VGI12,HGI12)+IKG(VGI12,NGI12)+IAHPG(VGI12,CAGI12)+ICAG(VGI12)+ILG(VGI12))+IAPPI-ISYNGI12-ISYNGGI12
dVGI13/dT=-(ITG(VGI13,RGI13)+INAG(VGI13,HGI13)+IKG(VGI13,NGI13)+IAHPG(VGI13,CAGI13)+ICAG(VGI13)+ILG(VGI13))+IAPPI-ISYNGI13-ISYNGGI13
dVGI14/dT=-(ITG(VGI14,RGI14)+INAG(VGI14,HGI14)+IKG(VGI14,NGI14)+IAHPG(VGI14,CAGI14)+ICAG(VGI14)+ILG(VGI14))+IAPPI-ISYNGI14-ISYNGGI14
dVGI15/dT=-(ITG(VGI15,RGI15)+INAG(VGI15,HGI15)+IKG(VGI15,NGI15)+IAHPG(VGI15,CAGI15)+ICAG(VGI15)+ILG(VGI15))+IAPPI-ISYNGI15-ISYNGGI15
dVGI16/dT=-(ITG(VGI16,RGI16)+INAG(VGI16,HGI16)+IKG(VGI16,NGI16)+IAHPG(VGI16,CAGI16)+ICAG(VGI16)+ILG(VGI16))+IAPPI-ISYNGI16-ISYNGGI16
dNGI1/dT=DELTANG*(NINFG(VGI1)-NGI1)/TAUNG(VGI1)
dNGI2/dT=DELTANG*(NINFG(VGI2)-NGI2)/TAUNG(VGI2)
dNGI3/dT=DELTANG*(NINFG(VGI3)-NGI3)/TAUNG(VGI3)
dNGI4/dT=DELTANG*(NINFG(VGI4)-NGI4)/TAUNG(VGI4)
dNGI5/dT=DELTANG*(NINFG(VGI5)-NGI5)/TAUNG(VGI5)
dNGI6/dT=DELTANG*(NINFG(VGI6)-NGI6)/TAUNG(VGI6)
dNGI7/dT=DELTANG*(NINFG(VGI7)-NGI7)/TAUNG(VGI7)
dNGI8/dT=DELTANG*(NINFG(VGI8)-NGI8)/TAUNG(VGI8)
dNGI9/dT=DELTANG*(NINFG(VGI9)-NGI9)/TAUNG(VGI9)
dNGI10/dT=DELTANG*(NINFG(VGI10)-NGI10)/TAUNG(VGI10)
dNGI11/dT=DELTANG*(NINFG(VGI11)-NGI11)/TAUNG(VGI11)
dNGI12/dT=DELTANG*(NINFG(VGI12)-NGI12)/TAUNG(VGI12)
dNGI13/dT=DELTANG*(NINFG(VGI13)-NGI13)/TAUNG(VGI13)
dNGI14/dT=DELTANG*(NINFG(VGI14)-NGI14)/TAUNG(VGI14)
dNGI15/dT=DELTANG*(NINFG(VGI15)-NGI15)/TAUNG(VGI15)
dNGI16/dT=DELTANG*(NINFG(VGI16)-NGI16)/TAUNG(VGI16)
dHGI1/dT=DELTAHG*(HINFG(VGI1)-HGI1)/TAUHG(VGI1)
dHGI2/dT=DELTAHG*(HINFG(VGI2)-HGI2)/TAUHG(VGI2)
dHGI3/dT=DELTAHG*(HINFG(VGI3)-HGI3)/TAUHG(VGI3)
dHGI4/dT=DELTAHG*(HINFG(VGI4)-HGI4)/TAUHG(VGI4)
dHGI5/dT=DELTAHG*(HINFG(VGI5)-HGI5)/TAUHG(VGI5)
dHGI6/dT=DELTAHG*(HINFG(VGI6)-HGI6)/TAUHG(VGI6)
dHGI7/dT=DELTAHG*(HINFG(VGI7)-HGI7)/TAUHG(VGI7)
dHGI8/dT=DELTAHG*(HINFG(VGI8)-HGI8)/TAUHG(VGI8)
dHGI9/dT=DELTAHG*(HINFG(VGI9)-HGI9)/TAUHG(VGI9)
dHGI10/dT=DELTAHG*(HINFG(VGI10)-HGI10)/TAUHG(VGI10)
dHGI11/dT=DELTAHG*(HINFG(VGI11)-HGI11)/TAUHG(VGI11)
dHGI12/dT=DELTAHG*(HINFG(VGI12)-HGI12)/TAUHG(VGI12)
dHGI13/dT=DELTAHG*(HINFG(VGI13)-HGI13)/TAUHG(VGI13)
dHGI14/dT=DELTAHG*(HINFG(VGI14)-HGI14)/TAUHG(VGI14)
dHGI15/dT=DELTAHG*(HINFG(VGI15)-HGI15)/TAUHG(VGI15)
dHGI16/dT=DELTAHG*(HINFG(VGI16)-HGI16)/TAUHG(VGI16)
dRGI1/dT=PHIG*(RINFG(VGI1)-RGI1)/TAURG
dRGI2/dT=PHIG*(RINFG(VGI2)-RGI2)/TAURG
dRGI3/dT=PHIG*(RINFG(VGI3)-RGI3)/TAURG
dRGI4/dT=PHIG*(RINFG(VGI4)-RGI4)/TAURG
dRGI5/dT=PHIG*(RINFG(VGI5)-RGI5)/TAURG
dRGI6/dT=PHIG*(RINFG(VGI6)-RGI6)/TAURG
dRGI7/dT=PHIG*(RINFG(VGI7)-RGI7)/TAURG
dRGI8/dT=PHIG*(RINFG(VGI8)-RGI8)/TAURG
dRGI9/dT=PHIG*(RINFG(VGI9)-RGI9)/TAURG
dRGI10/dT=PHIG*(RINFG(VGI10)-RGI10)/TAURG
dRGI11/dT=PHIG*(RINFG(VGI11)-RGI11)/TAURG
dRGI12/dT=PHIG*(RINFG(VGI12)-RGI12)/TAURG
dRGI13/dT=PHIG*(RINFG(VGI13)-RGI13)/TAURG
dRGI14/dT=PHIG*(RINFG(VGI14)-RGI14)/TAURG
dRGI15/dT=PHIG*(RINFG(VGI15)-RGI15)/TAURG
dRGI16/dT=PHIG*(RINFG(VGI16)-RGI16)/TAURG
dCAGI1/dT=EPSG*(-ICAG(VGI1)-ITG(VGI1,RGI1)-KCAGI*CAGI1)
dCAGI2/dT=EPSG*(-ICAG(VGI2)-ITG(VGI2,RGI2)-KCAGI*CAGI2)
dCAGI3/dT=EPSG*(-ICAG(VGI3)-ITG(VGI3,RGI3)-KCAGI*CAGI3)
dCAGI4/dT=EPSG*(-ICAG(VGI4)-ITG(VGI4,RGI4)-KCAGI*CAGI4)
dCAGI5/dT=EPSG*(-ICAG(VGI5)-ITG(VGI5,RGI5)-KCAGI*CAGI5)
dCAGI6/dT=EPSG*(-ICAG(VGI6)-ITG(VGI6,RGI6)-KCAGI*CAGI6)
dCAGI7/dT=EPSG*(-ICAG(VGI7)-ITG(VGI7,RGI7)-KCAGI*CAGI7)
dCAGI8/dT=EPSG*(-ICAG(VGI8)-ITG(VGI8,RGI8)-KCAGI*CAGI8)
dCAGI9/dT=EPSG*(-ICAG(VGI9)-ITG(VGI9,RGI9)-KCAGI*CAGI9)
dCAGI10/dT=EPSG*(-ICAG(VGI10)-ITG(VGI10,RGI10)-KCAGI*CAGI10)
dCAGI11/dT=EPSG*(-ICAG(VGI11)-ITG(VGI11,RGI11)-KCAGI*CAGI11)
dCAGI12/dT=EPSG*(-ICAG(VGI12)-ITG(VGI12,RGI12)-KCAGI*CAGI12)
dCAGI13/dT=EPSG*(-ICAG(VGI13)-ITG(VGI13,RGI13)-KCAGI*CAGI13)
dCAGI14/dT=EPSG*(-ICAG(VGI14)-ITG(VGI14,RGI14)-KCAGI*CAGI14)
dCAGI15/dT=EPSG*(-ICAG(VGI15)-ITG(VGI15,RGI15)-KCAGI*CAGI15)
dCAGI16/dT=EPSG*(-ICAG(VGI16)-ITG(VGI16,RGI16)-KCAGI*CAGI16)
dSGI1/dT=ALPHAG*(1-SGI1)*SINFG(VGI1+ABG)-BETAGI*SGI1
dSGI2/dT=ALPHAG*(1-SGI2)*SINFG(VGI2+ABG)-BETAGI*SGI2
dSGI3/dT=ALPHAG*(1-SGI3)*SINFG(VGI3+ABG)-BETAGI*SGI3
dSGI4/dT=ALPHAG*(1-SGI4)*SINFG(VGI4+ABG)-BETAGI*SGI4
dSGI5/dT=ALPHAG*(1-SGI5)*SINFG(VGI5+ABG)-BETAGI*SGI5
dSGI6/dT=ALPHAG*(1-SGI6)*SINFG(VGI6+ABG)-BETAGI*SGI6
dSGI7/dT=ALPHAG*(1-SGI7)*SINFG(VGI7+ABG)-BETAGI*SGI7
dSGI8/dT=ALPHAG*(1-SGI8)*SINFG(VGI8+ABG)-BETAGI*SGI8
dSGI9/dT=ALPHAG*(1-SGI9)*SINFG(VGI9+ABG)-BETAGI*SGI9
dSGI10/dT=ALPHAG*(1-SGI10)*SINFG(VGI10+ABG)-BETAGI*SGI10
dSGI11/dT=ALPHAG*(1-SGI11)*SINFG(VGI11+ABG)-BETAGI*SGI11
dSGI12/dT=ALPHAG*(1-SGI12)*SINFG(VGI12+ABG)-BETAGI*SGI12
dSGI13/dT=ALPHAG*(1-SGI13)*SINFG(VGI13+ABG)-BETAGI*SGI13
dSGI14/dT=ALPHAG*(1-SGI14)*SINFG(VGI14+ABG)-BETAGI*SGI14
dSGI15/dT=ALPHAG*(1-SGI15)*SINFG(VGI15+ABG)-BETAGI*SGI15
dSGI16/dT=ALPHAG*(1-SGI16)*SINFG(VGI16+ABG)-BETAGI*SGI16
dVT1/dT=-(ILF(VT1)+INAF(VT1,MINFTC(VT1),HT1)+IKF(VT1,.75*(1-HT1))+ITF(VT1,MTINF(VT1),HTT1))+FTC(T)-GSYNTC*SGA*(VT1-VSYNTC)
dHTT1/dT=QHT*(HTINF(VT1)-HTT1)/TAUHT(VT1)
dHT1/dT=TADJ*(HINFTC(VT1)-HT1)/TAUHTC(VT1)
dVT2/dT=-(ILF(VT2)+INAF(VT2,MINFTC(VT2),HT2)+IKF(VT2,.75*(1-HT2))+ITF(VT2,MTINF(VT2),HTT2))+FTC(T)-GSYNTC*SGB*(VT2-VSYNTC)
dHTT2/dT=QHT*(HTINF(VT2)-HTT2)/TAUHT(VT2)
dHT2/dT=TADJ*(HINFTC(VT2)-HT2)/TAUHTC(VT2)
dCOUN1/dT=0
dCOUN2/dT=0
dZ1/dT=1-HV(Z1-50)
dZ2/dT=1-HV(Z2-50)
dGOOD1/dT=0
dGOOD2/dT=0
SG=FTC(T)
SGT=FTC2(T)
SA=SGA
SB=SGB

where ...
ISYN1 = GSYN*(SG5+SG2)*(V1-VSYN)   
ISYN2 = GSYN*(SG6+SG1)*(V2-VSYN)   
ISYN3 = GSYN*(SG8+SG4)*(V3-VSYN)   
ISYN4 = GSYN*(SG7+SG3)*(V4-VSYN)   
ISYN5 = GSYN*(SG2+SG6)*(V5-VSYN)   
ISYN6 = GSYN*(SG1+SG5)*(V6-VSYN)   
ISYN7 = GSYN*(SG3+SG8)*(V7-VSYN)   
ISYN8 = GSYN*(SG4+SG7)*(V8-VSYN)   
ISYN9 = GSYN*(SG13+SG10)*(V9-VSYN)   
ISYN10 = GSYN*(SG14+SG9)*(V10-VSYN)   
ISYN11 = GSYN*(SG16+SG12)*(V11-VSYN)   
ISYN12 = GSYN*(SG15+SG11)*(V12-VSYN)   
ISYN13 = GSYN*(SG10+SG14)*(V13-VSYN)   
ISYN14 = GSYN*(SG9+SG13)*(V14-VSYN)   
ISYN15 = GSYN*(SG11+SG16)*(V15-VSYN)   
ISYN16 = GSYN*(SG12+SG15)*(V16-VSYN)   
ISYNG1 = GSYNG*(S8+S4+.2*S16)*(VG1-VSYNG)   
ISYNG2 = GSYNG*(S7+S3+.2*S15)*(VG2-VSYNG)   
ISYNG3 = GSYNG*(S1+S5+.2*S9)*(VG3-VSYNG)   
ISYNG4 = GSYNG*(S2+S6+.2*S10)*(VG4-VSYNG)   
ISYNG5 = GSYNG*(S4+S8+.2*S12)*(VG5-VSYNG)   
ISYNG6 = GSYNG*(S3+S7+.2*S11)*(VG6-VSYNG)   
ISYNG7 = GSYNG*(S5+S2+.2*S13)*(VG7-VSYNG)   
ISYNG8 = GSYNG*(S6+S1+.2*S14)*(VG8-VSYNG)   
ISYNG9 = GSYNG*(S16+S12+.2*S8)*(VG9-VSYNG)   
ISYNG10 = GSYNG*(S15+S11+.2*S7)*(VG10-VSYNG)   
ISYNG11 = GSYNG*(S9+S13+.2*S1)*(VG11-VSYNG)   
ISYNG12 = GSYNG*(S10+S14+.2*S2)*(VG12-VSYNG)   
ISYNG13 = GSYNG*(S12+S16+.2*S4)*(VG13-VSYNG)   
ISYNG14 = GSYNG*(S11+S15+.2*S3)*(VG14-VSYNG)   
ISYNG15 = GSYNG*(S13+S10+.2*S5)*(VG15-VSYNG)   
ISYNG16 = GSYNG*(S14+S9+.2*S6)*(VG16-VSYNG)   
ISYNGG1 = GSYNGG*(SG2+SG3)*(VG1-VSYNGG)   
ISYNGG2 = GSYNGG*(SG1+SG5)*(VG2-VSYNGG)   
ISYNGG3 = GSYNGG*(SG4+SG8)*(VG3-VSYNGG)   
ISYNGG4 = GSYNGG*(SG3+SG1)*(VG4-VSYNGG)   
ISYNGG5 = GSYNGG*(SG6+SG7)*(VG5-VSYNGG)   
ISYNGG6 = GSYNGG*(SG5+SG2)*(VG6-VSYNGG)   
ISYNGG7 = GSYNGG*(SG8+SG3)*(VG7-VSYNGG)   
ISYNGG8 = GSYNGG*(SG7+SG4)*(VG8-VSYNGG)   
ISYNGG9 = GSYNGG*(SG10+SG11)*(VG9-VSYNGG)   
ISYNGG10 = GSYNGG*(SG9+SG13)*(VG10-VSYNGG)   
ISYNGG11 = GSYNGG*(SG12+SG16)*(VG11-VSYNGG)   
ISYNGG12 = GSYNGG*(SG11+SG9)*(VG12-VSYNGG)   
ISYNGG13 = GSYNGG*(SG14+SG15)*(VG13-VSYNGG)   
ISYNGG14 = GSYNGG*(SG13+SG10)*(VG14-VSYNGG)   
ISYNGG15 = GSYNGG*(SG16+SG11)*(VG15-VSYNGG)   
ISYNGG16 = GSYNGG*(SG15+SG12)*(VG16-VSYNGG)   
ISYNGI1 = GSYNGI*SI1*(VGI1-VSYNG)   
ISYNGI2 = GSYNGI*SI2*(VGI2-VSYNG)   
ISYNGI3 = GSYNGI*SI3*(VGI3-VSYNG)   
ISYNGI4 = GSYNGI*SI4*(VGI4-VSYNG)   
ISYNGI5 = GSYNGI*SI5*(VGI5-VSYNG)   
ISYNGI6 = GSYNGI*SI6*(VGI6-VSYNG)   
ISYNGI7 = GSYNGI*SI7*(VGI7-VSYNG)   
ISYNGI8 = GSYNGI*SI8*(VGI8-VSYNG)   
ISYNGI9 = GSYNGI*SI9*(VGI9-VSYNG)   
ISYNGI10 = GSYNGI*SI10*(VGI10-VSYNG)   
ISYNGI11 = GSYNGI*SI11*(VGI11-VSYNG)   
ISYNGI12 = GSYNGI*SI12*(VGI12-VSYNG)   
ISYNGI13 = GSYNGI*SI13*(VGI13-VSYNG)   
ISYNGI14 = GSYNGI*SI14*(VGI14-VSYNG)   
ISYNGI15 = GSYNGI*SI15*(VGI15-VSYNG)   
ISYNGI16 = GSYNGI*SI16*(VGI16-VSYNG)   
ISYNGGI1 = GSYNGGI*SGGI1*(VGI1-VSYN)   
ISYNGGI2 = GSYNGGI*SGGI2*(VGI2-VSYN)   
ISYNGGI3 = GSYNGGI*SGGI3*(VGI3-VSYN)   
ISYNGGI4 = GSYNGGI*SGGI4*(VGI4-VSYN)   
ISYNGGI5 = GSYNGGI*SGGI5*(VGI5-VSYN)   
ISYNGGI6 = GSYNGGI*SGGI6*(VGI6-VSYN)   
ISYNGGI7 = GSYNGGI*SGGI7*(VGI7-VSYN)   
ISYNGGI8 = GSYNGGI*SGGI8*(VGI8-VSYN)   
ISYNGGI9 = GSYNGGI*SGGI9*(VGI9-VSYN)   
ISYNGGI10 = GSYNGGI*SGGI10*(VGI10-VSYN)   
ISYNGGI11 = GSYNGGI*SGGI11*(VGI11-VSYN)   
ISYNGGI12 = GSYNGGI*SGGI12*(VGI12-VSYN)   
ISYNGGI13 = GSYNGGI*SGGI13*(VGI13-VSYN)   
ISYNGGI14 = GSYNGGI*SGGI14*(VGI14-VSYN)   
ISYNGGI15 = GSYNGGI*SGGI15*(VGI15-VSYN)   
ISYNGGI16 = GSYNGGI*SGGI16*(VGI16-VSYN)   
SGA = SGI1+SGI2+SGI5+SGI6+SGI9+SGI10+SGI13+SGI14   
SGB = SGI3+SGI4+SGI7+SGI8+SGI11+SGI12+SGI15+SGI16   

User-defined functions:
SINF(V) = 1./(1.+EXP(-(V+THETAS)/SS))  
MINF(V) = 1./(1.+EXP(-(V+THETAM)/SM))  
HINF(V) = 1./(1.+EXP((V-THETAH)/SH))  
NINF(V) = 1./(1.+EXP((V-THETAN)/SN))  
TAUN(V) = TAUN0+TAUN1/(1+EXP((V+THN)/SIGMAN))  
TAUH(V) = TAUH0+TAUH1/(1+EXP((V+THH)/SIGMAH))  
RINF(V) = 1/(1+EXP((V-THETAR)/KR))  
TAUR(V) = TAUR0+TAUR1/(1+EXP((V+THR)/SIGMAR))  
TINF(V) = 1/(1+EXP((V-THETAT)/KT))  
RNEW(R) = 1/(1+EXP((R-RTH)/RSIG))-1/(1+EXP(-RTH/RSIG))  
SINFG(VG) = 1/(1+EXP(-(VG-THETASG)/KSG))  
SINF1G(VG) = 1/(1+EXP(-(VG-THETAS1G)/KS1G))  
RINFG(VG) = 1/(1+EXP(-(VG-THETARG)/KRG))  
MINFG(VG) = 1./(1.+EXP(-(VG-THETAMG)/SIGMAMG))  
NINFG(VG) = 1./(1.+EXP(-(VG-THETANG)/SIGMANG))  
TAUNG(VG) = TAUN0G+TAUN1G/(1+EXP(-(VG-THNG)/SNG))  
HINFG(VG) = 1./(1.+EXP(-(VG-THETAHG)/SIGMAHG))  
TAUHG(VG) = TAUH0G+TAUH1G/(1+EXP(-(VG-THHG)/SHG))  
HV(V) = 1./(1.+EXP(-V/.001))  
INAF(V,M,H) = GNABAR*M^3*H*(V-ENA)  
IKF(V,N) = GKBAR*N^4*(V-EK)  
ILF(V) = GLBAR*(V-ELEAK)  
ITF(V,MT,HT) = GTBAR*MT^2*HT*V  
MINFTC(V) =   1/(1+EXP(-(V+37)/7))  
MTINF(V) =  1/(1+EXP(-(V+60)/6.2))  
AH(V) =   0.128*EXP(-(46+V)/18)  
BH(V) =   APR/(1+EXP(-(23+V)/5))  
TAUHTC(V) =   1/(AH(V)+BH(V))  
HINFTC(V) =   1/(1+EXP((V+41)/4))  
HTINF(V) =  1/(1+EXP((V+84)/4))  
TAUHT(V) = (28+APT*EXP((V+25)/(-10.5)))  
FTC(T) = ITC*HV(SIN(6.2831853*(T+SHI)/PERIOD))*(1-HV(SIN(6.2831853*(T+SHI+DUR)/PERIOD)))  
FTC2(T) = HV(SIN(6.2831853*(T+SHI2)/PERIOD))*(1-HV(SIN(6.2831853*(T+SHI2+DUR2)/PERIOD)))  
IL(V) = GL*(V-VL)  
INA(V,H) = GNA*(MINF(V))^3*H*(V-VNA)  
IK(V,N) = GK*N^4*(V-VK)  
IAHP(V,CA) = GAHP*(V-VK)*CA/(CA+K1)  
ICA(V) = GCA*((SINF(V))^2)*(V-VCA)  
IT(V,R) = GT*(TINF(V)**3)*(RNEW(R)^2)*(V-VCA)  
ITG(VG,RG) = GTG*(SINFG(VG)^3)*RG*(VG-VCAG)  
INAG(VG,HG) = GNAG*(MINFG(VG)^3)*HG*(VG-VNAG)  
IKG(VG,NG) = GKG*(NG^4)*(VG-VKG)  
IAHPG(VG,CAG) = GAHPG*(VG-VKG)*CAG/(CAG+K1G)  
ICAG(VG) = GCAG*((SINF1G(VG))^2)*(VG-VCAG)  
ILG(VG) = GLG*(VG-VLG)  

Loading data, please wait...