Circadian clock model in mammals (detailed version) (Kim & Forger 2012)

 Download zip file 
Help downloading and running models
Accession:145801
"… To understand the biochemical mechanisms of this timekeeping, we have developed a detailed mathematical model of the mammalian circadian clock. Our model can accurately predict diverse experimental data including the phenotypes of mutations or knockdown of clock genes as well as the time courses and relative expression of clock transcripts and proteins. Using this model, we show how a universal motif of circadian timekeeping, where repressors tightly bind activators rather than directly binding to DNA, can generate oscillations when activators and repressors are in stoichiometric balance. …"
Reference:
1 . Kim JK, Forger DB (2012) A mechanism for robust circadian timekeeping via stoichiometric balance Molecular Systems Biology 8:630:1-14 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Molecular Network;
Brain Region(s)/Organism:
Cell Type(s): Suprachiasmatic nucleus (SCN) neuron;
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: XPP; MATLAB; Mathematica;
Model Concept(s): Oscillations; Brain Rhythms; Circadian Rhythms;
Implementer(s): Kim, Jae Kyoung [kimjack0 at kaist.ac.kr];
********** MODEL NAME
Detailed Modelh

********** MODEL NOTES


********** MODEL STATES
d/dt(GR) = -(unbin*GR) + bin*(1 - G - GR)*(x01011 + x02011)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
d/dt(G) = -(unbin*G) + bin*(1 - G - GR)*x00011                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(GrR) = -(unbinr*GrR) + binr*(1 - Gr - GrR)*(x01011 + x02011)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(Gr) = -(unbinr*Gr) + binr*(1 - Gr - GrR)*x00011                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
d/dt(GcR) = -(unbinc*GcR) + binc*(1 - Gc - GcR)*(x01011 + x02011)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(Gc) = -(unbinc*Gc) + binc*(1 - Gc - GcR)*x00011                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
d/dt(GBR) = -(unbinrev*GBR) + binrev*GB*(revn + revng + revngp + revnp)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
d/dt(GB) = unbinrev*GBR - binrev*GB*(revn + revng + revngp + revnp)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
d/dt(GBRb) = -(unbinrevb*GBRb) + binrevb*GBb*(revn + revng + revngp + revnp)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
d/dt(GBb) = unbinrevb*GBRb - binrevb*GBb*(revn + revng + revngp + revnp)d/dt(MnPo) = trPo*G - tmc*MnPo - umPo*MnPo                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(McPo) = -(umPo*McPo) + tmc*MnPod/dt(MnPt) = trPt*G - tmc*MnPt - umPt*MnPtd/dt(McPt) = -(umPt*McPt) + tmc*MnPt                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(MnRt) = trRt*Gc - tmc*MnRt - umRt*MnRtd/dt(McRt) = -(umRt*McRt) + tmc*MnRt                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
d/dt(MnRev) = -(tmcrev*MnRev) - umRev*MnRev + trRev*Gr*x00011d/dt(McRev) = -(umRev*McRev) + tmcrev*MnRev                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           
d/dt(MnRo) = trRo*G*GB - tmc*MnRo - umRo*MnRod/dt(McRo) = -(umRo*McRo) + tmc*MnRo                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(MnB) = trB*GBb - tmc*MnB - umB*MnBd/dt(McB) = -(umB*McB) + tmc*MnBd/dt(MnNp) = trNp*GB - tmc*MnNp - umNp*MnNp                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(McNp) = -(umNp*McNp) + tmc*MnNpd/dt(B) = -(ub*B) + uncbin*BC - cbin*B*Cl + tlb*McB                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
d/dt(Cl) = tlc + uncbin*BC - uc*Cl - cbin*B*Cl + tlnp*McNpd/dt(BC) = -(phos*BC) - ubc*BC - uncbin*BC + cbin*B*Cl                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
d/dt(cyrev) = -((nlrev + urev)*cyrev) + dg*cyrevg + tlrev*McRev + nerev*revn - ag*cyrev*x00200                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(revn) = nlrev*cyrev + (-nerev - urev)*revn + dg*revng - ag*Nf*revn*x00210                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(cyrevg) = -(cyrevg*(dg + nlrev + urev + gto)) + nerev*revng + ag*cyrev*x00200                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(revng) = nlrev*cyrevg - (dg + nerev + urev + gto)*revng + ag*Nf*revn*x00210                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
d/dt(cyrevgp) = -((dg + nlrev + uprev)*cyrevgp) + cyrevg*gto + nerev*revngp                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
d/dt(revngp) = nlrev*cyrevgp + gto*revng - (dg + nerev + uprev)*revngp                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
d/dt(cyrevp) = dg*cyrevgp - (nlrev + uprev)*cyrevp + nerev*revnpd/dt(revnp) = nlrev*cyrevp + dg*revngp - (nerev + uprev)*revnp                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(gto) = trgto*G*GB - ugto*gto                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x00001) = phos*BC - nlbc*x00001 - ubc*x00001                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x00011) = nlbc*x00001 - ubc*x00011 + uro*x01011 - cbbin*Nf*x00011*(x01010 + x02010) + urt*x02011 +  uncbbin*(x01011 + x02011) + upu*(x50011 + x50111 + x50211 + x50311) +up*(x20011 + x20111 +  x40011 + x40111 + x40211 + x40311 + x60011 + x60111 + x60211 + x60311) - bbin*Nf*x00011*(  x20010 + x20110  + x21010 + x21110 +  x22010 + x22110  +  x40010 + x40110 + x40210 + x40310 + x41010 + x41110 + x41210 + x41310 + x42010 + x42110 + x42210 + x42310 + x50010 + x50110 + x50210 + x50310 + x51010 + x51110 + x51210 + x51310 + x52010 + x52110 + x52210 + x52310 + x60010 + x60110 + x60210 + x60310 + x61010 + x61110 + x61210 + x61310 + x62010 + x62110 + x62210 + x62310) + unbbin*(  x20011 + x20111 + x21011 + x21111 +  x22011 + x22111  + x40011 + x40111 + x40211 + x40311 + x41011 + x41111 + x41211 + x41311 + x42011 + x42111 + x42211 + x42311 + x50011 + x50111 + x50211 + x50311 + x51011 + x51111 + x51211 + x51311 + x52011 + x52111 + x52211 + x52311 + x60011 + x60111 + x60211 + x60311 + x61011 + x61111 + x61211 + x61311 + x62011 + x62111 + x62211 + x62311)
d/dt(x00100) = lne*x00110 + upu*(x10100 + x30100 +  x30300 + x50100 +  x50300) +up*(x20100 + x40100 + x40300 + x60100 +  x60300) - ac*x00100*(x10000 +  x20000 + x21000 + x22000 + x30000 +  x40000 + x41000 + x42000 + x50000 + x51000 + x52000 + x60000 + x61000 + x62000) +dc*(x10100 + x20100 + x21100 + x22100 + x30100 +  x40100 + x41100 + x42100 + x50100 + x51100 + x52100 + x60100 + x61100 + x62100) - ac*x00100*(x30200 +  x40200 + x41200 + x42200 + x50200 + x51200 + x52200 + x60200 + x61200 + x62200) +dc*(x30300 + x40300 + x41300 + x42300 + x50300 + x51300 + x52300 + x60300 + x61300 + x62300)                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
d/dt(x00110) = -(lne*x00110) + upu*(  x50110 + x50111 + x50310 + x50311) +up*(x20110 + x20111 + x40110 + x40111 + x40310 + x40311 + x60110 + x60111 + x60310 + x60311) - ac*Nf*x00110*(x20010 + x21010 + x22010 +  x40010 + x41010 + x42010 + x50010 + x51010 + x52010 + x60010 + x61010 + x62010) - ac*Nf*x00110*( x20011 + x21011 + x22011 +  x40011 + x41011 + x42011 + x50011 + x51011 + x52011 + x60011 + x61011 + x62011) +dc*(  x20110 + x21110 + x22110 + x40110 + x41110 + x42110 + x50110 + x51110 + x52110 + x60110 + x61110 + x62110) +dc*( x20111 + x21111 + x22111 + x40111 + x41111 + x42111 + x50111 + x51111 + x52111 + x60111 + x61111 + x62111) - ac*Nf*x00110*(  x40210 + x41210 + x42210 + x50210 + x51210 + x52210 + x60210 + x61210 + x62210) - ac*Nf*x00110*( x40211 + x41211 + x42211 + x50211 + x51211 + x52211 + x60211 + x61211 + x62211) +dc*(  x40310 + x41310 + x42310 + x50310 + x51310 + x52310 + x60310 + x61310 + x62310) +dc*(x40311 + x41311 + x42311 + x50311 + x51311 + x52311 + x60311 + x61311 + x62311)                                                  
d/dt(x00200) =dg*cyrevg + urev*cyrevg + dg*cyrevgp + uprev*cyrevgp - ag*cyrev*x00200 + lne*x00210 +upu*(  x30200 + x30300 +  x50200 + x50300 ) +up*(x40200 + x40300 +  x60200 + x60300 ) - agp*x00200*(x30000 +  x30100  +  x40000 + x40100 + x41000 +  x41100 + x42000 +  x42100 +  x50000 + x50100 +  x51000 + x51100 + x52000  + x52100 + x60000 + x60100 +  x61000 + x61100  + x62000 +  x62100 ) + dg*(x30200 +  x30300 +  x40200 +  x40300 +  x41200 + x41300 + x42200 + x42300 +x50200 + x50300 + x51200  + x51300 + x52200 +  x52300 + x60200 + x60300 +  x61200 + x61300 +  x62200 +x62300 )                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
d/dt(x00210) =dg*revng + urev*revng + dg*revngp + uprev*revngp - lne*x00210 - ag*Nf*revn*x00210 +upu*(x50210 + x50211 + x50310 + x50311) +up*( x40210 + x40211 + x40310 + x40311 + x60210 + x60211 + x60310 + x60311) - agp*Nf*x00210*(   x40010 + x40011 + x40110 + x40111 + x41010 + x41011 + x41110 + x41111 + x42010 + x42011 + x42110 + x42111 + x50010 + x50011 + x50110 + x50111 + x51010 + x51011 + x51110 + x51111 + x52010 + x52011 + x52110 + x52111 + x60010 + x60011 + x60110 + x60111 + x61010 + x61011 + x61110 + x61111 + x62010 + x62011 + x62110 + x62111) + dg*( x40210 + x40211 + x40310 + x40311 + x41210 + x41211 + x41310 + x41311 + x42210 + x42211 + x42310 + x42311 + x50210 + x50211 + x50310 + x50311 + x51210 + x51211 + x51310 + x51311 + x52210 + x52211 + x52310 + x52311 + x60210 + x60211 + x60310 + x60311 + x61210 + x61211 + x61310 + x61311 + x62210 + x62211 + x62310 + x62311)                                                                                                                                                                             
d/dt(x01000) =tlr*McRo - uro*x01000 - ar*x01000*(x20000 + x20100 +  x40000 + x40100 + x40200 + x40300 + x50000 + x50100 + x50200 + x50300 + x60000 + x60100 + x60200 + x60300) +dr*(x21000 + x21100 +  x41000 + x41100 + x41200 + x41300 + x51000 + x51100 + x51200 + x51300 + x61000 + x61100 + x61200 + x61300)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x01010) =-(uro*x01010) - cbbin*Nf*x00011*x01010 + uncbbin*x01011 - ar*Nf*x01010*(x20010 + x20110 + x40010 + x40110 + x40210 + x40310 + x50010 + x50110 + x50210 + x50310 + x60010 + x60110 + x60210 + x60310) -ar*Nf*x01010* (x20011 + x20111 + x40011 + x40111 + x40211 + x40311 + x50011 + x50111 + x50211 + x50311 + x60011 + x60111 + x60211 + x60311)+ dr*(x21010 + x21110 +  x41010 + x41110 + x41210 + x41310 + x51010 + x51110 + x51210 + x51310 + x61010 + x61110 + x61210 + x61310) + dr*(x21011 + x21111 + x41011 + x41111 + x41211 + x41311 + x51011 + x51111 + x51211 + x51311 + x61011 + x61111 + x61211 + x61311)                                                                                                                                                                                                                                                                                                                                                                                                                                                              
d/dt(x01011) =cbbin*Nf*x00011*x01010 - uncbbin*x01011 - uro*x01011 - ar*Nf*x01011*(x20010 + x20110  + x40010 + x40110 + x40210 + x40310 + x50010 + x50110 + x50210 + x50310 + x60010 + x60110 + x60210 + x60310) + dr*(x21011 + x21111 +  x41011 + x41111 + x41211 + x41311 + x51011 + x51111 + x51211 + x51311 + x61011 + x61111 + x61211 + x61311)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
d/dt(x02000) =tlr*McRt - urt*x02000 - ar*x02000*(x20000 + x20100 +  x40000 + x40100 + x40200 + x40300 + x50000 + x50100 + x50200 + x50300 + x60000 + x60100 + x60200 + x60300) +dr*(x22000 + x22100 + x42000 + x42100 + x42200 + x42300 + x52000 + x52100 + x52200 + x52300 + x62000 + x62100 + x62200 + x62300)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
d/dt(x02010) =-(urt*x02010) - cbbin*Nf*x00011*x02010 + uncbbin*x02011 - ar*Nf*x02010*(x20010 + x20110 + x40010 + x40110 + x40210 + x40310 + x50010 + x50110 + x50210 + x50310 + x60010 + x60110 + x60210 + x60310) - ar*Nf*x02010*(x20011 + x20111 + x40011 + x40111 + x40211 + x40311 + x50011 + x50111 + x50211 + x50311 + x60011 + x60111 + x60211 + x60311)+ dr*(x22010 + x22110 + x42010 + x42110 + x42210 + x42310 + x52010 + x52110 + x52210 + x52310 + x62010 + x62110 + x62210 + x62310) + dr*(x22011 + x22111 +  x42011 + x42111 + x42211 + x42311 + x52011 + x52111 + x52211 + x52311 + x62011 + x62111 + x62211 + x62311)                                                                                                                                                                                                                                                                                                                                                                                                                                                              
d/dt(x02011) =cbbin*Nf*x00011*x02010 - uncbbin*x02011 - urt*x02011 - ar*Nf*x02011*(x20010 + x20110  + x40010 + x40110 + x40210 + x40310 + x50010 + x50110 + x50210 + x50310 + x60010 + x60110 + x60210 + x60310) + dr*(x22011 + x22111  + x42011 + x42111 + x42211 + x42311 + x52011 + x52111 + x52211 + x52311 + x62011 + x62111 + x62211 + x62311)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
d/dt(x10000) = tlp*McPo - upu*x10000 - ac*x00100*x10000 + dc*x10100                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
d/dt(x10100) =ac*x00100*x10000 - dc*x10100 - hoo*x10100 - upu*x10100                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
d/dt(x20000) =-(nl*x20000) - up*x20000 - ac*x00100*x20000 - ar*(x01000 + x02000)*x20000 +ne*x20010 + dc*x20100 + dr*(x21000 + x22000)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
d/dt(x20010) =nl*x20000 - ne*x20010 - up*x20010 - bbin*Nf*x00011*x20010 - ac*Nf*x00110*x20010 - ar*Nf*(x01010 + x02010)*x20010 - ar*Nf*(x01011 + x02011)*x20010 + ubc*x20011 +unbbin*x20011 + dc*x20110 + dr*(x21010 + x22010) + dr*(x21011 + x22011)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
d/dt(x20011) =bbin*Nf*x00011*x20010 - ubc*x20011 - unbbin*x20011 - up*x20011 - ac*Nf*x00110*x20011 - ar*Nf*(x01010 + x02010)*x20011 + dc*x20111 + dr*(x21011 + x22011)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
d/dt(x20100) =hoo*x10100 + ac*x00100*x20000 - dc*x20100 - nl*x20100 - up*x20100 - ar*(x01000 + x02000)*x20100 + ne*x20110 + dr*(x21100 + x22100)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
d/dt(x20110) =ac*Nf*x00110*x20010 + nl*x20100 - dc*x20110 - ne*x20110 - up*x20110 - bbin*Nf*x00011*x20110 - ar*Nf*(x01010 + x02010)*x20110 - ar*Nf*(x01011 + x02011)*x20110 + ubc*x20111 + unbbin*x20111 + dr*(x21110 + x22110) +dr*(x21111 + x22111)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
d/dt(x20111) =ac*Nf*x00110*x20011 + bbin*Nf*x00011*x20110 - dc*x20111 - ubc*x20111 - unbbin*x20111 - up*x20111 - ar*Nf*(x01010 + x02010)*x20111 + dr*(x21111 + x22111)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
d/dt(x21000) =ar*x01000*x20000 - dr*x21000 - nl*x21000 - ac*x00100*x21000 + ne*x21010 +dc*x21100                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
d/dt(x21010) =ar*Nf*x01010*x20010 + nl*x21000 - dr*x21010 - ne*x21010 - bbin*Nf*x00011*x21010 - ac*Nf*x00110*x21010 + unbbin*x21011 + dc*x21110                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
d/dt(x21011) =ar*Nf*x01011*x20010 + ar*Nf*x01010*x20011 + bbin*Nf*x00011*x21010 - 2*dr*x21011 - unbbin*x21011 - ac*Nf*x00110*x21011 + dc*x21111                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
d/dt(x21100) =ar*x01000*x20100 + ac*x00100*x21000 - dc*x21100 - dr*x21100 - nl*x21100 +ne*x21110                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
d/dt(x21110) =ar*Nf*x01010*x20110 + ac*Nf*x00110*x21010 + nl*x21100 - dc*x21110 - dr*x21110 - ne*x21110 - bbin*Nf*x00011*x21110 + unbbin*x21111                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
d/dt(x21111) =ar*Nf*x01011*x20110 + ar*Nf*x01010*x20111 + ac*Nf*x00110*x21011 +bbin*Nf*x00011*x21110 - dc*x21111 - 2*dr*x21111 - unbbin*x21111                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x22000) =ar*x02000*x20000 - dr*x22000 - nl*x22000 - ac*x00100*x22000 + ne*x22010 +dc*x22100                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
d/dt(x22010) =ar*Nf*x02010*x20010 + nl*x22000 - dr*x22010 - ne*x22010 - bbin*Nf*x00011*x22010 - ac*Nf*x00110*x22010 + unbbin*x22011 + dc*x22110                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
d/dt(x22011) =ar*Nf*x02011*x20010 + ar*Nf*x02010*x20011 + bbin*Nf*x00011*x22010 - 2*dr*x22011 - unbbin*x22011 - ac*Nf*x00110*x22011 + dc*x22111                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
d/dt(x22100) =ar*x02000*x20100 + ac*x00100*x22000 - dc*x22100 - dr*x22100 - nl*x22100 +ne*x22110                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
d/dt(x22110) =ar*Nf*x02010*x20110 + ac*Nf*x00110*x22010 + nl*x22100 - dc*x22110 - dr*x22110 - ne*x22110 - bbin*Nf*x00011*x22110 + unbbin*x22111                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
d/dt(x22111) =ar*Nf*x02011*x20110 + ar*Nf*x02010*x20111 + ac*Nf*x00110*x22011 +bbin*Nf*x00011*x22110 - dc*x22111 - 2*dr*x22111 - unbbin*x22111                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x30000) =tlp*McPt - upu*x30000 - ac*x00100*x30000 - agp*x00200*x30000 + dc*x30100 + dg*x30200                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x30100) =ac*x00100*x30000 - dc*x30100 - hto*x30100 - upu*x30100 - agp*x00200*x30100 +dg*x30300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
d/dt(x30200) =agp*x00200*x30000 - dg*x30200 - upu*x30200 - gto*x30200 - ac*x00100*x30200 +dc*x30300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
d/dt(x30300) =agp*x00200*x30100 + ac*x00100*x30200 - dc*x30300 - dg*x30300 - hto*x30300 - upu*x30300 - gto*x30300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x40000) =-(nl*x40000) - up*x40000 - ac*x00100*x40000 - agp*x00200*x40000 - ar*(x01000 + x02000)*x40000 + ne*x40010 + dc*x40100 + dg*x40200 +dr*(x41000 + x42000)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
d/dt(x40010) =nl*x40000 - ne*x40010 - up*x40010 - bbin*Nf*x00011*x40010 - ac*Nf*x00110*x40010 - agp*Nf*x00210*x40010 - ar*Nf*(x01010 + x02010)*x40010 - ar*Nf*(x01011 + x02011)*x40010 + ubc*x40011 + unbbin*x40011 + dc*x40110 + dg*x40210 +dr*(x41010 + x42010) + dr*(x41011 + x42011)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           
d/dt(x40011) =bbin*Nf*x00011*x40010 - ubc*x40011 - unbbin*x40011 - up*x40011 - ac*Nf*x00110*x40011 - agp*Nf*x00210*x40011 - ar*Nf*(x01010 + x02010)*x40011 + dc*x40111 + dg*x40211 +dr*(x41011 + x42011)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           
d/dt(x40100) =hto*x30100 + ac*x00100*x40000 - dc*x40100 - nl*x40100 - up*x40100 - agp*x00200*x40100 - ar*(x01000 + x02000)*x40100 + ne*x40110 + dg*x40300 +dr*(x41100 + x42100)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
d/dt(x40110) = ac*Nf*x00110*x40010 + nl*x40100 - dc*x40110 - ne*x40110 - up*x40110 - bbin*Nf*x00011*x40110 - agp*Nf*x00210*x40110 - ar*Nf*(x01010 + x02010)*x40110 - ar*Nf*(x01011 + x02011)*x40110 + ubc*x40111 + unbbin*x40111 + dg*x40310 +dr*(x41110 + x42110) + dr*(x41111 + x42111)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          
d/dt(x40111) = ac*Nf*x00110*x40011 + bbin*Nf*x00011*x40110 - dc*x40111 - ubc*x40111 - unbbin*x40111 - up*x40111 - agp*Nf*x00210*x40111 - ar*Nf*(x01010 + x02010)*x40111 +dg*x40311 + dr*(x41111 + x42111)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          
d/dt(x40200) =agp*x00200*x40000 - dg*x40200 - nl*x40200 - up*x40200 - gto*x40200 - ac*x00100*x40200 - ar*(x01000 + x02000)*x40200 + ne*x40210 + dc*x40300 +dr*(x41200 + x42200)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
d/dt(x40210) =agp*Nf*x00210*x40010 + nl*x40200 - dg*x40210 - ne*x40210 - up*x40210 - gto*x40210 - bbin*Nf*x00011*x40210 - ac*Nf*x00110*x40210 - ar*Nf*(x01010 + x02010)*x40210 - ar*Nf*(x01011 + x02011)*x40210 + ubc*x40211 + unbbin*x40211 + dc*x40310 +dr*(x41210 + x42210) + dr*(x41211 + x42211)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
d/dt(x40211) =agp*Nf*x00210*x40011 + bbin*Nf*x00011*x40210 - dg*x40211 - ubc*x40211 - unbbin*x40211 - up*x40211 - gto*x40211 - ac*Nf*x00110*x40211 - ar*Nf*(x01010 + x02010)*x40211 +dc*x40311 + dr*(x41211 + x42211)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
d/dt(x40300) =hto*x30300 + agp*x00200*x40100 + ac*x00100*x40200 - dc*x40300 - dg*x40300 - nl*x40300 - up*x40300 - gto*x40300 - ar*(x01000 + x02000)*x40300 + ne*x40310 +dr*(x41300 + x42300)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
d/dt(x40310) = agp*Nf*x00210*x40110 + ac*Nf*x00110*x40210 + nl*x40300 - dc*x40310 - dg*x40310 - ne*x40310 - up*x40310 - gto*x40310 - bbin*Nf*x00011*x40310 - ar*Nf*(x01010 + x02010)*x40310 - ar*Nf*(x01011 + x02011)*x40310 + ubc*x40311 +unbbin*x40311 + dr*(x41310 + x42310) + dr*(x41311 + x42311)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
d/dt(x40311) = agp*Nf*x00210*x40111 + ac*Nf*x00110*x40211 + bbin*Nf*x00011*x40310 - dc*x40311 - dg*x40311 - ubc*x40311 - unbbin*x40311 - up*x40311 - gto*x40311 - ar*Nf*(x01010 + x02010)*x40311 + dr*(x41311 + x42311)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
d/dt(x41000) =ar*x01000*x40000 - dr*x41000 - nl*x41000 - ac*x00100*x41000 - agp*x00200*x41000 + ne*x41010 + dc*x41100 + dg*x41200                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x41010) =ar*Nf*x01010*x40010 + nl*x41000 - dr*x41010 - ne*x41010 - bbin*Nf*x00011*x41010 - ac*Nf*x00110*x41010 - agp*Nf*x00210*x41010 + unbbin*x41011 + dc*x41110 + dg*x41210                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x41011) =ar*Nf*x01011*x40010 + ar*Nf*x01010*x40011 + bbin*Nf*x00011*x41010 - 2*dr*x41011 - unbbin*x41011 - ac*Nf*x00110*x41011 - agp*Nf*x00210*x41011 + dc*x41111 + dg*x41211                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x41100) =ar*x01000*x40100 + ac*x00100*x41000 - dc*x41100 - dr*x41100 - nl*x41100 - agp*x00200*x41100 + ne*x41110 + dg*x41300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x41110) =ar*Nf*x01010*x40110 + ac*Nf*x00110*x41010 + nl*x41100 - dc*x41110 - dr*x41110 - ne*x41110 - bbin*Nf*x00011*x41110 - agp*Nf*x00210*x41110 + unbbin*x41111 + dg*x41310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x41111) =ar*Nf*x01011*x40110 + ar*Nf*x01010*x40111 + ac*Nf*x00110*x41011 +bbin*Nf*x00011*x41110 - dc*x41111 - 2*dr*x41111 - unbbin*x41111 - agp*Nf*x00210*x41111 + dg*x41311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x41200) =ar*x01000*x40200 + agp*x00200*x41000 - dg*x41200 - dr*x41200 - nl*x41200 - gto*x41200 - ac*x00100*x41200 + ne*x41210 + dc*x41300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x41210) =ar*Nf*x01010*x40210 + agp*Nf*x00210*x41010 + nl*x41200 - dg*x41210 - dr*x41210 - ne*x41210 - gto*x41210 - bbin*Nf*x00011*x41210 - ac*Nf*x00110*x41210 +unbbin*x41211 + dc*x41310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x41211) =ar*Nf*x01011*x40210 + ar*Nf*x01010*x40211 + agp*Nf*x00210*x41011 +bbin*Nf*x00011*x41210 - dg*x41211 - 2*dr*x41211 - unbbin*x41211 - gto*x41211 - ac*Nf*x00110*x41211 + dc*x41311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x41300) =ar*x01000*x40300 + agp*x00200*x41100 + ac*x00100*x41200 - dc*x41300 - dg*x41300 - dr*x41300 - nl*x41300 - gto*x41300 + ne*x41310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x41310) =ar*Nf*x01010*x40310 + agp*Nf*x00210*x41110 + ac*Nf*x00110*x41210 + nl*x41300 - dc*x41310 - dg*x41310 - dr*x41310 - ne*x41310 - gto*x41310 - bbin*Nf*x00011*x41310 +unbbin*x41311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x41311) =ar*Nf*x01011*x40310 + ar*Nf*x01010*x40311 + agp*Nf*x00210*x41111 +ac*Nf*x00110*x41211 + bbin*Nf*x00011*x41310 - dc*x41311 - dg*x41311 - 2*dr*x41311 - unbbin*x41311 - gto*x41311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x42000) =ar*x02000*x40000 - dr*x42000 - nl*x42000 - ac*x00100*x42000 - agp*x00200*x42000 + ne*x42010 + dc*x42100 + dg*x42200                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x42010) =ar*Nf*x02010*x40010 + nl*x42000 - dr*x42010 - ne*x42010 - bbin*Nf*x00011*x42010 - ac*Nf*x00110*x42010 - agp*Nf*x00210*x42010 + unbbin*x42011 + dc*x42110 + dg*x42210                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x42011) =ar*Nf*x02011*x40010 + ar*Nf*x02010*x40011 + bbin*Nf*x00011*x42010 - 2*dr*x42011 - unbbin*x42011 - ac*Nf*x00110*x42011 - agp*Nf*x00210*x42011 + dc*x42111 + dg*x42211                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x42100) =ar*x02000*x40100 + ac*x00100*x42000 - dc*x42100 - dr*x42100 - nl*x42100 - agp*x00200*x42100 + ne*x42110 + dg*x42300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x42110) =ar*Nf*x02010*x40110 + ac*Nf*x00110*x42010 + nl*x42100 - dc*x42110 - dr*x42110 - ne*x42110 - bbin*Nf*x00011*x42110 - agp*Nf*x00210*x42110 + unbbin*x42111 + dg*x42310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x42111) =ar*Nf*x02011*x40110 + ar*Nf*x02010*x40111 + ac*Nf*x00110*x42011 +bbin*Nf*x00011*x42110 - dc*x42111 - 2*dr*x42111 - unbbin*x42111 - agp*Nf*x00210*x42111 + dg*x42311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x42200) =ar*x02000*x40200 + agp*x00200*x42000 - dg*x42200 - dr*x42200 - nl*x42200 - gto*x42200 - ac*x00100*x42200 + ne*x42210 + dc*x42300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x42210) =ar*Nf*x02010*x40210 + agp*Nf*x00210*x42010 + nl*x42200 - dg*x42210 - dr*x42210 - ne*x42210 - gto*x42210 - bbin*Nf*x00011*x42210 - ac*Nf*x00110*x42210 +unbbin*x42211 + dc*x42310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x42211) =ar*Nf*x02011*x40210 + ar*Nf*x02010*x40211 + agp*Nf*x00210*x42011 +bbin*Nf*x00011*x42210 - dg*x42211 - 2*dr*x42211 - unbbin*x42211 - gto*x42211 - ac*Nf*x00110*x42211 + dc*x42311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x42300) =ar*x02000*x40300 + agp*x00200*x42100 + ac*x00100*x42200 - dc*x42300 - dg*x42300 - dr*x42300 - nl*x42300 - gto*x42300 + ne*x42310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x42310) =ar*Nf*x02010*x40310 + agp*Nf*x00210*x42110 + ac*Nf*x00110*x42210 + nl*x42300 - dc*x42310 - dg*x42310 - dr*x42310 - ne*x42310 - gto*x42310 - bbin*Nf*x00011*x42310 +unbbin*x42311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x42311) =ar*Nf*x02011*x40310 + ar*Nf*x02010*x40311 + agp*Nf*x00210*x42111 +ac*Nf*x00110*x42211 + bbin*Nf*x00011*x42310 - dc*x42311 - dg*x42311 - 2*dr*x42311 - unbbin*x42311 - gto*x42311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x50000) =-(nl*x50000) - upu*x50000 - ac*x00100*x50000 - agp*x00200*x50000 - ar*(x01000 + x02000)*x50000 + ne*x50010 + dc*x50100 + dg*x50200 +dr*(x51000 + x52000)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
d/dt(x50010) =nl*x50000 - ne*x50010 - upu*x50010 - bbin*Nf*x00011*x50010 - ac*Nf*x00110*x50010 - agp*Nf*x00210*x50010 - ar*Nf*(x01010 + x02010)*x50010 - ar*Nf*(x01011 + x02011)*x50010 + ubc*x50011 + unbbin*x50011 + dc*x50110 + dg*x50210 +dr*(x51010 + x52010) + dr*(x51011 + x52011)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          
d/dt(x50011) =bbin*Nf*x00011*x50010 - ubc*x50011 - unbbin*x50011 - upu*x50011 - ac*Nf*x00110*x50011 - agp*Nf*x00210*x50011 - ar*Nf*(x01010 + x02010)*x50011 + dc*x50111 + dg*x50211 +dr*(x51011 + x52011)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          
d/dt(x50100) =ac*x00100*x50000 - dc*x50100 - hto*x50100 - nl*x50100 - upu*x50100 - agp*x00200*x50100 - ar*(x01000 + x02000)*x50100 + ne*x50110 + dg*x50300 +dr*(x51100 + x52100)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
d/dt(x50110) =ac*Nf*x00110*x50010 + nl*x50100 - dc*x50110 - hto*x50110 - ne*x50110 - upu*x50110 - bbin*Nf*x00011*x50110 - agp*Nf*x00210*x50110 - ar*Nf*(x01010 + x02010)*x50110 - ar*Nf*(x01011 + x02011)*x50110 + ubc*x50111 + unbbin*x50111 + dg*x50310 +dr*(x51110 + x52110) + dr*(x51111 + x52111)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
d/dt(x50111) =ac*Nf*x00110*x50011 + bbin*Nf*x00011*x50110 - dc*x50111 - hto*x50111 - ubc*x50111 - unbbin*x50111 - upu*x50111 - agp*Nf*x00210*x50111 - ar*Nf*(x01010 + x02010)*x50111 +dg*x50311 + dr*(x51111 + x52111)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
d/dt(x50200) =gto*x30200 + agp*x00200*x50000 - dg*x50200 - nl*x50200 - upu*x50200 - ac*x00100*x50200 - ar*(x01000 + x02000)*x50200 + ne*x50210 + dc*x50300 +dr*(x51200 + x52200)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
d/dt(x50210) =  agp*Nf*x00210*x50010 + nl*x50200 - dg*x50210 - ne*x50210 - upu*x50210 - bbin*Nf*x00011*x50210 - ac*Nf*x00110*x50210 - ar*Nf*(x01010 + x02010)*x50210 - ar*Nf*(x01011 + x02011)*x50210 + ubc*x50211 + unbbin*x50211 + dc*x50310 +dr*(x51210 + x52210) + dr*(x51211 + x52211)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
d/dt(x50211) =agp*Nf*x00210*x50011 + bbin*Nf*x00011*x50210 - dg*x50211 - ubc*x50211 - unbbin*x50211 - upu*x50211 - ac*Nf*x00110*x50211 - ar*Nf*(x01010 + x02010)*x50211 +dc*x50311 + dr*(x51211 + x52211)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          
d/dt(x50300) =gto*x30300 + agp*x00200*x50100 + ac*x00100*x50200 - dc*x50300 - dg*x50300 - hto*x50300 - nl*x50300 - upu*x50300 - ar*(x01000 + x02000)*x50300 + ne*x50310 +dr*(x51300 + x52300)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
d/dt(x50310) = agp*Nf*x00210*x50110 + ac*Nf*x00110*x50210 + nl*x50300 - dc*x50310 - dg*x50310 - hto*x50310 - ne*x50310 - upu*x50310 - bbin*Nf*x00011*x50310 - ar*Nf*(x01010 + x02010)*x50310 - ar*Nf*(x01011 + x02011)*x50310 + ubc*x50311 +unbbin*x50311 + dr*(x51310 + x52310) + dr*(x51311 + x52311)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
d/dt(x50311) = agp*Nf*x00210*x50111 + ac*Nf*x00110*x50211 + bbin*Nf*x00011*x50310 - dc*x50311 - dg*x50311 - hto*x50311 - ubc*x50311 - unbbin*x50311 - upu*x50311 - ar*Nf*(x01010 + x02010)*x50311 + dr*(x51311 + x52311)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           
d/dt(x51000) =ar*x01000*x50000 - dr*x51000 - nl*x51000 - ac*x00100*x51000 - agp*x00200*x51000 + ne*x51010 + dc*x51100 + dg*x51200                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x51010) =ar*Nf*x01010*x50010 + nl*x51000 - dr*x51010 - ne*x51010 - bbin*Nf*x00011*x51010 - ac*Nf*x00110*x51010 - agp*Nf*x00210*x51010 + unbbin*x51011 + dc*x51110 + dg*x51210                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x51011) =ar*Nf*x01011*x50010 + ar*Nf*x01010*x50011 + bbin*Nf*x00011*x51010 - 2*dr*x51011 - unbbin*x51011 - ac*Nf*x00110*x51011 - agp*Nf*x00210*x51011 + dc*x51111 + dg*x51211                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x51100) =ar*x01000*x50100 + ac*x00100*x51000 - dc*x51100 - dr*x51100 - nl*x51100 - agp*x00200*x51100 + ne*x51110 + dg*x51300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x51110) =ar*Nf*x01010*x50110 + ac*Nf*x00110*x51010 + nl*x51100 - dc*x51110 - dr*x51110 - ne*x51110 - bbin*Nf*x00011*x51110 - agp*Nf*x00210*x51110 + unbbin*x51111 + dg*x51310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x51111) =ar*Nf*x01011*x50110 + ar*Nf*x01010*x50111 + ac*Nf*x00110*x51011 +bbin*Nf*x00011*x51110 - dc*x51111 - 2*dr*x51111 - unbbin*x51111 - agp*Nf*x00210*x51111 + dg*x51311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x51200) =ar*x01000*x50200 + agp*x00200*x51000 - dg*x51200 - dr*x51200 - nl*x51200 - ac*x00100*x51200 + ne*x51210 + dc*x51300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x51210) = ar*Nf*x01010*x50210 + agp*Nf*x00210*x51010 + nl*x51200 - dg*x51210 - dr*x51210 - ne*x51210 - bbin*Nf*x00011*x51210 - ac*Nf*x00110*x51210 + unbbin*x51211 +dc*x51310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x51211) =ar*Nf*x01011*x50210 + ar*Nf*x01010*x50211 + agp*Nf*x00210*x51011 +bbin*Nf*x00011*x51210 - dg*x51211 - 2*dr*x51211 - unbbin*x51211 - ac*Nf*x00110*x51211 +dc*x51311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
d/dt(x51300) =ar*x01000*x50300 + agp*x00200*x51100 + ac*x00100*x51200 - dc*x51300 - dg*x51300 - dr*x51300 - nl*x51300 + ne*x51310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x51310) =ar*Nf*x01010*x50310 + agp*Nf*x00210*x51110 + ac*Nf*x00110*x51210 +nl*x51300 - dc*x51310 - dg*x51310 - dr*x51310 - ne*x51310 - bbin*Nf*x00011*x51310 +unbbin*x51311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
d/dt(x51311) = ar*Nf*x01011*x50310 + ar*Nf*x01010*x50311 + agp*Nf*x00210*x51111 +ac*Nf*x00110*x51211 + bbin*Nf*x00011*x51310 - dc*x51311 - dg*x51311 - 2*dr*x51311 - unbbin*x51311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x52000) =ar*x02000*x50000 - dr*x52000 - nl*x52000 - ac*x00100*x52000 - agp*x00200*x52000 + ne*x52010 + dc*x52100 + dg*x52200                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x52010) =ar*Nf*x02010*x50010 + nl*x52000 - dr*x52010 - ne*x52010 - bbin*Nf*x00011*x52010 - ac*Nf*x00110*x52010 - agp*Nf*x00210*x52010 + unbbin*x52011 + dc*x52110 + dg*x52210                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x52011) =ar*Nf*x02011*x50010 + ar*Nf*x02010*x50011 + bbin*Nf*x00011*x52010 - 2*dr*x52011 - unbbin*x52011 - ac*Nf*x00110*x52011 - agp*Nf*x00210*x52011 + dc*x52111 + dg*x52211                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x52100) =ar*x02000*x50100 + ac*x00100*x52000 - dc*x52100 - dr*x52100 - nl*x52100 - agp*x00200*x52100 + ne*x52110 + dg*x52300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x52110) =ar*Nf*x02010*x50110 + ac*Nf*x00110*x52010 + nl*x52100 - dc*x52110 - dr*x52110 - ne*x52110 - bbin*Nf*x00011*x52110 - agp*Nf*x00210*x52110 + unbbin*x52111 + dg*x52310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x52111) =ar*Nf*x02011*x50110 + ar*Nf*x02010*x50111 + ac*Nf*x00110*x52011 +bbin*Nf*x00011*x52110 - dc*x52111 - 2*dr*x52111 - unbbin*x52111 - agp*Nf*x00210*x52111 + dg*x52311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x52200) = ar*x02000*x50200 + agp*x00200*x52000 - dg*x52200 - dr*x52200 - nl*x52200 - ac*x00100*x52200 + ne*x52210 + dc*x52300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x52210) =ar*Nf*x02010*x50210 + agp*Nf*x00210*x52010 + nl*x52200 - dg*x52210 - dr*x52210 - ne*x52210 - bbin*Nf*x00011*x52210 - ac*Nf*x00110*x52210 + unbbin*x52211 +dc*x52310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x52211) =ar*Nf*x02011*x50210 + ar*Nf*x02010*x50211 + agp*Nf*x00210*x52011 +bbin*Nf*x00011*x52210 - dg*x52211 - 2*dr*x52211 - unbbin*x52211 - ac*Nf*x00110*x52211 +dc*x52311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
d/dt(x52300) = ar*x02000*x50300 + agp*x00200*x52100 + ac*x00100*x52200 - dc*x52300 - dg*x52300 - dr*x52300 - nl*x52300 + ne*x52310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x52310) = ar*Nf*x02010*x50310 + agp*Nf*x00210*x52110 + ac*Nf*x00110*x52210 +nl*x52300 - dc*x52310 - dg*x52310 - dr*x52310 - ne*x52310 - bbin*Nf*x00011*x52310 +unbbin*x52311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x52311) = ar*Nf*x02011*x50310 + ar*Nf*x02010*x50311 + agp*Nf*x00210*x52111 +ac*Nf*x00110*x52211 + bbin*Nf*x00011*x52310 - dc*x52311 - dg*x52311 - 2*dr*x52311 - unbbin*x52311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x60000) =-(nl*x60000) - up*x60000 - ac*x00100*x60000 - agp*x00200*x60000 - ar*(x01000 + x02000)*x60000 + ne*x60010 + dc*x60100 + dg*x60200 +dr*(x61000 + x62000)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
d/dt(x60010) =nl*x60000 - ne*x60010 - up*x60010 - bbin*Nf*x00011*x60010 - ac*Nf*x00110*x60010 - agp*Nf*x00210*x60010 - ar*Nf*(x01010 + x02010)*x60010 - ar*Nf*(x01011 + x02011)*x60010 + ubc*x60011 + unbbin*x60011 + dc*x60110 + dg*x60210 +dr*(x61010 + x62010) + dr*(x61011 + x62011)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           
d/dt(x60011) =bbin*Nf*x00011*x60010 - ubc*x60011 - unbbin*x60011 - up*x60011 - ac*Nf*x00110*x60011 - agp*Nf*x00210*x60011 - ar*Nf*(x01010 + x02010)*x60011 + dc*x60111 + dg*x60211 +dr*(x61011 + x62011)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           
d/dt(x60100) =hto*x50100 + ac*x00100*x60000 - dc*x60100 - nl*x60100 - up*x60100 - agp*x00200*x60100 - ar*(x01000 + x02000)*x60100 + ne*x60110 + dg*x60300 +dr*(x61100 + x62100)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
d/dt(x60110) =hto*x50110 + ac*Nf*x00110*x60010 + nl*x60100 - dc*x60110 - ne*x60110 - up*x60110 - bbin*Nf*x00011*x60110 - agp*Nf*x00210*x60110 - ar*Nf*(x01010 + x02010)*x60110 - ar*Nf*(x01011 + x02011)*x60110 + ubc*x60111 + unbbin*x60111 + dg*x60310 +dr*(x61110 + x62110) + dr*(x61111 + x62111)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
d/dt(x60111) =hto*x50111 + ac*Nf*x00110*x60011 + bbin*Nf*x00011*x60110 - dc*x60111 - ubc*x60111 - unbbin*x60111 - up*x60111 - agp*Nf*x00210*x60111 - ar*Nf*(x01010 + x02010)*x60111 +dg*x60311 + dr*(x61111 + x62111)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
d/dt(x60200) =gto*x40200 + agp*x00200*x60000 - dg*x60200 - nl*x60200 - up*x60200 - ac*x00100*x60200 - ar*(x01000 + x02000)*x60200 + ne*x60210 + dc*x60300 +dr*(x61200 + x62200)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
d/dt(x60210) =gto*x40210 + agp*Nf*x00210*x60010 + nl*x60200 - dg*x60210 - ne*x60210 - up*x60210 - bbin*Nf*x00011*x60210 - ac*Nf*x00110*x60210 - ar*Nf*(x01010 + x02010)*x60210 - ar*Nf*(x01011 + x02011)*x60210 + ubc*x60211 + unbbin*x60211 + dc*x60310 +dr*(x61210 + x62210) + dr*(x61211 + x62211)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
d/dt(x60211) =gto*x40211 + agp*Nf*x00210*x60011 + bbin*Nf*x00011*x60210 - dg*x60211 - ubc*x60211 - unbbin*x60211 - up*x60211 - ac*Nf*x00110*x60211 - ar*Nf*(x01010 + x02010)*x60211 +dc*x60311 + dr*(x61211 + x62211)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
d/dt(x60300) =gto*x40300 + hto*x50300 + agp*x00200*x60100 + ac*x00100*x60200 - dc*x60300 - dg*x60300 - nl*x60300 - up*x60300 - ar*(x01000 + x02000)*x60300 + ne*x60310 +dr*(x61300 + x62300)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
d/dt(x60310) =gto*x40310 + hto*x50310 + agp*Nf*x00210*x60110 + ac*Nf*x00110*x60210 + nl*x60300 - dc*x60310 - dg*x60310 - ne*x60310 - up*x60310 - bbin*Nf*x00011*x60310 - ar*Nf*(x01010 + x02010)*x60310 - ar*Nf*(x01011 + x02011)*x60310 + ubc*x60311 +unbbin*x60311 + dr*(x61310 + x62310) + dr*(x61311 + x62311)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x60311) =gto*x40311 + hto*x50311 + agp*Nf*x00210*x60111 + ac*Nf*x00110*x60211 +bbin*Nf*x00011*x60310 - dc*x60311 - dg*x60311 - ubc*x60311 - unbbin*x60311 - up*x60311 - ar*Nf*(x01010 + x02010)*x60311 + dr*(x61311 + x62311)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x61000) =ar*x01000*x60000 - dr*x61000 - nl*x61000 - ac*x00100*x61000 - agp*x00200*x61000 + ne*x61010 + dc*x61100 + dg*x61200                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x61010) =ar*Nf*x01010*x60010 + nl*x61000 - dr*x61010 - ne*x61010 - bbin*Nf*x00011*x61010 - ac*Nf*x00110*x61010 - agp*Nf*x00210*x61010 + unbbin*x61011 + dc*x61110 + dg*x61210                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x61011) =ar*Nf*x01011*x60010 + ar*Nf*x01010*x60011 + bbin*Nf*x00011*x61010 - 2*dr*x61011 - unbbin*x61011 - ac*Nf*x00110*x61011 - agp*Nf*x00210*x61011 + dc*x61111 + dg*x61211                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x61100) =ar*x01000*x60100 + ac*x00100*x61000 - dc*x61100 - dr*x61100 - nl*x61100 - agp*x00200*x61100 + ne*x61110 + dg*x61300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x61110) =ar*Nf*x01010*x60110 + ac*Nf*x00110*x61010 + nl*x61100 - dc*x61110 - dr*x61110 - ne*x61110 - bbin*Nf*x00011*x61110 - agp*Nf*x00210*x61110 + unbbin*x61111 + dg*x61310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x61111) =ar*Nf*x01011*x60110 + ar*Nf*x01010*x60111 + ac*Nf*x00110*x61011 +bbin*Nf*x00011*x61110 - dc*x61111 - 2*dr*x61111 - unbbin*x61111 - agp*Nf*x00210*x61111 + dg*x61311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x61200) =gto*x41200 + ar*x01000*x60200 + agp*x00200*x61000 - dg*x61200 - dr*x61200 - nl*x61200 - ac*x00100*x61200 + ne*x61210 + dc*x61300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x61210) =gto*x41210 + ar*Nf*x01010*x60210 + agp*Nf*x00210*x61010 + nl*x61200 - dg*x61210 - dr*x61210 - ne*x61210 - bbin*Nf*x00011*x61210 - ac*Nf*x00110*x61210 + unbbin*x61211 +dc*x61310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x61211) =gto*x41211 + ar*Nf*x01011*x60210 + ar*Nf*x01010*x60211 + agp*Nf*x00210*x61011 +bbin*Nf*x00011*x61210 - dg*x61211 - 2*dr*x61211 - unbbin*x61211 - ac*Nf*x00110*x61211 +dc*x61311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
d/dt(x61300) =gto*x41300 + ar*x01000*x60300 + agp*x00200*x61100 + ac*x00100*x61200 - dc*x61300 - dg*x61300 - dr*x61300 - nl*x61300 + ne*x61310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x61310) =gto*x41310 + ar*Nf*x01010*x60310 + agp*Nf*x00210*x61110 + ac*Nf*x00110*x61210 +nl*x61300 - dc*x61310 - dg*x61310 - dr*x61310 - ne*x61310 - bbin*Nf*x00011*x61310 +unbbin*x61311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
d/dt(x61311) =gto*x41311 + ar*Nf*x01011*x60310 + ar*Nf*x01010*x60311 + agp*Nf*x00210*x61111 +ac*Nf*x00110*x61211 + bbin*Nf*x00011*x61310 - dc*x61311 - dg*x61311 - 2*dr*x61311 - unbbin*x61311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x62000) =ar*x02000*x60000 - dr*x62000 - nl*x62000 - ac*x00100*x62000 - agp*x00200*x62000 + ne*x62010 + dc*x62100 + dg*x62200                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x62010) =ar*Nf*x02010*x60010 + nl*x62000 - dr*x62010 - ne*x62010 - bbin*Nf*x00011*x62010 - ac*Nf*x00110*x62010 - agp*Nf*x00210*x62010 + unbbin*x62011 + dc*x62110 + dg*x62210                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x62011) =ar*Nf*x02011*x60010 + ar*Nf*x02010*x60011 + bbin*Nf*x00011*x62010 - 2*dr*x62011 - unbbin*x62011 - ac*Nf*x00110*x62011 - agp*Nf*x00210*x62011 + dc*x62111 + dg*x62211                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x62100) =ar*x02000*x60100 + ac*x00100*x62000 - dc*x62100 - dr*x62100 - nl*x62100 - agp*x00200*x62100 + ne*x62110 + dg*x62300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x62110) =ar*Nf*x02010*x60110 + ac*Nf*x00110*x62010 + nl*x62100 - dc*x62110 - dr*x62110 - ne*x62110 - bbin*Nf*x00011*x62110 - agp*Nf*x00210*x62110 + unbbin*x62111 + dg*x62310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
d/dt(x62111) =ar*Nf*x02011*x60110 + ar*Nf*x02010*x60111 + ac*Nf*x00110*x62011 +bbin*Nf*x00011*x62110 - dc*x62111 - 2*dr*x62111 - unbbin*x62111 - agp*Nf*x00210*x62111 + dg*x62311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
d/dt(x62200) =gto*x42200 + ar*x02000*x60200 + agp*x00200*x62000 - dg*x62200 - dr*x62200 - nl*x62200 - ac*x00100*x62200 + ne*x62210 + dc*x62300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x62210) =gto*x42210 + ar*Nf*x02010*x60210 + agp*Nf*x00210*x62010 + nl*x62200 - dg*x62210 - dr*x62210 - ne*x62210 - bbin*Nf*x00011*x62210 - ac*Nf*x00110*x62210 + unbbin*x62211 +dc*x62310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x62211) =gto*x42211 + ar*Nf*x02011*x60210 + ar*Nf*x02010*x60211 + agp*Nf*x00210*x62011 +bbin*Nf*x00011*x62210 - dg*x62211 - 2*dr*x62211 - unbbin*x62211 - ac*Nf*x00110*x62211 +dc*x62311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
d/dt(x62300) =gto*x42300 + ar*x02000*x60300 + agp*x00200*x62100 + ac*x00100*x62200 - dc*x62300 - dg*x62300 - dr*x62300 - nl*x62300 + ne*x62310                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
d/dt(x62310) =gto*x42310 + ar*Nf*x02010*x60310 + agp*Nf*x00210*x62110 + ac*Nf*x00110*x62210 +nl*x62300 - dc*x62310 - dg*x62310 - dr*x62310 - ne*x62310 - bbin*Nf*x00011*x62310 +unbbin*x62311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
d/dt(x62311) =gto*x42311 + ar*Nf*x02011*x60310 + ar*Nf*x02010*x60311 + agp*Nf*x00210*x62111 +ac*Nf*x00110*x62211 + bbin*Nf*x00011*x62310 - dc*x62311 - dg*x62311 - 2*dr*x62311 - unbbin*x62311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
G(0)=0.210048                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
GR(0)=0.787955                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
gto(0)=1.13837                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
MnRo(0)=3.76085                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
McRo(0)=1.09547                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
MnRt(0)=12.3735                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
McRt(0)=4.47054                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
MnPo(0)=4.18947                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
McPo(0)=0.579108                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
MnPt(0)=8.34963                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
McPt(0)=1.41168                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
MnB(0)=35.1842                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
McB(0)=8.56043                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
B(0)=23.7741                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
GBb(0)=0.589308                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
GBRb(0)=0.404247                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
Cl(0)=48.0977                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
MnRev(0)=31.869                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
McRev(0)=161.144                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
cyrev(0)=123.95                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
revn(0)=467.844                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
cyrevg(0)=0.193289                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
revng(0)=58.2391                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
cyrevgp(0)=0.0397137                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
revngp(0)=19.0662                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
GB(0)=0.599965                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
GBR(0)=0.39359                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
cyrevp(0)=0.11012                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
revnp(0)=67.6175                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
BC(0)=6.5815                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
Gr(0)=0.416217                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
GrR(0)=0.558565                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
McNp(0)=0.25596                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
MnNp(0)=0.520367                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
Gc(0)=0.193886                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
GcR(0)=0.804404                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
x00001(0)=0.342061                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
x00011(0)=8.17977                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
x00100(0)=50.5794                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
x00110(0)=3.08486                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
x00200(0)=0.0985802                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x00210(0)=1.30846                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
x01000(0)=41.1658                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
x01010(0)=0.0156084                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x01011(0)=8.72089                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
x02000(0)=46.5721                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
x02010(0)=0.0031953                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x02011(0)=1.36799                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
x10000(0)=0.422898                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
x10100(0)=1.09254                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
x20000(0)=0.00181722                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x20010(0)=0.0000519698                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x20011(0)=0.00512223                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x20100(0)=0.106575                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
x20110(0)=0.000833073                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x20111(0)=0.0566636                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x21000(0)=0.00329174                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x21010(0)=0.000014996                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x21011(0)=0.00543134                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x21100(0)=0.0876909                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x21110(0)=0.000333744                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x21111(0)=0.0548927                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x22000(0)=0.00348546                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x22010(0)=0.0000154312                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x22011(0)=0.00510457                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x22100(0)=0.0938758                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x22110(0)=0.000353611                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x22111(0)=0.0535918                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x30000(0)=0.914365                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
x30100(0)=0.664039                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
x30200(0)=0.0226162                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x30300(0)=0.0241531                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x40000(0)=0.00356643                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x40010(0)=0.0000194469                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x40011(0)=0.000959363                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x40100(0)=0.258442                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
x40110(0)=0.000985386                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x40111(0)=0.0249991                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x40200(0)=0.000157103                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x40210(0)=0.000011242                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x40211(0)=0.00132337                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x40300(0)=0.0105914                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x40310(0)=0.00023062                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x40311(0)=0.025311                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
x41000(0)=0.00390694                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x41010(0)=0.0000134882                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x41011(0)=0.000809309                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x41100(0)=0.124712                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
x41110(0)=0.000416813                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x41111(0)=0.0148919                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x41200(0)=0.000188833                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x41210(0)=0.00000179979                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
x41211(0)=0.00109293                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x41300(0)=0.00584769                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x41310(0)=0.0000435837                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x41311(0)=0.0166107                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x42000(0)=0.00428454                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x42010(0)=0.0000147716                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x42011(0)=0.00086132                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x42100(0)=0.137696                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
x42110(0)=0.000459811                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x42111(0)=0.0160349                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x42200(0)=0.0002068                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x42210(0)=0.00000194018                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
x42211(0)=0.0011516                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x42300(0)=0.00644626                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x42310(0)=0.0000476349                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x42311(0)=0.0177021                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x50000(0)=0.00569806                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x50010(0)=0.000190832                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x50011(0)=0.0566834                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x50100(0)=0.0287569                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x50110(0)=0.000413832                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x50111(0)=0.0494651                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x50200(0)=0.00414505                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x50210(0)=0.000386223                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x50211(0)=0.124464                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
x50300(0)=0.00740796                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x50310(0)=0.000654378                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x50311(0)=0.104089                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
x51000(0)=0.00666974                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x51010(0)=0.0000254221                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x51011(0)=0.00484524                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x51100(0)=0.0828541                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x51110(0)=0.000294564                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x51111(0)=0.0323573                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x51200(0)=0.0012142                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x51210(0)=0.000011508                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x51211(0)=0.00976858                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x51300(0)=0.0087474                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x51310(0)=0.0000798494                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x51311(0)=0.0607878                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x52000(0)=0.00706385                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x52010(0)=0.0000264292                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x52011(0)=0.00442251                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x52100(0)=0.0852703                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x52110(0)=0.000300906                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x52111(0)=0.0301143                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x52200(0)=0.00132641                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x52210(0)=0.0000112497                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x52211(0)=0.00878991                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x52300(0)=0.0091505                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x52310(0)=0.0000778605                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x52311(0)=0.055586                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
x60000(0)=0.000792825                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x60010(0)=0.0000257804                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x60011(0)=0.00334582                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x60100(0)=0.0274125                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x60110(0)=0.000385498                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x60111(0)=0.0589937                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x60200(0)=0.000125512                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x60210(0)=0.0000499769                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x60211(0)=0.00735186                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x60300(0)=0.00550233                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x60310(0)=0.000631987                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x60311(0)=0.12768                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
x61000(0)=0.00252101                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x61010(0)=0.0000100779                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x61011(0)=0.00248103                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x61100(0)=0.0581233                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x61110(0)=0.000207429                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x61111(0)=0.0239137                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x61200(0)=0.000305307                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x61210(0)=0.00000490217                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
x61211(0)=0.00521949                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x61300(0)=0.00753943                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x61310(0)=0.0000639585                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x61311(0)=0.0471435                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x62000(0)=0.00259772                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x62010(0)=0.0000101709                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x62011(0)=0.00224807                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x62100(0)=0.0602479                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x62110(0)=0.000213449                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x62111(0)=0.0225783                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
x62200(0)=0.000320067                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
x62210(0)=0.00000460372                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
x62211(0)=0.004701                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
x62300(0)=0.00797183                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               
x62310(0)=0.0000635261                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
x62311(0)=0.0440794

********** MODEL PARAMETERS
trPo=25.9201                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
trPt=44.854                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
trRo=23.0747                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
trRt=39.9409                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
trB=46.1038                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
trRev=102.92299999999999                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           
trNp=0.3297485941817571                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
tlp=1.81031                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
tlr=5.03882                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
tlb=0.530436                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
tlrev=8.90744                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
tlc=4.645890                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
tlnp=1.25099                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
agp=1.3962                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         
dg=2.93521                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         
ac=0.0456572                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
dc=0.108072                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
ar=0.0235285                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
dr=0.605268                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
cbin=0.0454894                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
uncbin=7.27215                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
bbin=6.92686                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
unbbin=0.130196                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
cbbin=6.59924                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
uncbbin=0.304176                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
ag=0.162392                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
bin=6.97166                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
unbin=0.255032                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
binrev=0.0120525                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
unbinrev=10.9741                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
binr=6.15445                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
unbinr=2.91009                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
binc=0.280863                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
unbinc=0.00886752                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
binrevb=0.00626588                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
unbinrevb=5.30559                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
tmc=0.16426                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
tmcrev=9.2631                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
nl=0.643086                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
ne=0.0269078                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
nlrev=9.63702                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
nerev=0.0152514                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
lne=0.594609                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
nlbc=5.26501                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
hoo=0.527453                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
hto=2.45584                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
phos=0.291429                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
lono=0.205813                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
lont=0.396392                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
lta=0.607387                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
ltb=0.013                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          
trgto=0.644602                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
ugto=0.0625777                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
Nf=3.35063                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         
up=3.537                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           
uro=0.17491                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
urt=0.481895                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
umNp=0.369493                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
umPo=0.7669622                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
umPt=0.58891980                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
umRo=0.40342526000000006                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           
umRt=0.45554418                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
ub=0.018800209999999998                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
uc=0.025165115999999998                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
ubc=0.34882936000000003                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
upu=0.0700322                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
urev=1.6487614725107194                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
uprev=0.5173026723236928                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           
umB=0.7954019999999999                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
umRev=1.510194

********** MODEL VARIABLES


********** MODEL REACTIONS


********** MODEL FUNCTIONS


********** MODEL EVENTS


********** MODEL MATLAB FUNCTIONS


Loading data, please wait...