##########################################################################
# Detailed Model of Mammalian Circadian Clock
# Generated: 11-Sep-2012 20:26:23
# Generated by Jae Kyoung Kim and Daniel Forger by using SBtoolbox2.  
# This file can be used by XPPAUT for simulation and
# bifurcation analysis.
##########################################################################

########################################################
# DIFFERENTIAL EQUATIONS
########################################################

#Since XPP does not distinguish GR and Gr, Gr is changed to Grn 

GR'=-(unbin*GR)+bin*(1-G-GR)*(x01011+x02011)
G'=-(unbin*G)+bin*(1-G-GR)*x00011
GrR'=-(unbinr*GrR)+binr*(1-Grn-GrR)*(x01011+x02011)
Grn'=-(unbinr*Grn)+binr*(1-Grn-GrR)*x00011
GcR'=-(unbinc*GcR)+binc*(1-Gc-GcR)*(x01011+x02011)
Gc'=-(unbinc*Gc)+binc*(1-Gc-GcR)*x00011
GBR'=-(unbinrev*GBR)+binrev*GB*(revn+revng+revngp+revnp)
GB'=unbinrev*GBR-binrev*GB*(revn+revng+revngp+revnp)
GBRb'=-(unbinrevb*GBRb)+binrevb*GBb*(revn+revng+revngp+revnp)
GBb'=unbinrevb*GBRb-binrevb*GBb*(revn+revng+revngp+revnp)
MnPo'=trPo*G-tmc*MnPo-umPo*MnPo
McPo'=-(umPo*McPo)+tmc*MnPo
MnPt'=trPt*G-tmc*MnPt-umPt*MnPt
McPt'=-(umPt*McPt)+tmc*MnPt
MnRt'=trRt*Gc-tmc*MnRt-umRt*MnRt
McRt'=-(umRt*McRt)+tmc*MnRt
MnRev'=-(tmcrev*MnRev)-umRev*MnRev+trRev*Grn*x00011
McRev'=-(umRev*McRev)+tmcrev*MnRev
MnRo'=trRo*G*GB-tmc*MnRo-umRo*MnRo
McRo'=-(umRo*McRo)+tmc*MnRo
MnB'=trB*GBb-tmc*MnB-umB*MnB
McB'=-(umB*McB)+tmc*MnB
MnNp'=trNp*GB-tmc*MnNp-umNp*MnNp
McNp'=-(umNp*McNp)+tmc*MnNp
B'=-(ub*B)+uncbin*BC-cbin*B*Cl+tlb*McB
Cl'=tlc+uncbin*BC-uc*Cl-cbin*B*Cl+tlnp*McNp
BC'=-(phos*BC)-ubc*BC-uncbin*BC+cbin*B*Cl
cyrev'=-((nlrev+urev)*cyrev)+dg*cyrevg+tlrev*McRev+nerev*revn-ag*cyrev*x00200
revn'=nlrev*cyrev+(-nerev-urev)*revn+dg*revng-ag*Nf*revn*x00210
cyrevg'=-(cyrevg*(dg+nlrev+urev+gto))+nerev*revng+ag*cyrev*x00200
revng'=nlrev*cyrevg-(dg+nerev+urev+gto)*revng+ag*Nf*revn*x00210
cyrevgp'=-((dg+nlrev+uprev)*cyrevgp)+cyrevg*gto+nerev*revngp
revngp'=nlrev*cyrevgp+gto*revng-(dg+nerev+uprev)*revngp
cyrevp'=dg*cyrevgp-(nlrev+uprev)*cyrevp+nerev*revnp
revnp'=nlrev*cyrevp+dg*revngp-(nerev+uprev)*revnp
gto'=trgto*G*GB-ugto*gto
x00001'=phos*BC-nlbc*x00001-ubc*x00001
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)
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)
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)
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)
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)
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)
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)
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)
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)
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)
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)
x10000'=tlp*McPo-upu*x10000-ac*x00100*x10000+dc*x10100
x10100'=ac*x00100*x10000-dc*x10100-hoo*x10100-upu*x10100
x20000'=-(nl*x20000)-up*x20000-ac*x00100*x20000-ar*(x01000+x02000)*x20000+ne*x20010+dc*x20100+dr*(x21000+x22000)
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)
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)
x20100'=hoo*x10100+ac*x00100*x20000-dc*x20100-nl*x20100-up*x20100-ar*(x01000+x02000)*x20100+ne*x20110+dr*(x21100+x22100)
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)
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)
x21000'=ar*x01000*x20000-dr*x21000-nl*x21000-ac*x00100*x21000+ne*x21010+dc*x21100
x21010'=ar*Nf*x01010*x20010+nl*x21000-dr*x21010-ne*x21010-bbin*Nf*x00011*x21010-ac*Nf*x00110*x21010+unbbin*x21011+dc*x21110
x21011'=ar*Nf*x01011*x20010+ar*Nf*x01010*x20011+bbin*Nf*x00011*x21010-2*dr*x21011-unbbin*x21011-ac*Nf*x00110*x21011+dc*x21111
x21100'=ar*x01000*x20100+ac*x00100*x21000-dc*x21100-dr*x21100-nl*x21100+ne*x21110
x21110'=ar*Nf*x01010*x20110+ac*Nf*x00110*x21010+nl*x21100-dc*x21110-dr*x21110-ne*x21110-bbin*Nf*x00011*x21110+unbbin*x21111
x21111'=ar*Nf*x01011*x20110+ar*Nf*x01010*x20111+ac*Nf*x00110*x21011+bbin*Nf*x00011*x21110-dc*x21111-2*dr*x21111-unbbin*x21111
x22000'=ar*x02000*x20000-dr*x22000-nl*x22000-ac*x00100*x22000+ne*x22010+dc*x22100
x22010'=ar*Nf*x02010*x20010+nl*x22000-dr*x22010-ne*x22010-bbin*Nf*x00011*x22010-ac*Nf*x00110*x22010+unbbin*x22011+dc*x22110
x22011'=ar*Nf*x02011*x20010+ar*Nf*x02010*x20011+bbin*Nf*x00011*x22010-2*dr*x22011-unbbin*x22011-ac*Nf*x00110*x22011+dc*x22111
x22100'=ar*x02000*x20100+ac*x00100*x22000-dc*x22100-dr*x22100-nl*x22100+ne*x22110
x22110'=ar*Nf*x02010*x20110+ac*Nf*x00110*x22010+nl*x22100-dc*x22110-dr*x22110-ne*x22110-bbin*Nf*x00011*x22110+unbbin*x22111
x22111'=ar*Nf*x02011*x20110+ar*Nf*x02010*x20111+ac*Nf*x00110*x22011+bbin*Nf*x00011*x22110-dc*x22111-2*dr*x22111-unbbin*x22111
x30000'=tlp*McPt-upu*x30000-ac*x00100*x30000-agp*x00200*x30000+dc*x30100+dg*x30200
x30100'=ac*x00100*x30000-dc*x30100-hto*x30100-upu*x30100-agp*x00200*x30100+dg*x30300
x30200'=agp*x00200*x30000-dg*x30200-upu*x30200-gto*x30200-ac*x00100*x30200+dc*x30300
x30300'=agp*x00200*x30100+ac*x00100*x30200-dc*x30300-dg*x30300-hto*x30300-upu*x30300-gto*x30300
x40000'=-(nl*x40000)-up*x40000-ac*x00100*x40000-agp*x00200*x40000-ar*(x01000+x02000)*x40000+ne*x40010+dc*x40100+dg*x40200+dr*(x41000+x42000)
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)
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)
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)
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)
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)
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)
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)
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)
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)
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)
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)
x41000'=ar*x01000*x40000-dr*x41000-nl*x41000-ac*x00100*x41000-agp*x00200*x41000+ne*x41010+dc*x41100+dg*x41200
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
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
x41100'=ar*x01000*x40100+ac*x00100*x41000-dc*x41100-dr*x41100-nl*x41100-agp*x00200*x41100+ne*x41110+dg*x41300
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
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
x41200'=ar*x01000*x40200+agp*x00200*x41000-dg*x41200-dr*x41200-nl*x41200-gto*x41200-ac*x00100*x41200+ne*x41210+dc*x41300
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
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
x41300'=ar*x01000*x40300+agp*x00200*x41100+ac*x00100*x41200-dc*x41300-dg*x41300-dr*x41300-nl*x41300-gto*x41300+ne*x41310
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
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
x42000'=ar*x02000*x40000-dr*x42000-nl*x42000-ac*x00100*x42000-agp*x00200*x42000+ne*x42010+dc*x42100+dg*x42200
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
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
x42100'=ar*x02000*x40100+ac*x00100*x42000-dc*x42100-dr*x42100-nl*x42100-agp*x00200*x42100+ne*x42110+dg*x42300
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
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
x42200'=ar*x02000*x40200+agp*x00200*x42000-dg*x42200-dr*x42200-nl*x42200-gto*x42200-ac*x00100*x42200+ne*x42210+dc*x42300
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
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
x42300'=ar*x02000*x40300+agp*x00200*x42100+ac*x00100*x42200-dc*x42300-dg*x42300-dr*x42300-nl*x42300-gto*x42300+ne*x42310
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
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
x50000'=-(nl*x50000)-upu*x50000-ac*x00100*x50000-agp*x00200*x50000-ar*(x01000+x02000)*x50000+ne*x50010+dc*x50100+dg*x50200+dr*(x51000+x52000)
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)
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)
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)
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)
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)
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)
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)
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)
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)
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)
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)
x51000'=ar*x01000*x50000-dr*x51000-nl*x51000-ac*x00100*x51000-agp*x00200*x51000+ne*x51010+dc*x51100+dg*x51200
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
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
x51100'=ar*x01000*x50100+ac*x00100*x51000-dc*x51100-dr*x51100-nl*x51100-agp*x00200*x51100+ne*x51110+dg*x51300
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
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
x51200'=ar*x01000*x50200+agp*x00200*x51000-dg*x51200-dr*x51200-nl*x51200-ac*x00100*x51200+ne*x51210+dc*x51300
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
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
x51300'=ar*x01000*x50300+agp*x00200*x51100+ac*x00100*x51200-dc*x51300-dg*x51300-dr*x51300-nl*x51300+ne*x51310
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
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
x52000'=ar*x02000*x50000-dr*x52000-nl*x52000-ac*x00100*x52000-agp*x00200*x52000+ne*x52010+dc*x52100+dg*x52200
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
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
x52100'=ar*x02000*x50100+ac*x00100*x52000-dc*x52100-dr*x52100-nl*x52100-agp*x00200*x52100+ne*x52110+dg*x52300
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
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
x52200'=ar*x02000*x50200+agp*x00200*x52000-dg*x52200-dr*x52200-nl*x52200-ac*x00100*x52200+ne*x52210+dc*x52300
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
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
x52300'=ar*x02000*x50300+agp*x00200*x52100+ac*x00100*x52200-dc*x52300-dg*x52300-dr*x52300-nl*x52300+ne*x52310
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
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
x60000'=-(nl*x60000)-up*x60000-ac*x00100*x60000-agp*x00200*x60000-ar*(x01000+x02000)*x60000+ne*x60010+dc*x60100+dg*x60200+dr*(x61000+x62000)
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)
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)
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)
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)
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)
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)
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)
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)
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)
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)
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)
x61000'=ar*x01000*x60000-dr*x61000-nl*x61000-ac*x00100*x61000-agp*x00200*x61000+ne*x61010+dc*x61100+dg*x61200
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
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
x61100'=ar*x01000*x60100+ac*x00100*x61000-dc*x61100-dr*x61100-nl*x61100-agp*x00200*x61100+ne*x61110+dg*x61300
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
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
x61200'=gto*x41200+ar*x01000*x60200+agp*x00200*x61000-dg*x61200-dr*x61200-nl*x61200-ac*x00100*x61200+ne*x61210+dc*x61300
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
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
x61300'=gto*x41300+ar*x01000*x60300+agp*x00200*x61100+ac*x00100*x61200-dc*x61300-dg*x61300-dr*x61300-nl*x61300+ne*x61310
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
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
x62000'=ar*x02000*x60000-dr*x62000-nl*x62000-ac*x00100*x62000-agp*x00200*x62000+ne*x62010+dc*x62100+dg*x62200
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
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
x62100'=ar*x02000*x60100+ac*x00100*x62000-dc*x62100-dr*x62100-nl*x62100-agp*x00200*x62100+ne*x62110+dg*x62300
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
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
x62200'=gto*x42200+ar*x02000*x60200+agp*x00200*x62000-dg*x62200-dr*x62200-nl*x62200-ac*x00100*x62200+ne*x62210+dc*x62300
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
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
x62300'=gto*x42300+ar*x02000*x60300+agp*x00200*x62100+ac*x00100*x62200-dc*x62300-dg*x62300-dr*x62300-nl*x62300+ne*x62310
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
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

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

########################################################
# FUNCTIONS
########################################################
power(x,y)=x^y

########################################################
# INITIAL CONDITIONS
########################################################
GR(0)=0.787955
G(0)=0.210048
GrR(0)=0.558565
Grn(0)=0.416217
GcR(0)=0.804404
Gc(0)=0.193886
GBR(0)=0.39359
GB(0)=0.599965
GBRb(0)=0.404247
GBb(0)=0.589308
MnPo(0)=4.18947
McPo(0)=0.579108
MnPt(0)=8.34963
McPt(0)=1.41168
MnRt(0)=12.3735
McRt(0)=4.47054
MnRev(0)=31.869
McRev(0)=161.144
MnRo(0)=3.76085
McRo(0)=1.09547
MnB(0)=35.1842
McB(0)=8.56043
MnNp(0)=0.520367
McNp(0)=0.25596
B(0)=23.7741
Cl(0)=48.0977
BC(0)=6.5815
cyrev(0)=123.95
revn(0)=467.844
cyrevg(0)=0.193289
revng(0)=58.2391
cyrevgp(0)=0.0397137
revngp(0)=19.0662
cyrevp(0)=0.11012
revnp(0)=67.6175
gto(0)=1.13837
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)=5.19698e-05
x20011(0)=0.00512223
x20100(0)=0.106575
x20110(0)=0.000833073
x20111(0)=0.0566636
x21000(0)=0.00329174
x21010(0)=1.4996e-05
x21011(0)=0.00543134
x21100(0)=0.0876909
x21110(0)=0.000333744
x21111(0)=0.0548927
x22000(0)=0.00348546
x22010(0)=1.54312e-05
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)=1.94469e-05
x40011(0)=0.000959363
x40100(0)=0.258442
x40110(0)=0.000985386
x40111(0)=0.0249991
x40200(0)=0.000157103
x40210(0)=1.1242e-05
x40211(0)=0.00132337
x40300(0)=0.0105914
x40310(0)=0.00023062
x40311(0)=0.025311
x41000(0)=0.00390694
x41010(0)=1.34882e-05
x41011(0)=0.000809309
x41100(0)=0.124712
x41110(0)=0.000416813
x41111(0)=0.0148919
x41200(0)=0.000188833
x41210(0)=1.79979e-06
x41211(0)=0.00109293
x41300(0)=0.00584769
x41310(0)=4.35837e-05
x41311(0)=0.0166107
x42000(0)=0.00428454
x42010(0)=1.47716e-05
x42011(0)=0.00086132
x42100(0)=0.137696
x42110(0)=0.000459811
x42111(0)=0.0160349
x42200(0)=0.0002068
x42210(0)=1.94018e-06
x42211(0)=0.0011516
x42300(0)=0.00644626
x42310(0)=4.76349e-05
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)=2.54221e-05
x51011(0)=0.00484524
x51100(0)=0.0828541
x51110(0)=0.000294564
x51111(0)=0.0323573
x51200(0)=0.0012142
x51210(0)=1.1508e-05
x51211(0)=0.00976858
x51300(0)=0.0087474
x51310(0)=7.98494e-05
x51311(0)=0.0607878
x52000(0)=0.00706385
x52010(0)=2.64292e-05
x52011(0)=0.00442251
x52100(0)=0.0852703
x52110(0)=0.000300906
x52111(0)=0.0301143
x52200(0)=0.00132641
x52210(0)=1.12497e-05
x52211(0)=0.00878991
x52300(0)=0.0091505
x52310(0)=7.78605e-05
x52311(0)=0.055586
x60000(0)=0.000792825
x60010(0)=2.57804e-05
x60011(0)=0.00334582
x60100(0)=0.0274125
x60110(0)=0.000385498
x60111(0)=0.0589937
x60200(0)=0.000125512
x60210(0)=4.99769e-05
x60211(0)=0.00735186
x60300(0)=0.00550233
x60310(0)=0.000631987
x60311(0)=0.12768
x61000(0)=0.00252101
x61010(0)=1.00779e-05
x61011(0)=0.00248103
x61100(0)=0.0581233
x61110(0)=0.000207429
x61111(0)=0.0239137
x61200(0)=0.000305307
x61210(0)=4.90217e-06
x61211(0)=0.00521949
x61300(0)=0.00753943
x61310(0)=6.39585e-05
x61311(0)=0.0471435
x62000(0)=0.00259772
x62010(0)=1.01709e-05
x62011(0)=0.00224807
x62100(0)=0.0602479
x62110(0)=0.000213449
x62111(0)=0.0225783
x62200(0)=0.000320067
x62210(0)=4.60372e-06
x62211(0)=0.004701
x62300(0)=0.00797183
x62310(0)=6.35261e-05
x62311(0)=0.0440794

########################################################
# INTEGRATOR SETTINGS AND DONE
########################################################
@ method=stiff
@ bounds=100000
@ maxstor=20000
done