In [1]:

    

import pandas as pd
import numpy as np
import re


    

In [2]:

    

species_list = ['ADIO', 'ADIOM2', 'ALD3', 'C10H2', 'C10H2M2', 'C10H2M4', 'C2H6', 'C3H4O', 'C3H4O2', 'C3H6', 'C3H6O2', 'C3H8O2', 
                'CH2CO', 'CH3CHO', 'CH3OH', 'CH4', 'CHAR', 'CO', 'CO2', 'COUMARYL', 'ETOH', 'H2', 'H2O', 'KET', 'KETD', 'KETDM2', 'KETM2', 
                'LIG', 'LIGC', 'LIGH', 'LIGM2', 'LIGO', 'MGUAI', 'OH', 'PADIO', 'PADIOM2', 'PC2H2', 'PCH2OH', 'PCH2P', 'PCH3', 'PCHO', 
                'PCHOHP', 'PCHP2', 'PCOH', 'PCOHP2', 'PCOS', 'PFET3', 'PFET3M2', 'PH2', 'PHENOL', 'PKETM2', 'PLIG', 'PLIGC', 'PLIGH', 
                'PLIGM2', 'PLIGO', 'PRADIO', 'PRADIOM2', 'PRFET3', 'PRFET3M2', 'PRKETM2', 'PRLIGH', 'PRLIGH2', 'PRLIGM2A', 'RADIO', 
                'RADIOM2', 'RC3H3O', 'RC3H5O2', 'RC3H7O2', 'RCH3', 'RCH3O', 'RKET', 'RKETM2', 'RLIGA', 'RLIGB', 'RLIGH', 'RLIGM2A', 
                'RLIGM2B', 'RMGUAI', 'RPHENOL', 'RPHENOX', 'RPHENOXM2', 'SYNAPYL', 'VADIO', 'VADIOM2', 'VCOUMARYL', 'VKET', 'VKETD', 
                'VKETDM2', 'VKETM2', 'VMGUAI', 'VPHENOL', 'VSYNAPYL']


    

In [3]:

    

reaction_file = open('reaction_list','r').read().replace(' ', '')


    

In [4]:

    

reaction_file


    

    
Out[4]:





'/*ODEforPLIGC=>PLIGC=0*/\ndydt[PLIGC]=-A30*T**n30*np.exp(-1*E30/R/T)*PLIGC-A105*T**n105*np.exp(-1*E105/R/T)*PLIGC;\n/*ODEforPLIGH=>PLIGH=1*/\ndydt[PLIGH]=-A153*T**n153*np.exp(-1*E153/R/T)*RADIO*PLIGH-A137*T**n137*np.exp(-1*E137/R/T)*PLIGH*PRFET3M2+A201*T**n201*np.exp(-1*E201/R/T)*PRLIGH*LIGM2-A133*T**n133*np.exp(-1*E133/R/T)*PLIGH*RC3H5O2-A142*T**n142*np.exp(-1*E142/R/T)*PRKETM2*PLIGH-\nA135*T**n135*np.exp(-1*E135/R/T)*RC3H7O2*PLIGH+A306*T**n306*np.exp(-1*E306/R/T)*LIG*PRLIGH+A327*T**n327*np.exp(-1*E327/R/T)*LIG*PRLIGH-A104*T**n104*np.exp(-1*E104/R/T)*PLIGH-A141*T**n141*np.exp(-1*E141/R/T)*PLIGH*RCH3-A139*T**n139*np.exp(-1*E139/R/T)*RLIGM2B*PLIGH+\nA243*T**n243*np.exp(-1*E243/R/T)*PRLIGH*ADIOM2+A180*T**n180*np.exp(-1*E180/R/T)*PRLIGH*LIGM2+A285*T**n285*np.exp(-1*E285/R/T)*PRLIGH*C10H2-A145*T**n145*np.exp(-1*E145/R/T)*RC3H3O*PLIGH-A143*T**n143*np.exp(-1*E143/R/T)*PLIGH*RKET-A147*T**n147*np.exp(-1*E147/R/T)*PLIGH*RLIGA+\nA159*T**n159*np.exp(-1*E159/R/T)*PLIGM2*PRLIGH-A152*T**n152*np.exp(-1*E152/R/T)*PLIGH*RPHENOL-A151*T**n151*np.exp(-1*E151/R/T)*PLIGH*RCH3O-A144*T**n144*np.exp(-1*E144/R/T)*PLIGH*PRADIO-A150*T**n150*np.exp(-1*E150/R/T)*PLIGH*OH-A148*T**n148*np.exp(-1*E148/R/T)*PRADIOM2*PLIGH-\nA146*T**n146*np.exp(-1*E146/R/T)*PLIGH*RLIGB-A136*T**n136*np.exp(-1*E136/R/T)*PLIGH*RADIOM2+A348*T**n348*np.exp(-1*E348/R/T)*PFET3*PRLIGH-A140*T**n140*np.exp(-1*E140/R/T)*RLIGM2A*PLIGH-A0*T**n0*np.exp(-1*E0/R/T)*PLIGH-A134*T**n134*np.exp(-1*E134/R/T)*PLIGH*PRFET3+\nA222*T**n222*np.exp(-1*E222/R/T)*PRLIGH*PFET3M2+A264*T**n264*np.exp(-1*E264/R/T)*KETM2*PRLIGH+A390*T**n390*np.exp(-1*E390/R/T)*PRLIGH*KET+A117*T**n117*np.exp(-1*E117/R/T)*PRLIGH*LIGH-A149*T**n149*np.exp(-1*E149/R/T)*RMGUAI*PLIGH+A369*T**n369*np.exp(-1*E369/R/T)*PRLIGH*ADIO;\n/*ODEforPLIGO=>PLIGO=2*/\ndydt[PLIGO]=-A31*T**n31*np.exp(-1*E31/R/T)*PLIGO-A106*T**n106*np.exp(-1*E106/R/T)*PLIGO;\n/*ODEforPRLIGH=>PRLIGH=3*/\ndydt[PRLIGH]=-A243*T**n243*np.exp(-1*E243/R/T)*PRLIGH*ADIOM2-A390*T**n390*np.exp(-1*E390/R/T)*PRLIGH*KET-A180*T**n180*np.exp(-1*E180/R/T)*PRLIGH*LIGM2-A285*T**n285*np.exp(-1*E285/R/T)*PRLIGH*C10H2-A201*T**n201*np.exp(-1*E201/R/T)*PRLIGH*LIGM2-\nA348*T**n348*np.exp(-1*E348/R/T)*PFET3*PRLIGH-A306*T**n306*np.exp(-1*E306/R/T)*LIG*PRLIGH-A327*T**n327*np.exp(-1*E327/R/T)*LIG*PRLIGH+A0*T**n0*np.exp(-1*E0/R/T)*PLIGH-A159*T**n159*np.exp(-1*E159/R/T)*PLIGM2*PRLIGH-A117*T**n117*np.exp(-1*E117/R/T)*PRLIGH*LIGH-\nA138*T**n138*np.exp(-1*E138/R/T)*PLIGH*PRLIGH-A222*T**n222*np.exp(-1*E222/R/T)*PRLIGH*PFET3M2-A264*T**n264*np.exp(-1*E264/R/T)*KETM2*PRLIGH-A369*T**n369*np.exp(-1*E369/R/T)*PRLIGH*ADIO;\n/*ODEforLIGH=>LIGH=4*/\ndydt[LIGH]=-A123*T**n123*np.exp(-1*E123/R/T)*LIGH*PRADIO-A132*T**n132*np.exp(-1*E132/R/T)*RADIO*LIGH-A129*T**n129*np.exp(-1*E129/R/T)*LIGH*OH-A115*T**n115*np.exp(-1*E115/R/T)*LIGH*RADIOM2-A122*T**n122*np.exp(-1*E122/R/T)*LIGH*RKET-\nA127*T**n127*np.exp(-1*E127/R/T)*PRADIOM2*LIGH+A104*T**n104*np.exp(-1*E104/R/T)*PLIGH-A120*T**n120*np.exp(-1*E120/R/T)*RCH3*LIGH-A114*T**n114*np.exp(-1*E114/R/T)*RC3H7O2*LIGH-A126*T**n126*np.exp(-1*E126/R/T)*LIGH*RLIGA-A116*T**n116*np.exp(-1*E116/R/T)*PRFET3M2*LIGH-\nA118*T**n118*np.exp(-1*E118/R/T)*RLIGM2B*LIGH-A121*T**n121*np.exp(-1*E121/R/T)*PRKETM2*LIGH-A124*T**n124*np.exp(-1*E124/R/T)*RC3H3O*LIGH-A112*T**n112*np.exp(-1*E112/R/T)*LIGH*RC3H5O2-A125*T**n125*np.exp(-1*E125/R/T)*RLIGB*LIGH-A113*T**n113*np.exp(-1*E113/R/T)*LIGH*PRFET3-\nA119*T**n119*np.exp(-1*E119/R/T)*RLIGM2A*LIGH-A128*T**n128*np.exp(-1*E128/R/T)*RMGUAI*LIGH-A1*T**n1*np.exp(-1*E1/R/T)*LIGH-A131*T**n131*np.exp(-1*E131/R/T)*LIGH*RPHENOL-A117*T**n117*np.exp(-1*E117/R/T)*PRLIGH*LIGH-A130*T**n130*np.exp(-1*E130/R/T)*RCH3O*LIGH;\n/*ODEforC3H6=>C3H6=5*/\ndydt[C3H6]=+A1*T**n1*np.exp(-1*E1/R/T)*LIGH;\n/*ODEforOH=>OH=6*/\ndydt[OH]=+A99*T**n99*np.exp(-1*E99/R/T)*PCHOHP+A98*T**n98*np.exp(-1*E98/R/T)*PCH2OH+A26*T**n26*np.exp(-1*E26/R/T)*RKET-A234*T**n234*np.exp(-1*E234/R/T)*OH*PFET3M2-A339*T**n339*np.exp(-1*E339/R/T)*LIG*OH-\nA129*T**n129*np.exp(-1*E129/R/T)*LIGH*OH-A192*T**n192*np.exp(-1*E192/R/T)*LIGM2*OH-A297*T**n297*np.exp(-1*E297/R/T)*C10H2*OH+A100*T**n100*np.exp(-1*E100/R/T)*PCOHP2-A318*T**n318*np.exp(-1*E318/R/T)*LIG*OH-A276*T**n276*np.exp(-1*E276/R/T)*KETM2*OH+\nA9*T**n9*np.exp(-1*E9/R/T)*PRKETM2-A171*T**n171*np.exp(-1*E171/R/T)*PLIGM2*OH-A213*T**n213*np.exp(-1*E213/R/T)*LIGM2*OH-A69*T**n69*np.exp(-1*E69/R/T)*RCH3*OH+A25*T**n25*np.exp(-1*E25/R/T)*RADIO+A23*T**n23*np.exp(-1*E23/R/T)*RADIOM2-\nA381*T**n381*np.exp(-1*E381/R/T)*ADIO*OH+A24*T**n24*np.exp(-1*E24/R/T)*RKETM2+A1*T**n1*np.exp(-1*E1/R/T)*LIGH-A360*T**n360*np.exp(-1*E360/R/T)*PFET3*OH-A402*T**n402*np.exp(-1*E402/R/T)*KET*OH-A255*T**n255*np.exp(-1*E255/R/T)*OH*ADIOM2-\nA150*T**n150*np.exp(-1*E150/R/T)*PLIGH*OH;\n/*ODEforRLIGM2A=>RLIGM2A=7*/\ndydt[RLIGM2A]=+A191*T**n191*np.exp(-1*E191/R/T)*RMGUAI*LIGM2-A308*T**n308*np.exp(-1*E308/R/T)*LIG*RLIGM2A-A287*T**n287*np.exp(-1*E287/R/T)*RLIGM2A*C10H2+A194*T**n194*np.exp(-1*E194/R/T)*LIGM2*RPHENOL+A189*T**n189*np.exp(-1*E189/R/T)*RLIGA*LIGM2\n-A392*T**n392*np.exp(-1*E392/R/T)*RLIGM2A*KET+A180*T**n180*np.exp(-1*E180/R/T)*PRLIGH*LIGM2-A350*T**n350*np.exp(-1*E350/R/T)*RLIGM2A*PFET3+A185*T**n185*np.exp(-1*E185/R/T)*RKET*LIGM2-A119*T**n119*np.exp(-1*E119/R/T)*RLIGM2A*LIGH-A371*T**n371*np.exp(-1*E371/R/T)*RLIGM2A*ADIO+\nA1*T**n1*np.exp(-1*E1/R/T)*LIGH-2.000000e+00*A62*T**n62*np.exp(-1*E62/R/T)*RLIGM2A*RLIGM2A-A15*T**n15*np.exp(-1*E15/R/T)*RLIGM2A-A266*T**n266*np.exp(-1*E266/R/T)*KETM2*RLIGM2A-A245*T**n245*np.exp(-1*E245/R/T)*RLIGM2A*ADIOM2+\nA176*T**n176*np.exp(-1*E176/R/T)*PRFET3*LIGM2+A195*T**n195*np.exp(-1*E195/R/T)*RADIO*LIGM2+A178*T**n178*np.exp(-1*E178/R/T)*LIGM2*RADIOM2-A140*T**n140*np.exp(-1*E140/R/T)*RLIGM2A*PLIGH+A193*T**n193*np.exp(-1*E193/R/T)*RCH3O*LIGM2-A161*T**n161*np.exp(-1*E161/R/T)*PLIGM2*RLIGM2A+\nA187*T**n187*np.exp(-1*E187/R/T)*RC3H3O*LIGM2+A192*T**n192*np.exp(-1*E192/R/T)*LIGM2*OH+A186*T**n186*np.exp(-1*E186/R/T)*PRADIO*LIGM2-A224*T**n224*np.exp(-1*E224/R/T)*RLIGM2A*PFET3M2+A184*T**n184*np.exp(-1*E184/R/T)*PRKETM2*LIGM2+A12*T**n12*np.exp(-1*E12/R/T)*RLIGH+\nA188*T**n188*np.exp(-1*E188/R/T)*RLIGB*LIGM2+A175*T**n175*np.exp(-1*E175/R/T)*RC3H5O2*LIGM2+A183*T**n183*np.exp(-1*E183/R/T)*RCH3*LIGM2-A203*T**n203*np.exp(-1*E203/R/T)*RLIGM2A*LIGM2-A329*T**n329*np.exp(-1*E329/R/T)*LIG*RLIGM2A+A181*T**n181*np.exp(-1*E181/R/T)*RLIGM2B*LIGM2+\nA179*T**n179*np.exp(-1*E179/R/T)*PRFET3M2*LIGM2+A177*T**n177*np.exp(-1*E177/R/T)*RC3H7O2*LIGM2+A190*T**n190*np.exp(-1*E190/R/T)*PRADIOM2*LIGM2;\n/*ODEforLIGM2=>LIGM2=8*/\ndydt[LIGM2]=-A191*T**n191*np.exp(-1*E191/R/T)*RMGUAI*LIGM2+A308*T**n308*np.exp(-1*E308/R/T)*LIG*RLIGM2A-A210*T**n210*np.exp(-1*E210/R/T)*RLIGA*LIGM2-A194*T**n194*np.exp(-1*E194/R/T)*LIGM2*RPHENOL+A392*T**n392*np.exp(-1*E392/R/T)*RLIGM2A*KET-\nA212*T**n212*np.exp(-1*E212/R/T)*RMGUAI*LIGM2+A139*T**n139*np.exp(-1*E139/R/T)*RLIGM2B*PLIGH-A199*T**n199*np.exp(-1*E199/R/T)*LIGM2*RADIOM2+A350*T**n350*np.exp(-1*E350/R/T)*RLIGM2A*PFET3+A307*T**n307*np.exp(-1*E307/R/T)*RLIGM2B*LIG-A2*T**n2*np.exp(-1*E2/R/T)*LIGM2-\nA196*T**n196*np.exp(-1*E196/R/T)*RC3H5O2*LIGM2+A107*T**n107*np.exp(-1*E107/R/T)*LIGO+A371*T**n371*np.exp(-1*E371/R/T)*RLIGM2A*ADIO+A328*T**n328*np.exp(-1*E328/R/T)*RLIGM2B*LIG+A265*T**n265*np.exp(-1*E265/R/T)*RLIGM2B*KETM2+A161*T**n161*np.exp(-1*E161/R/T)*PLIGM2*RLIGM2A-\nA192*T**n192*np.exp(-1*E192/R/T)*LIGM2*OH+A224*T**n224*np.exp(-1*E224/R/T)*RLIGM2A*PFET3M2-A213*T**n213*np.exp(-1*E213/R/T)*LIGM2*OH-A208*T**n208*np.exp(-1*E208/R/T)*RC3H3O*LIGM2-A175*T**n175*np.exp(-1*E175/R/T)*RC3H5O2*LIGM2-A183*T**n183*np.exp(-1*E183/R/T)*RCH3*LIGM2-\nA200*T**n200*np.exp(-1*E200/R/T)*PRFET3M2*LIGM2+A329*T**n329*np.exp(-1*E329/R/T)*LIG*RLIGM2A-A179*T**n179*np.exp(-1*E179/R/T)*PRFET3M2*LIGM2-A209*T**n209*np.exp(-1*E209/R/T)*RLIGB*LIGM2-A207*T**n207*np.exp(-1*E207/R/T)*PRADIO*LIGM2-A205*T**n205*np.exp(-1*E205/R/T)*PRKETM2*LIGM2-\nA198*T**n198*np.exp(-1*E198/R/T)*RC3H7O2*LIGM2-A177*T**n177*np.exp(-1*E177/R/T)*RC3H7O2*LIGM2-A214*T**n214*np.exp(-1*E214/R/T)*RCH3O*LIGM2-A190*T**n190*np.exp(-1*E190/R/T)*PRADIOM2*LIGM2-A216*T**n216*np.exp(-1*E216/R/T)*RADIO*LIGM2-A201*T**n201*np.exp(-1*E201/R/T)*PRLIGH*LIGM2+\nA287*T**n287*np.exp(-1*E287/R/T)*RLIGM2A*C10H2-A211*T**n211*np.exp(-1*E211/R/T)*PRADIOM2*LIGM2-A189*T**n189*np.exp(-1*E189/R/T)*RLIGA*LIGM2+A286*T**n286*np.exp(-1*E286/R/T)*RLIGM2B*C10H2-A204*T**n204*np.exp(-1*E204/R/T)*RCH3*LIGM2-A180*T**n180*np.exp(-1*E180/R/T)*PRLIGH*LIGM2-\nA215*T**n215*np.exp(-1*E215/R/T)*LIGM2*RPHENOL-A185*T**n185*np.exp(-1*E185/R/T)*RKET*LIGM2+A119*T**n119*np.exp(-1*E119/R/T)*RLIGM2A*LIGH+A244*T**n244*np.exp(-1*E244/R/T)*RLIGM2B*ADIOM2-A206*T**n206*np.exp(-1*E206/R/T)*RKET*LIGM2+A266*T**n266*np.exp(-1*E266/R/T)*KETM2*RLIGM2A+\nA223*T**n223*np.exp(-1*E223/R/T)*RLIGM2B*PFET3M2+A245*T**n245*np.exp(-1*E245/R/T)*RLIGM2A*ADIOM2-A195*T**n195*np.exp(-1*E195/R/T)*RADIO*LIGM2-A176*T**n176*np.exp(-1*E176/R/T)*PRFET3*LIGM2-A178*T**n178*np.exp(-1*E178/R/T)*LIGM2*RADIOM2-A193*T**n193*np.exp(-1*E193/R/T)*RCH3O*LIGM2+\nA140*T**n140*np.exp(-1*E140/R/T)*RLIGM2A*PLIGH-A187*T**n187*np.exp(-1*E187/R/T)*RC3H3O*LIGM2-A186*T**n186*np.exp(-1*E186/R/T)*PRADIO*LIGM2-A184*T**n184*np.exp(-1*E184/R/T)*PRKETM2*LIGM2-A188*T**n188*np.exp(-1*E188/R/T)*RLIGB*LIGM2+A370*T**n370*np.exp(-1*E370/R/T)*RLIGM2B*ADIO-\nA197*T**n197*np.exp(-1*E197/R/T)*PRFET3*LIGM2+A118*T**n118*np.exp(-1*E118/R/T)*RLIGM2B*LIGH+A349*T**n349*np.exp(-1*E349/R/T)*RLIGM2B*PFET3+A391*T**n391*np.exp(-1*E391/R/T)*RLIGM2B*KET+A160*T**n160*np.exp(-1*E160/R/T)*RLIGM2B*PLIGM2;\n/*ODEforRPHENOXM2=>RPHENOXM2=9*/\ndydt[RPHENOXM2]=-A10*T**n10*np.exp(-1*E10/R/T)*RPHENOXM2-A49*T**n49*np.exp(-1*E49/R/T)*C10H2M2*RPHENOXM2+A16*T**n16*np.exp(-1*E16/R/T)*PRLIGM2A-A48*T**n48*np.exp(-1*E48/R/T)*C10H2M4*RPHENOXM2-A45*T**n45*np.exp(-1*E45/R/T)*KET*RPHENOXM2-\nA32*T**n32*np.exp(-1*E32/R/T)*ADIOM2*RPHENOXM2-A83*T**n83*np.exp(-1*E83/R/T)*CHAR*RPHENOXM2-A46*T**n46*np.exp(-1*E46/R/T)*KETD*RPHENOXM2-A44*T**n44*np.exp(-1*E44/R/T)*ADIO*RPHENOXM2-A33*T**n33*np.exp(-1*E33/R/T)*KETM2*RPHENOXM2+A3*T**n3*np.exp(-1*E3/R/T)*PLIGM2-A35*T**n35*np.exp(-1*E35/R/T)*SYNAPYL*RPHENOXM2\n-A34*T**n34*np.exp(-1*E34/R/T)*RPHENOXM2*KETDM2+A2*T**n2*np.exp(-1*E2/R/T)*LIGM2-A47*T**n47*np.exp(-1*E47/R/T)*COUMARYL*RPHENOXM2-A54*T**n54*np.exp(-1*E54/R/T)*RCH3*RPHENOXM2-A53*T**n53*np.exp(-1*E53/R/T)*RCH3O*RPHENOXM2+A15*T**n15*np.exp(-1*E15/R/T)*RLIGM2A+\nA18*T**n18*np.exp(-1*E18/R/T)*PRFET3M2;\n/*ODEforRADIOM2=>RADIOM2=10*/\ndydt[RADIOM2]=-A367*T**n367*np.exp(-1*E367/R/T)*ADIO*RADIOM2-2.000000e+00*A60*T**n60*np.exp(-1*E60/R/T)*RADIOM2*RADIOM2-A283*T**n283*np.exp(-1*E283/R/T)*C10H2*RADIOM2+\nA256*T**n256*np.exp(-1*E256/R/T)*RCH3O*ADIOM2-A115*T**n115*np.exp(-1*E115/R/T)*LIGH*RADIOM2+A249*T**n249*np.exp(-1*E249/R/T)*PRADIO*ADIOM2-A199*T**n199*np.exp(-1*E199/R/T)*LIGM2*RADIOM2+A243*T**n243*np.exp(-1*E243/R/T)*PRLIGH*ADIOM2-A304*T**n304*np.exp(-1*E304/R/T)*LIG*RADIOM2-\nA23*T**n23*np.exp(-1*E23/R/T)*RADIOM2+A244*T**n244*np.exp(-1*E244/R/T)*RLIGM2B*ADIOM2+A2*T**n2*np.exp(-1*E2/R/T)*LIGM2+A257*T**n257*np.exp(-1*E257/R/T)*ADIOM2*RPHENOL-A262*T**n262*np.exp(-1*E262/R/T)*KETM2*RADIOM2-A220*T**n220*np.exp(-1*E220/R/T)*PFET3M2*RADIOM2+\nA242*T**n242*np.exp(-1*E242/R/T)*PRFET3M2*ADIOM2+A238*T**n238*np.exp(-1*E238/R/T)*RC3H5O2*ADIOM2+A255*T**n255*np.exp(-1*E255/R/T)*OH*ADIOM2+A247*T**n247*np.exp(-1*E247/R/T)*PRKETM2*ADIOM2-A346*T**n346*np.exp(-1*E346/R/T)*PFET3*RADIOM2-A14*T**n14*np.exp(-1*E14/R/T)*RADIOM2+\nA245*T**n245*np.exp(-1*E245/R/T)*RLIGM2A*ADIOM2-A178*T**n178*np.exp(-1*E178/R/T)*LIGM2*RADIOM2+A246*T**n246*np.exp(-1*E246/R/T)*RCH3*ADIOM2-A136*T**n136*np.exp(-1*E136/R/T)*PLIGH*RADIOM2+A253*T**n253*np.exp(-1*E253/R/T)*PRADIOM2*ADIOM2-A388*T**n388*np.exp(-1*E388/R/T)*KET*RADIOM2-\nA325*T**n325*np.exp(-1*E325/R/T)*LIG*RADIOM2+A240*T**n240*np.exp(-1*E240/R/T)*RC3H7O2*ADIOM2-A157*T**n157*np.exp(-1*E157/R/T)*PLIGM2*RADIOM2+A239*T**n239*np.exp(-1*E239/R/T)*PRFET3*ADIOM2+A252*T**n252*np.exp(-1*E252/R/T)*RLIGA*ADIOM2+A254*T**n254*np.exp(-1*E254/R/T)*RMGUAI*ADIOM2+\nA250*T**n250*np.exp(-1*E250/R/T)*RC3H3O*ADIOM2+A248*T**n248*np.exp(-1*E248/R/T)*RKET*ADIOM2+A258*T**n258*np.exp(-1*E258/R/T)*RADIO*ADIOM2+A251*T**n251*np.exp(-1*E251/R/T)*RLIGB*ADIOM2;\n/*ODEforPLIGM2=>PLIGM2=11*/\ndydt[PLIGM2]=-A167*T**n167*np.exp(-1*E167/R/T)*PLIGM2*RLIGB-A173*T**n173*np.exp(-1*E173/R/T)*PLIGM2*RPHENOL-A161*T**n161*np.exp(-1*E161/R/T)*PLIGM2*RLIGM2A-A168*T**n168*np.exp(-1*E168/R/T)*PLIGM2*RLIGA-A172*T**n172*np.exp(-1*E172/R/T)*PLIGM2*RCH3O-\nA169*T**n169*np.exp(-1*E169/R/T)*PLIGM2*PRADIOM2-A158*T**n158*np.exp(-1*E158/R/T)*PLIGM2*PRFET3M2-A171*T**n171*np.exp(-1*E171/R/T)*PLIGM2*OH-A157*T**n157*np.exp(-1*E157/R/T)*PLIGM2*RADIOM2-A170*T**n170*np.exp(-1*E170/R/T)*PLIGM2*RMGUAI-A3*T**n3*np.exp(-1*E3/R/T)*PLIGM2+\nA31*T**n31*np.exp(-1*E31/R/T)*PLIGO-A155*T**n155*np.exp(-1*E155/R/T)*PLIGM2*PRFET3-A156*T**n156*np.exp(-1*E156/R/T)*PLIGM2*RC3H7O2-A165*T**n165*np.exp(-1*E165/R/T)*PLIGM2*PRADIO-A174*T**n174*np.exp(-1*E174/R/T)*RADIO*PLIGM2-A164*T**n164*np.exp(-1*E164/R/T)*PLIGM2*RKET-\nA154*T**n154*np.exp(-1*E154/R/T)*PLIGM2*RC3H5O2-A162*T**n162*np.exp(-1*E162/R/T)*PLIGM2*RCH3-A160*T**n160*np.exp(-1*E160/R/T)*RLIGM2B*PLIGM2-A159*T**n159*np.exp(-1*E159/R/T)*PLIGM2*PRLIGH-A166*T**n166*np.exp(-1*E166/R/T)*PLIGM2*RC3H3O-A163*T**n163*np.exp(-1*E163/R/T)*PLIGM2*PRKETM2;\n/*ODEforPRADIOM2=>PRADIOM2=12*/\ndydt[PRADIOM2]=-A148*T**n148*np.exp(-1*E148/R/T)*PRADIOM2*PLIGH-A379*T**n379*np.exp(-1*E379/R/T)*PRADIOM2*ADIO-A316*T**n316*np.exp(-1*E316/R/T)*LIG*PRADIOM2-A253*T**n253*np.exp(-1*E253/R/T)*PRADIOM2*ADIOM2-A211*T**n211*np.exp(-1*E211/R/T)*PRADIOM2*LIGM2\n-A127*T**n127*np.exp(-1*E127/R/T)*PRADIOM2*LIGH-A295*T**n295*np.exp(-1*E295/R/T)*PRADIOM2*C10H2-A232*T**n232*np.exp(-1*E232/R/T)*PRADIOM2*PFET3M2-A71*T**n71*np.exp(-1*E71/R/T)*PRADIOM2*RCH3-A169*T**n169*np.exp(-1*E169/R/T)*PLIGM2*PRADIOM2+A3*T**n3*np.exp(-1*E3/R/T)*PLIGM2-\nA274*T**n274*np.exp(-1*E274/R/T)*KETM2*PRADIOM2+A6*T**n6*np.exp(-1*E6/R/T)*PADIOM2-A58*T**n58*np.exp(-1*E58/R/T)*PRADIOM2*RCH3O-2.000000e+00*A110*T**n110*np.exp(-1*E110/R/T)*PRADIOM2*PRADIOM2-A400*T**n400*np.exp(-1*E400/R/T)*PRADIOM2*KET-\nA337*T**n337*np.exp(-1*E337/R/T)*LIG*PRADIOM2-A358*T**n358*np.exp(-1*E358/R/T)*PRADIOM2*PFET3-A190*T**n190*np.exp(-1*E190/R/T)*PRADIOM2*LIGM2;\n/*ODEforLIG=>LIG=13*/\ndydt[LIG]=+A230*T**n230*np.exp(-1*E230/R/T)*RLIGB*PFET3M2+A356*T**n356*np.exp(-1*E356/R/T)*PFET3*RLIGB-A339*T**n339*np.exp(-1*E339/R/T)*LIG*OH-A308*T**n308*np.exp(-1*E308/R/T)*LIG*RLIGM2A-A316*T**n316*np.exp(-1*E316/R/T)*LIG*PRADIOM2-\nA306*T**n306*np.exp(-1*E306/R/T)*LIG*PRLIGH-A327*T**n327*np.exp(-1*E327/R/T)*LIG*PRLIGH+A210*T**n210*np.exp(-1*E210/R/T)*RLIGA*LIGM2-A318*T**n318*np.exp(-1*E318/R/T)*LIG*OH-A310*T**n310*np.exp(-1*E310/R/T)*LIG*PRKETM2-A319*T**n319*np.exp(-1*E319/R/T)*LIG*RCH3O-\nA304*T**n304*np.exp(-1*E304/R/T)*LIG*RADIOM2-A307*T**n307*np.exp(-1*E307/R/T)*RLIGM2B*LIG-A317*T**n317*np.exp(-1*E317/R/T)*LIG*RMGUAI-A305*T**n305*np.exp(-1*E305/R/T)*LIG*PRFET3M2+A273*T**n273*np.exp(-1*E273/R/T)*KETM2*RLIGA-A312*T**n312*np.exp(-1*E312/R/T)*LIG*PRADIO+\nA147*T**n147*np.exp(-1*E147/R/T)*PLIGH*RLIGA-A342*T**n342*np.exp(-1*E342/R/T)*RADIO*LIG-A328*T**n328*np.exp(-1*E328/R/T)*RLIGM2B*LIG-A321*T**n321*np.exp(-1*E321/R/T)*RADIO*LIG-A330*T**n330*np.exp(-1*E330/R/T)*LIG*RCH3-A331*T**n331*np.exp(-1*E331/R/T)*LIG*PRKETM2-\nA332*T**n332*np.exp(-1*E332/R/T)*LIG*RKET-A324*T**n324*np.exp(-1*E324/R/T)*LIG*RC3H7O2+A252*T**n252*np.exp(-1*E252/R/T)*RLIGA*ADIOM2+A108*T**n108*np.exp(-1*E108/R/T)*LIGC+A126*T**n126*np.exp(-1*E126/R/T)*LIGH*RLIGA+A125*T**n125*np.exp(-1*E125/R/T)*RLIGB*LIGH-\nA329*T**n329*np.exp(-1*E329/R/T)*LIG*RLIGM2A-A341*T**n341*np.exp(-1*E341/R/T)*LIG*RPHENOL+A209*T**n209*np.exp(-1*E209/R/T)*RLIGB*LIGM2-A313*T**n313*np.exp(-1*E313/R/T)*LIG*RC3H3O+A272*T**n272*np.exp(-1*E272/R/T)*KETM2*RLIGB+A231*T**n231*np.exp(-1*E231/R/T)*RLIGA*PFET3M2+\nA293*T**n293*np.exp(-1*E293/R/T)*RLIGB*C10H2-A337*T**n337*np.exp(-1*E337/R/T)*LIG*PRADIOM2+A251*T**n251*np.exp(-1*E251/R/T)*RLIGB*ADIOM2+A189*T**n189*np.exp(-1*E189/R/T)*RLIGA*LIGM2+A377*T**n377*np.exp(-1*E377/R/T)*RLIGB*ADIO+A357*T**n357*np.exp(-1*E357/R/T)*PFET3*RLIGA+\nA398*T**n398*np.exp(-1*E398/R/T)*RLIGB*KET-A326*T**n326*np.exp(-1*E326/R/T)*LIG*PRFET3M2-A311*T**n311*np.exp(-1*E311/R/T)*LIG*RKET-A334*T**n334*np.exp(-1*E334/R/T)*LIG*RC3H3O-A303*T**n303*np.exp(-1*E303/R/T)*LIG*RC3H7O2-A309*T**n309*np.exp(-1*E309/R/T)*LIG*RCH3-\nA301*T**n301*np.exp(-1*E301/R/T)*LIG*RC3H5O2+A146*T**n146*np.exp(-1*E146/R/T)*PLIGH*RLIGB-A323*T**n323*np.exp(-1*E323/R/T)*LIG*PRFET3+A167*T**n167*np.exp(-1*E167/R/T)*PLIGM2*RLIGB-A4*T**n4*np.exp(-1*E4/R/T)*LIG+A168*T**n168*np.exp(-1*E168/R/T)*PLIGM2*RLIGA-\nA333*T**n333*np.exp(-1*E333/R/T)*LIG*PRADIO-A320*T**n320*np.exp(-1*E320/R/T)*LIG*RPHENOL+A378*T**n378*np.exp(-1*E378/R/T)*ADIO*RLIGA+A294*T**n294*np.exp(-1*E294/R/T)*RLIGA*C10H2-A325*T**n325*np.exp(-1*E325/R/T)*LIG*RADIOM2+A188*T**n188*np.exp(-1*E188/R/T)*RLIGB*LIGM2-\nA338*T**n338*np.exp(-1*E338/R/T)*LIG*RMGUAI-A322*T**n322*np.exp(-1*E322/R/T)*LIG*RC3H5O2-A302*T**n302*np.exp(-1*E302/R/T)*LIG*PRFET3-A340*T**n340*np.exp(-1*E340/R/T)*LIG*RCH3O+A399*T**n399*np.exp(-1*E399/R/T)*KET*RLIGA;\n/*ODEforRPHENOX=>RPHENOX=14*/\ndydt[RPHENOX]=-A39*T**n39*np.exp(-1*E39/R/T)*SYNAPYL*RPHENOX+A4*T**n4*np.exp(-1*E4/R/T)*LIG-A42*T**n42*np.exp(-1*E42/R/T)*KETD*RPHENOX-A52*T**n52*np.exp(-1*E52/R/T)*RCH3O*RPHENOX-A81*T**n81*np.exp(-1*E81/R/T)*RC3H3O*RPHENOX-\nA36*T**n36*np.exp(-1*E36/R/T)*ADIOM2*RPHENOX-A41*T**n41*np.exp(-1*E41/R/T)*KET*RPHENOX-A50*T**n50*np.exp(-1*E50/R/T)*C10H2M4*RPHENOX-A40*T**n40*np.exp(-1*E40/R/T)*ADIO*RPHENOX-A80*T**n80*np.exp(-1*E80/R/T)*RPHENOX*RPHENOL-A37*T**n37*np.exp(-1*E37/R/T)*KETM2*RPHENOX+\nA22*T**n22*np.exp(-1*E22/R/T)*PRFET3+A28*T**n28*np.exp(-1*E28/R/T)*C10H2M4+A5*T**n5*np.exp(-1*E5/R/T)*PLIG-A43*T**n43*np.exp(-1*E43/R/T)*COUMARYL*RPHENOX-A11*T**n11*np.exp(-1*E11/R/T)*RPHENOX-A51*T**n51*np.exp(-1*E51/R/T)*C10H2M2*RPHENOX-A38*T**n38*np.exp(-1*E38/R/T)*RPHENOX*KETDM2\n-A82*T**n82*np.exp(-1*E82/R/T)*CHAR*RPHENOX-A73*T**n73*np.exp(-1*E73/R/T)*RLIGB*RPHENOX-A55*T**n55*np.exp(-1*E55/R/T)*RCH3*RPHENOX+A29*T**n29*np.exp(-1*E29/R/T)*C10H2M2+A20*T**n20*np.exp(-1*E20/R/T)*RLIGA;\n/*ODEforRADIO=>RADIO=15*/\ndydt[RADIO]=+A367*T**n367*np.exp(-1*E367/R/T)*ADIO*RADIOM2-2.000000e+00*A74*T**n74*np.exp(-1*E74/R/T)*RADIO*RADIO-A216*T**n216*np.exp(-1*E216/R/T)*RADIO*LIGM2-A153*T**n153*np.exp(-1*E153/R/T)*RADIO*PLIGH\n-A132*T**n132*np.exp(-1*E132/R/T)*RADIO*LIGH-A279*T**n279*np.exp(-1*E279/R/T)*RADIO*KETM2+A372*T**n372*np.exp(-1*E372/R/T)*RCH3*ADIO+A364*T**n364*np.exp(-1*E364/R/T)*ADIO*RC3H5O2+A377*T**n377*np.exp(-1*E377/R/T)*RLIGB*ADIO+A373*T**n373*np.exp(-1*E373/R/T)*PRKETM2*ADIO-\nA405*T**n405*np.exp(-1*E405/R/T)*RADIO*KET+A381*T**n381*np.exp(-1*E381/R/T)*ADIO*OH-A174*T**n174*np.exp(-1*E174/R/T)*RADIO*PLIGM2+A371*T**n371*np.exp(-1*E371/R/T)*RLIGM2A*ADIO+A380*T**n380*np.exp(-1*E380/R/T)*RMGUAI*ADIO-A321*T**n321*np.exp(-1*E321/R/T)*RADIO*LIG-\nA342*T**n342*np.exp(-1*E342/R/T)*RADIO*LIG-A19*T**n19*np.exp(-1*E19/R/T)*RADIO+A379*T**n379*np.exp(-1*E379/R/T)*PRADIOM2*ADIO-A195*T**n195*np.exp(-1*E195/R/T)*RADIO*LIGM2+A365*T**n365*np.exp(-1*E365/R/T)*ADIO*PRFET3+A366*T**n366*np.exp(-1*E366/R/T)*RC3H7O2*ADIO+\nA4*T**n4*np.exp(-1*E4/R/T)*LIG+A383*T**n383*np.exp(-1*E383/R/T)*ADIO*RPHENOL+A382*T**n382*np.exp(-1*E382/R/T)*RCH3O*ADIO+A378*T**n378*np.exp(-1*E378/R/T)*ADIO*RLIGA+A370*T**n370*np.exp(-1*E370/R/T)*RLIGM2B*ADIO-A25*T**n25*np.exp(-1*E25/R/T)*RADIO-\nA237*T**n237*np.exp(-1*E237/R/T)*RADIO*PFET3M2+A375*T**n375*np.exp(-1*E375/R/T)*ADIO*PRADIO+A368*T**n368*np.exp(-1*E368/R/T)*PRFET3M2*ADIO+A376*T**n376*np.exp(-1*E376/R/T)*RC3H3O*ADIO-A258*T**n258*np.exp(-1*E258/R/T)*RADIO*ADIOM2-A363*T**n363*np.exp(-1*E363/R/T)*RADIO*PFET3-\nA300*T**n300*np.exp(-1*E300/R/T)*RADIO*C10H2+A369*T**n369*np.exp(-1*E369/R/T)*PRLIGH*ADIO+A374*T**n374*np.exp(-1*E374/R/T)*ADIO*RKET;\n/*ODEforPLIG=>PLIG=16*/\ndydt[PLIG]=+A30*T**n30*np.exp(-1*E30/R/T)*PLIGC-A5*T**n5*np.exp(-1*E5/R/T)*PLIG;\n/*ODEforPRADIO=>PRADIO=17*/\ndydt[PRADIO]=-A123*T**n123*np.exp(-1*E123/R/T)*LIGH*PRADIO-A270*T**n270*np.exp(-1*E270/R/T)*KETM2*PRADIO-A333*T**n333*np.exp(-1*E333/R/T)*LIG*PRADIO-A186*T**n186*np.exp(-1*E186/R/T)*PRADIO*LIGM2-A396*T**n396*np.exp(-1*E396/R/T)*KET*PRADIO-\nA249*T**n249*np.exp(-1*E249/R/T)*PRADIO*ADIOM2+A5*T**n5*np.exp(-1*E5/R/T)*PLIG-2.000000e+00*A109*T**n109*np.exp(-1*E109/R/T)*PRADIO*PRADIO-A165*T**n165*np.exp(-1*E165/R/T)*PLIGM2*PRADIO-A354*T**n354*np.exp(-1*E354/R/T)*PFET3*PRADIO-\nA207*T**n207*np.exp(-1*E207/R/T)*PRADIO*LIGM2-A312*T**n312*np.exp(-1*E312/R/T)*LIG*PRADIO-A375*T**n375*np.exp(-1*E375/R/T)*ADIO*PRADIO+A7*T**n7*np.exp(-1*E7/R/T)*PADIO-A228*T**n228*np.exp(-1*E228/R/T)*PRADIO*PFET3M2-A291*T**n291*np.exp(-1*E291/R/T)*C10H2*PRADIO-\nA144*T**n144*np.exp(-1*E144/R/T)*PLIGH*PRADIO;\n/*ODEforPADIOM2=>PADIOM2=18*/\ndydt[PADIOM2]=+A148*T**n148*np.exp(-1*E148/R/T)*PRADIOM2*PLIGH+A379*T**n379*np.exp(-1*E379/R/T)*PRADIOM2*ADIO+A316*T**n316*np.exp(-1*E316/R/T)*LIG*PRADIOM2+A253*T**n253*np.exp(-1*E253/R/T)*PRADIOM2*ADIOM2+A211*T**n211*np.exp(-1*E211/R/T)*PRADIOM2*LIGM2\n+A127*T**n127*np.exp(-1*E127/R/T)*PRADIOM2*LIGH+A295*T**n295*np.exp(-1*E295/R/T)*PRADIOM2*C10H2+A232*T**n232*np.exp(-1*E232/R/T)*PRADIOM2*PFET3M2+A71*T**n71*np.exp(-1*E71/R/T)*PRADIOM2*RCH3+A169*T**n169*np.exp(-1*E169/R/T)*PLIGM2*PRADIOM2+A274*T**n274*np.exp(-1*E274/R/T)*KETM2*PRADIOM2-\nA6*T**n6*np.exp(-1*E6/R/T)*PADIOM2+A58*T**n58*np.exp(-1*E58/R/T)*PRADIOM2*RCH3O+A400*T**n400*np.exp(-1*E400/R/T)*PRADIOM2*KET+A337*T**n337*np.exp(-1*E337/R/T)*LIG*PRADIOM2+A358*T**n358*np.exp(-1*E358/R/T)*PRADIOM2*PFET3+A190*T**n190*np.exp(-1*E190/R/T)*PRADIOM2*LIGM2;\n/*ODEforPADIO=>PADIO=19*/\ndydt[PADIO]=+A123*T**n123*np.exp(-1*E123/R/T)*LIGH*PRADIO+A165*T**n165*np.exp(-1*E165/R/T)*PLIGM2*PRADIO+A354*T**n354*np.exp(-1*E354/R/T)*PFET3*PRADIO+A375*T**n375*np.exp(-1*E375/R/T)*ADIO*PRADIO+A312*T**n312*np.exp(-1*E312/R/T)*LIG*PRADIO+\nA270*T**n270*np.exp(-1*E270/R/T)*KETM2*PRADIO+A207*T**n207*np.exp(-1*E207/R/T)*PRADIO*LIGM2-A7*T**n7*np.exp(-1*E7/R/T)*PADIO+A333*T**n333*np.exp(-1*E333/R/T)*LIG*PRADIO+A228*T**n228*np.exp(-1*E228/R/T)*PRADIO*PFET3M2+A186*T**n186*np.exp(-1*E186/R/T)*PRADIO*LIGM2+\nA396*T**n396*np.exp(-1*E396/R/T)*KET*PRADIO+A291*T**n291*np.exp(-1*E291/R/T)*C10H2*PRADIO+A249*T**n249*np.exp(-1*E249/R/T)*PRADIO*ADIOM2+A144*T**n144*np.exp(-1*E144/R/T)*PLIGH*PRADIO;\n/*ODEforPKETM2=>PKETM2=20*/\ndydt[PKETM2]=+A16*T**n16*np.exp(-1*E16/R/T)*PRLIGM2A+A394*T**n394*np.exp(-1*E394/R/T)*PRKETM2*KET+A142*T**n142*np.exp(-1*E142/R/T)*PRKETM2*PLIGH+A289*T**n289*np.exp(-1*E289/R/T)*PRKETM2*C10H2+A310*T**n310*np.exp(-1*E310/R/T)*LIG*PRKETM2+\nA184*T**n184*np.exp(-1*E184/R/T)*PRKETM2*LIGM2+A373*T**n373*np.exp(-1*E373/R/T)*PRKETM2*ADIO+A331*T**n331*np.exp(-1*E331/R/T)*LIG*PRKETM2+A59*T**n59*np.exp(-1*E59/R/T)*PRKETM2*RCH3O+A268*T**n268*np.exp(-1*E268/R/T)*KETM2*PRKETM2+A352*T**n352*np.exp(-1*E352/R/T)*PRKETM2*PFET3+\nA121*T**n121*np.exp(-1*E121/R/T)*PRKETM2*LIGH+A72*T**n72*np.exp(-1*E72/R/T)*PRKETM2*RCH3+A226*T**n226*np.exp(-1*E226/R/T)*PRKETM2*PFET3M2-A8*T**n8*np.exp(-1*E8/R/T)*PKETM2+A205*T**n205*np.exp(-1*E205/R/T)*PRKETM2*LIGM2+A163*T**n163*np.exp(-1*E163/R/T)*PLIGM2*PRKETM2+\nA247*T**n247*np.exp(-1*E247/R/T)*PRKETM2*ADIOM2;\n/*ODEforPRKETM2=>PRKETM2=21*/\ndydt[PRKETM2]=-A394*T**n394*np.exp(-1*E394/R/T)*PRKETM2*KET-A142*T**n142*np.exp(-1*E142/R/T)*PRKETM2*PLIGH-A289*T**n289*np.exp(-1*E289/R/T)*PRKETM2*C10H2-A310*T**n310*np.exp(-1*E310/R/T)*LIG*PRKETM2-A184*T**n184*np.exp(-1*E184/R/T)*PRKETM2*LIGM2\n-A373*T**n373*np.exp(-1*E373/R/T)*PRKETM2*ADIO-A331*T**n331*np.exp(-1*E331/R/T)*LIG*PRKETM2-A9*T**n9*np.exp(-1*E9/R/T)*PRKETM2-A59*T**n59*np.exp(-1*E59/R/T)*PRKETM2*RCH3O-A268*T**n268*np.exp(-1*E268/R/T)*KETM2*PRKETM2-A352*T**n352*np.exp(-1*E352/R/T)*PRKETM2*PFET3-\nA121*T**n121*np.exp(-1*E121/R/T)*PRKETM2*LIGH-A72*T**n72*np.exp(-1*E72/R/T)*PRKETM2*RCH3-A226*T**n226*np.exp(-1*E226/R/T)*PRKETM2*PFET3M2+A8*T**n8*np.exp(-1*E8/R/T)*PKETM2-A205*T**n205*np.exp(-1*E205/R/T)*PRKETM2*LIGM2-A163*T**n163*np.exp(-1*E163/R/T)*PLIGM2*PRKETM2-\nA247*T**n247*np.exp(-1*E247/R/T)*PRKETM2*ADIOM2;\n/*ODEforKETDM2=>KETDM2=22*/\ndydt[KETDM2]=-A34*T**n34*np.exp(-1*E34/R/T)*RPHENOXM2*KETDM2-A86*T**n86*np.exp(-1*E86/R/T)*KETDM2-A38*T**n38*np.exp(-1*E38/R/T)*RPHENOX*KETDM2+A24*T**n24*np.exp(-1*E24/R/T)*RKETM2+A9*T**n9*np.exp(-1*E9/R/T)*PRKETM2;\n/*ODEforCO=>CO=23*/\ndydt[CO]=+A11*T**n11*np.exp(-1*E11/R/T)*RPHENOX+A10*T**n10*np.exp(-1*E10/R/T)*RPHENOXM2+A95*T**n95*np.exp(-1*E95/R/T)*PCOH+A111*T**n111*np.exp(-1*E111/R/T)*PCHO+A94*T**n94*np.exp(-1*E94/R/T)*PCOS;\n/*ODEforC10H2M4=>C10H2M4=24*/\ndydt[C10H2M4]=+A60*T**n60*np.exp(-1*E60/R/T)*RADIOM2*RADIOM2+5.000000e-01*A10*T**n10*np.exp(-1*E10/R/T)*RPHENOXM2+5.000000e-01*\nA49*T**n49*np.exp(-1*E49/R/T)*C10H2M2*RPHENOXM2+5.000000e-01*A39*T**n39*np.exp(-1*E39/R/T)*SYNAPYL*RPHENOX+A63*T**n63*np.exp(-1*E63/R/T)*RMGUAI*RMGUAI+5.000000e-01*A45*T**n45*np.exp(-1*E45/R/T)*KET*RPHENOXM2+\nA32*T**n32*np.exp(-1*E32/R/T)*ADIOM2*RPHENOXM2+5.000000e-01*A83*T**n83*np.exp(-1*E83/R/T)*CHAR*RPHENOXM2+5.000000e-01*A36*T**n36*np.exp(-1*E36/R/T)*ADIOM2*RPHENOX-\n5.000000e-01*A50*T**n50*np.exp(-1*E50/R/T)*C10H2M4*RPHENOX+5.000000e-01*A37*T**n37*np.exp(-1*E37/R/T)*KETM2*RPHENOX+5.000000e-01*\nA46*T**n46*np.exp(-1*E46/R/T)*KETD*RPHENOXM2-5.000000e-01*A28*T**n28*np.exp(-1*E28/R/T)*C10H2M4+5.000000e-01*A44*T**n44*np.exp(-1*E44/R/T)*ADIO*RPHENOXM2+A33*T**n33*np.exp(-1*E33/R/T)*KETM2*RPHENOXM2+\nA35*T**n35*np.exp(-1*E35/R/T)*SYNAPYL*RPHENOXM2+A65*T**n65*np.exp(-1*E65/R/T)*PRFET3M2*PRFET3M2+A34*T**n34*np.exp(-1*E34/R/T)*RPHENOXM2*KETDM2+A64*T**n64*np.exp(-1*E64/R/T)*RKETM2*RKETM2+2.000000e+00*A61*T**n61*np.exp(-1*E61/R/T)*RLIGM2B*RLIGM2B+\n5.000000e-01*A47*T**n47*np.exp(-1*E47/R/T)*COUMARYL*RPHENOXM2+2.000000e+00*A79*T**n79*np.exp(-1*E79/R/T)*RLIGH*RLIGH+5.000000e-01*\nA54*T**n54*np.exp(-1*E54/R/T)*RCH3*RPHENOXM2+A110*T**n110*np.exp(-1*E110/R/T)*PRADIOM2*PRADIOM2+5.000000e-01*A53*T**n53*np.exp(-1*E53/R/T)*RCH3O*RPHENOXM2+2.000000e+00*A62*T**n62*np.exp(-1*E62/R/T)*RLIGM2A*RLIGM2A\n+5.000000e-01*A38*T**n38*np.exp(-1*E38/R/T)*RPHENOX*KETDM2;\n/*ODEforH2=>H2=25*/\ndydt[H2]=+A78*T**n78*np.exp(-1*E78/R/T)*PRFET3*PRFET3+A74*T**n74*np.exp(-1*E74/R/T)*RADIO*RADIO+5.000000e-01*A49*T**n49*np.exp(-1*E49/R/T)*C10H2M2*RPHENOXM2+A97*T**n97*np.exp(-1*E97/R/T)*PC2H2+\n1.500000e+00*A52*T**n52*np.exp(-1*E52/R/T)*RCH3O*RPHENOX+5.000000e-01*A48*T**n48*np.exp(-1*E48/R/T)*C10H2M4*RPHENOXM2+1.500000e+00*\nA45*T**n45*np.exp(-1*E45/R/T)*KET*RPHENOXM2+5.000000e-01*A32*T**n32*np.exp(-1*E32/R/T)*ADIOM2*RPHENOXM2+A96*T**n96*np.exp(-1*E96/R/T)*PH2+5.000000e-01*A83*T**n83*np.exp(-1*E83/R/T)*CHAR*RPHENOXM2+\n1.500000e+00*A41*T**n41*np.exp(-1*E41/R/T)*KET*RPHENOX+1.500000e+00*A40*T**n40*np.exp(-1*E40/R/T)*ADIO*RPHENOX+1.500000e+00*\nA80*T**n80*np.exp(-1*E80/R/T)*RPHENOX*RPHENOL+2.000000e+00*A75*T**n75*np.exp(-1*E75/R/T)*RLIGB*RLIGB+1.500000e+00*A44*T**n44*np.exp(-1*E44/R/T)*ADIO*RPHENOXM2+\n2.500000e+00*A43*T**n43*np.exp(-1*E43/R/T)*COUMARYL*RPHENOX+5.000000e-01*A33*T**n33*np.exp(-1*E33/R/T)*KETM2*RPHENOXM2+1.500000e+00*\nA47*T**n47*np.exp(-1*E47/R/T)*COUMARYL*RPHENOXM2+5.000000e-01*A102*T**n102*np.exp(-1*E102/R/T)*PCH2P+5.000000e-01*A53*T**n53*np.exp(-1*E53/R/T)*RCH3O*RPHENOXM2+\n5.000000e-01*A38*T**n38*np.exp(-1*E38/R/T)*RPHENOX*KETDM2+2.000000e+00*A73*T**n73*np.exp(-1*E73/R/T)*RLIGB*RPHENOX+5.000000e-01*\nA10*T**n10*np.exp(-1*E10/R/T)*RPHENOXM2+5.000000e-01*A103*T**n103*np.exp(-1*E103/R/T)*PCHP2+1.500000e+00*A39*T**n39*np.exp(-1*E39/R/T)*SYNAPYL*RPHENOX+\n2.000000e+00*A77*T**n77*np.exp(-1*E77/R/T)*RKET*RKET+1.500000e+00*A42*T**n42*np.exp(-1*E42/R/T)*KETD*RPHENOX+A81*T**n81*np.exp(-1*E81/R/T)*RC3H3O*RPHENOX+\n1.500000e+00*A50*T**n50*np.exp(-1*E50/R/T)*C10H2M4*RPHENOX+5.000000e-01*A37*T**n37*np.exp(-1*E37/R/T)*KETM2*RPHENOX+1.500000e+00*\nA46*T**n46*np.exp(-1*E46/R/T)*KETD*RPHENOXM2+A35*T**n35*np.exp(-1*E35/R/T)*SYNAPYL*RPHENOXM2+1.500000e+00*A11*T**n11*np.exp(-1*E11/R/T)*RPHENOX+A109*T**n109*np.exp(-1*E109/R/T)*PRADIO*PRADIO+5.000000e-01\n*A34*T**n34*np.exp(-1*E34/R/T)*RPHENOXM2*KETDM2+1.500000e+00*A51*T**n51*np.exp(-1*E51/R/T)*C10H2M2*RPHENOX+A79*T**n79*np.exp(-1*E79/R/T)*RLIGH*RLIGH+2.000000e+00*A76*T**n76*np.exp(-1*E76/R/T)*RLIGA*RLIGA\n+1.500000e+00*A82*T**n82*np.exp(-1*E82/R/T)*CHAR*RPHENOX;\n/*ODEforC10H2=>C10H2=26*/\ndydt[C10H2]=+A78*T**n78*np.exp(-1*E78/R/T)*PRFET3*PRFET3+A74*T**n74*np.exp(-1*E74/R/T)*RADIO*RADIO-5.000000e-01*A284*T**n284*np.exp(-1*E284/R/T)*PRFET3M2*C10H2+\n5.000000e-01*A49*T**n49*np.exp(-1*E49/R/T)*C10H2M2*RPHENOXM2-5.000000e-01*A292*T**n292*np.exp(-1*E292/R/T)*RC3H3O*C10H2-5.000000e-01*\nA283*T**n283*np.exp(-1*E283/R/T)*C10H2*RADIOM2-5.000000e-01*A282*T**n282*np.exp(-1*E282/R/T)*RC3H7O2*C10H2+5.000000e-01*A52*T**n52*np.exp(-1*E52/R/T)*RCH3O*RPHENOX+\n5.000000e-01*A45*T**n45*np.exp(-1*E45/R/T)*KET*RPHENOXM2-5.000000e-01*A287*T**n287*np.exp(-1*E287/R/T)*RLIGM2A*C10H2+A41*T**n41*np.exp(-1*E41/R/T)*KET*RPHENOX-\n5.000000e-01*A286*T**n286*np.exp(-1*E286/R/T)*RLIGM2B*C10H2+A40*T**n40*np.exp(-1*E40/R/T)*ADIO*RPHENOX-5.000000e-01*A299*T**n299*np.exp(-1*E299/R/T)*C10H2*RPHENOL+A80*T**n80*np.exp(-1*E80/R/T)*RPHENOX*RPHENOL\n+2.000000e+00*A75*T**n75*np.exp(-1*E75/R/T)*RLIGB*RLIGB+5.000000e-01*A44*T**n44*np.exp(-1*E44/R/T)*ADIO*RPHENOXM2+A43*T**n43*np.exp(-1*E43/R/T)*COUMARYL*RPHENOX-\n5.000000e-01*A285*T**n285*np.exp(-1*E285/R/T)*PRLIGH*C10H2+5.000000e-01*A47*T**n47*np.exp(-1*E47/R/T)*COUMARYL*RPHENOXM2-5.000000e-01*\nA288*T**n288*np.exp(-1*E288/R/T)*RCH3*C10H2+5.000000e-01*A38*T**n38*np.exp(-1*E38/R/T)*RPHENOX*KETDM2+1.500000e+00*A73*T**n73*np.exp(-1*E73/R/T)*RLIGB*RPHENOX+\n5.000000e-01*A29*T**n29*np.exp(-1*E29/R/T)*C10H2M2-5.000000e-01*A290*T**n290*np.exp(-1*E290/R/T)*RKET*C10H2+5.000000e-01*\nA39*T**n39*np.exp(-1*E39/R/T)*SYNAPYL*RPHENOX+A77*T**n77*np.exp(-1*E77/R/T)*RKET*RKET+A42*T**n42*np.exp(-1*E42/R/T)*KETD*RPHENOX-5.000000e-01*A297*T**n297*np.exp(-1*E297/R/T)*C10H2*OH+5.000000e-01\n*A36*T**n36*np.exp(-1*E36/R/T)*ADIOM2*RPHENOX+5.000000e-01*A81*T**n81*np.exp(-1*E81/R/T)*RC3H3O*RPHENOX-5.000000e-01*A289*T**n289*np.exp(-1*E289/R/T)*PRKETM2*C10H2+\n5.000000e-01*A50*T**n50*np.exp(-1*E50/R/T)*C10H2M4*RPHENOX-5.000000e-01*A280*T**n280*np.exp(-1*E280/R/T)*C10H2*RC3H5O2-5.000000e-01*\nA295*T**n295*np.exp(-1*E295/R/T)*PRADIOM2*C10H2+5.000000e-01*A37*T**n37*np.exp(-1*E37/R/T)*KETM2*RPHENOX-5.000000e-01*A294*T**n294*np.exp(-1*E294/R/T)*RLIGA*C10H2+\n5.000000e-01*A46*T**n46*np.exp(-1*E46/R/T)*KETD*RPHENOXM2-5.000000e-01*A296*T**n296*np.exp(-1*E296/R/T)*RMGUAI*C10H2+5.000000e-01*\nA11*T**n11*np.exp(-1*E11/R/T)*RPHENOX+A109*T**n109*np.exp(-1*E109/R/T)*PRADIO*PRADIO-5.000000e-01*A298*T**n298*np.exp(-1*E298/R/T)*RCH3O*C10H2-5.000000e-01*\nA281*T**n281*np.exp(-1*E281/R/T)*C10H2*PRFET3+A51*T**n51*np.exp(-1*E51/R/T)*C10H2M2*RPHENOX-5.000000e-01*A293*T**n293*np.exp(-1*E293/R/T)*RLIGB*C10H2+2.000000e+00*A76*T**n76*np.exp(-1*E76/R/T)*RLIGA*RLIGA\n-5.000000e-01*A291*T**n291*np.exp(-1*E291/R/T)*C10H2*PRADIO+5.000000e-01*A55*T**n55*np.exp(-1*E55/R/T)*RCH3*RPHENOX+5.000000e-01*\nA82*T**n82*np.exp(-1*E82/R/T)*CHAR*RPHENOX-5.000000e-01*A300*T**n300*np.exp(-1*E300/R/T)*RADIO*C10H2;\n/*ODEforRLIGH=>RLIGH=27*/\ndydt[RLIGH]=+A123*T**n123*np.exp(-1*E123/R/T)*LIGH*PRADIO+A132*T**n132*np.exp(-1*E132/R/T)*RADIO*LIGH+A129*T**n129*np.exp(-1*E129/R/T)*LIGH*OH+A115*T**n115*np.exp(-1*E115/R/T)*LIGH*RADIOM2+A122*T**n122*np.exp(-1*E122/R/T)*LIGH*RKET+\nA127*T**n127*np.exp(-1*E127/R/T)*PRADIOM2*LIGH+A120*T**n120*np.exp(-1*E120/R/T)*RCH3*LIGH-A12*T**n12*np.exp(-1*E12/R/T)*RLIGH+A114*T**n114*np.exp(-1*E114/R/T)*RC3H7O2*LIGH+A126*T**n126*np.exp(-1*E126/R/T)*LIGH*RLIGA+A116*T**n116*np.exp(-1*E116/R/T)*PRFET3M2*LIGH+\nA118*T**n118*np.exp(-1*E118/R/T)*RLIGM2B*LIGH+A121*T**n121*np.exp(-1*E121/R/T)*PRKETM2*LIGH+A124*T**n124*np.exp(-1*E124/R/T)*RC3H3O*LIGH+A112*T**n112*np.exp(-1*E112/R/T)*LIGH*RC3H5O2+A125*T**n125*np.exp(-1*E125/R/T)*RLIGB*LIGH+A113*T**n113*np.exp(-1*E113/R/T)*LIGH*PRFET3+\nA119*T**n119*np.exp(-1*E119/R/T)*RLIGM2A*LIGH-2.000000e+00*A79*T**n79*np.exp(-1*E79/R/T)*RLIGH*RLIGH+A128*T**n128*np.exp(-1*E128/R/T)*RMGUAI*LIGH+A131*T**n131*np.exp(-1*E131/R/T)*LIGH*RPHENOL+A117*T**n117*np.exp(-1*E117/R/T)*PRLIGH*LIGH+\nA130*T**n130*np.exp(-1*E130/R/T)*RCH3O*LIGH;\n/*ODEforALD3=>ALD3=28*/\ndydt[ALD3]=+A13*T**n13*np.exp(-1*E13/R/T)*PRLIGH2+A12*T**n12*np.exp(-1*E12/R/T)*RLIGH;\n/*ODEforPRLIGH2=>PRLIGH2=29*/\ndydt[PRLIGH2]=+A148*T**n148*np.exp(-1*E148/R/T)*PRADIOM2*PLIGH+A146*T**n146*np.exp(-1*E146/R/T)*PLIGH*RLIGB+A153*T**n153*np.exp(-1*E153/R/T)*RADIO*PLIGH-A13*T**n13*np.exp(-1*E13/R/T)*PRLIGH2+A137*T**n137*np.exp(-1*E137/R/T)*PLIGH*PRFET3M2+\nA133*T**n133*np.exp(-1*E133/R/T)*PLIGH*RC3H5O2+A136*T**n136*np.exp(-1*E136/R/T)*PLIGH*RADIOM2+A142*T**n142*np.exp(-1*E142/R/T)*PRKETM2*PLIGH+A140*T**n140*np.exp(-1*E140/R/T)*RLIGM2A*PLIGH+A135*T**n135*np.exp(-1*E135/R/T)*RC3H7O2*PLIGH+A141*T**n141*np.exp(-1*E141/R/T)*PLIGH*RCH3+\nA134*T**n134*np.exp(-1*E134/R/T)*PLIGH*PRFET3+A138*T**n138*np.exp(-1*E138/R/T)*PLIGH*PRLIGH+A139*T**n139*np.exp(-1*E139/R/T)*RLIGM2B*PLIGH+A145*T**n145*np.exp(-1*E145/R/T)*RC3H3O*PLIGH+A143*T**n143*np.exp(-1*E143/R/T)*PLIGH*RKET+A147*T**n147*np.exp(-1*E147/R/T)*PLIGH*RLIGA+\nA152*T**n152*np.exp(-1*E152/R/T)*PLIGH*RPHENOL+A151*T**n151*np.exp(-1*E151/R/T)*PLIGH*RCH3O+A149*T**n149*np.exp(-1*E149/R/T)*RMGUAI*PLIGH+A144*T**n144*np.exp(-1*E144/R/T)*PLIGH*PRADIO+A150*T**n150*np.exp(-1*E150/R/T)*PLIGH*OH;\n/*ODEforPRLIGM2A=>PRLIGM2A=30*/\ndydt[PRLIGM2A]=+A167*T**n167*np.exp(-1*E167/R/T)*PLIGM2*RLIGB-A16*T**n16*np.exp(-1*E16/R/T)*PRLIGM2A+A13*T**n13*np.exp(-1*E13/R/T)*PRLIGH2+A173*T**n173*np.exp(-1*E173/R/T)*PLIGM2*RPHENOL+A161*T**n161*np.exp(-1*E161/R/T)*PLIGM2*RLIGM2A+\nA168*T**n168*np.exp(-1*E168/R/T)*PLIGM2*RLIGA+A172*T**n172*np.exp(-1*E172/R/T)*PLIGM2*RCH3O+A169*T**n169*np.exp(-1*E169/R/T)*PLIGM2*PRADIOM2+A158*T**n158*np.exp(-1*E158/R/T)*PLIGM2*PRFET3M2+A171*T**n171*np.exp(-1*E171/R/T)*PLIGM2*OH+A157*T**n157*np.exp(-1*E157/R/T)*PLIGM2*RADIOM2+\nA170*T**n170*np.exp(-1*E170/R/T)*PLIGM2*RMGUAI+A155*T**n155*np.exp(-1*E155/R/T)*PLIGM2*PRFET3+A156*T**n156*np.exp(-1*E156/R/T)*PLIGM2*RC3H7O2+A165*T**n165*np.exp(-1*E165/R/T)*PLIGM2*PRADIO+A174*T**n174*np.exp(-1*E174/R/T)*RADIO*PLIGM2+A164*T**n164*np.exp(-1*E164/R/T)*PLIGM2*RKET+\nA154*T**n154*np.exp(-1*E154/R/T)*PLIGM2*RC3H5O2+A162*T**n162*np.exp(-1*E162/R/T)*PLIGM2*RCH3+A160*T**n160*np.exp(-1*E160/R/T)*RLIGM2B*PLIGM2+A159*T**n159*np.exp(-1*E159/R/T)*PLIGM2*PRLIGH+A166*T**n166*np.exp(-1*E166/R/T)*PLIGM2*RC3H3O+A163*T**n163*np.exp(-1*E163/R/T)*PLIGM2*PRKETM2;\n/*ODEforRMGUAI=>RMGUAI=31*/\ndydt[RMGUAI]=-A191*T**n191*np.exp(-1*E191/R/T)*RMGUAI*LIGM2+A14*T**n14*np.exp(-1*E14/R/T)*RADIOM2-A401*T**n401*np.exp(-1*E401/R/T)*RMGUAI*KET-2.000000e+00*A63*T**n63*np.exp(-1*E63/R/T)*RMGUAI*RMGUAI\n-A233*T**n233*np.exp(-1*E233/R/T)*RMGUAI*PFET3M2-A212*T**n212*np.exp(-1*E212/R/T)*RMGUAI*LIGM2-A338*T**n338*np.exp(-1*E338/R/T)*LIG*RMGUAI-A170*T**n170*np.exp(-1*E170/R/T)*PLIGM2*RMGUAI-A296*T**n296*np.exp(-1*E296/R/T)*RMGUAI*C10H2-A254*T**n254*np.exp(-1*E254/R/T)*RMGUAI*ADIOM2-\nA317*T**n317*np.exp(-1*E317/R/T)*LIG*RMGUAI-A128*T**n128*np.exp(-1*E128/R/T)*RMGUAI*LIGH+A17*T**n17*np.exp(-1*E17/R/T)*RLIGM2B-A359*T**n359*np.exp(-1*E359/R/T)*PFET3*RMGUAI-A380*T**n380*np.exp(-1*E380/R/T)*RMGUAI*ADIO-A275*T**n275*np.exp(-1*E275/R/T)*KETM2*RMGUAI-\nA149*T**n149*np.exp(-1*E149/R/T)*RMGUAI*PLIGH;\n/*ODEforC3H6O2=>C3H6O2=32*/\ndydt[C3H6O2]=+A301*T**n301*np.exp(-1*E301/R/T)*LIG*RC3H5O2+A19*T**n19*np.exp(-1*E19/R/T)*RADIO+A14*T**n14*np.exp(-1*E14/R/T)*RADIOM2+A133*T**n133*np.exp(-1*E133/R/T)*PLIGH*RC3H5O2+A364*T**n364*np.exp(-1*E364/R/T)*ADIO*RC3H5O2+\nA343*T**n343*np.exp(-1*E343/R/T)*PFET3*RC3H5O2+A280*T**n280*np.exp(-1*E280/R/T)*C10H2*RC3H5O2+A217*T**n217*np.exp(-1*E217/R/T)*RC3H5O2*PFET3M2+A259*T**n259*np.exp(-1*E259/R/T)*KETM2*RC3H5O2+A175*T**n175*np.exp(-1*E175/R/T)*RC3H5O2*LIGM2+A112*T**n112*np.exp(-1*E112/R/T)*LIGH*RC3H5O2+\nA385*T**n385*np.exp(-1*E385/R/T)*KET*RC3H5O2+A196*T**n196*np.exp(-1*E196/R/T)*RC3H5O2*LIGM2+A322*T**n322*np.exp(-1*E322/R/T)*LIG*RC3H5O2+A154*T**n154*np.exp(-1*E154/R/T)*PLIGM2*RC3H5O2+A238*T**n238*np.exp(-1*E238/R/T)*RC3H5O2*ADIOM2;\n/*ODEforKETM2=>KETM2=33*/\ndydt[KETM2]=-A260*T**n260*np.exp(-1*E260/R/T)*KETM2*PRFET3-A278*T**n278*np.exp(-1*E278/R/T)*KETM2*RPHENOL-A269*T**n269*np.exp(-1*E269/R/T)*KETM2*RKET-A279*T**n279*np.exp(-1*E279/R/T)*RADIO*KETM2-A265*T**n265*np.exp(-1*E265/R/T)*RLIGM2B*KETM2-\nA270*T**n270*np.exp(-1*E270/R/T)*KETM2*PRADIO-A85*T**n85*np.exp(-1*E85/R/T)*KETM2-A37*T**n37*np.exp(-1*E37/R/T)*KETM2*RPHENOX-A264*T**n264*np.exp(-1*E264/R/T)*KETM2*PRLIGH-A276*T**n276*np.exp(-1*E276/R/T)*KETM2*OH-A277*T**n277*np.exp(-1*E277/R/T)*KETM2*RCH3O-\nA33*T**n33*np.exp(-1*E33/R/T)*KETM2*RPHENOXM2-A259*T**n259*np.exp(-1*E259/R/T)*KETM2*RC3H5O2-A268*T**n268*np.exp(-1*E268/R/T)*KETM2*PRKETM2-A274*T**n274*np.exp(-1*E274/R/T)*KETM2*PRADIOM2-A263*T**n263*np.exp(-1*E263/R/T)*KETM2*PRFET3M2-A273*T**n273*np.exp(-1*E273/R/T)*KETM2*RLIGA-\nA272*T**n272*np.exp(-1*E272/R/T)*KETM2*RLIGB-A262*T**n262*np.exp(-1*E262/R/T)*KETM2*RADIOM2-A267*T**n267*np.exp(-1*E267/R/T)*KETM2*RCH3+A15*T**n15*np.exp(-1*E15/R/T)*RLIGM2A-A271*T**n271*np.exp(-1*E271/R/T)*KETM2*RC3H3O-A275*T**n275*np.exp(-1*E275/R/T)*KETM2*RMGUAI-\nA261*T**n261*np.exp(-1*E261/R/T)*KETM2*RC3H7O2-A266*T**n266*np.exp(-1*E266/R/T)*KETM2*RLIGM2A;\n/*ODEforRLIGM2B=>RLIGM2B=34*/\ndydt[RLIGM2B]=+A216*T**n216*np.exp(-1*E216/R/T)*RADIO*LIGM2+A201*T**n201*np.exp(-1*E201/R/T)*PRLIGH*LIGM2+A211*T**n211*np.exp(-1*E211/R/T)*PRADIOM2*LIGM2+A210*T**n210*np.exp(-1*E210/R/T)*RLIGA*LIGM2-A286*T**n286*np.exp(-1*E286/R/T)*RLIGM2B*C10H2\n+A204*T**n204*np.exp(-1*E204/R/T)*RCH3*LIGM2+A212*T**n212*np.exp(-1*E212/R/T)*RMGUAI*LIGM2-A139*T**n139*np.exp(-1*E139/R/T)*RLIGM2B*PLIGH+A199*T**n199*np.exp(-1*E199/R/T)*LIGM2*RADIOM2+A215*T**n215*np.exp(-1*E215/R/T)*LIGM2*RPHENOL-A307*T**n307*np.exp(-1*E307/R/T)*RLIGM2B*LIG-\nA244*T**n244*np.exp(-1*E244/R/T)*RLIGM2B*ADIOM2-2.000000e+00*A61*T**n61*np.exp(-1*E61/R/T)*RLIGM2B*RLIGM2B+A196*T**n196*np.exp(-1*E196/R/T)*RC3H5O2*LIGM2+A206*T**n206*np.exp(-1*E206/R/T)*RKET*LIGM2-A328*T**n328*np.exp(-1*E328/R/T)*RLIGM2B*LIG-\nA223*T**n223*np.exp(-1*E223/R/T)*RLIGM2B*PFET3M2-A265*T**n265*np.exp(-1*E265/R/T)*RLIGM2B*KETM2-A370*T**n370*np.exp(-1*E370/R/T)*RLIGM2B*ADIO+A213*T**n213*np.exp(-1*E213/R/T)*LIGM2*OH+A197*T**n197*np.exp(-1*E197/R/T)*PRFET3*LIGM2+A208*T**n208*np.exp(-1*E208/R/T)*RC3H3O*LIGM2+\nA203*T**n203*np.exp(-1*E203/R/T)*RLIGM2A*LIGM2+A200*T**n200*np.exp(-1*E200/R/T)*PRFET3M2*LIGM2-A118*T**n118*np.exp(-1*E118/R/T)*RLIGM2B*LIGH-A181*T**n181*np.exp(-1*E181/R/T)*RLIGM2B*LIGM2-A349*T**n349*np.exp(-1*E349/R/T)*RLIGM2B*PFET3+A207*T**n207*np.exp(-1*E207/R/T)*PRADIO*LIGM2+\nA209*T**n209*np.exp(-1*E209/R/T)*RLIGB*LIGM2-A391*T**n391*np.exp(-1*E391/R/T)*RLIGM2B*KET+A205*T**n205*np.exp(-1*E205/R/T)*PRKETM2*LIGM2-A17*T**n17*np.exp(-1*E17/R/T)*RLIGM2B+A198*T**n198*np.exp(-1*E198/R/T)*RC3H7O2*LIGM2-A160*T**n160*np.exp(-1*E160/R/T)*RLIGM2B*PLIGM2+\nA214*T**n214*np.exp(-1*E214/R/T)*RCH3O*LIGM2;\n/*ODEforPFET3M2=>PFET3M2=35*/\ndydt[PFET3M2]=-A229*T**n229*np.exp(-1*E229/R/T)*RC3H3O*PFET3M2-A219*T**n219*np.exp(-1*E219/R/T)*RC3H7O2*PFET3M2+A284*T**n284*np.exp(-1*E284/R/T)*PRFET3M2*C10H2-A230*T**n230*np.exp(-1*E230/R/T)*RLIGB*PFET3M2+A137*T**n137*np.exp(-1*E137/R/T)*PLIGH*PRFET3M2\n+A347*T**n347*np.exp(-1*E347/R/T)*PFET3*PRFET3M2-A236*T**n236*np.exp(-1*E236/R/T)*PFET3M2*RPHENOL-A232*T**n232*np.exp(-1*E232/R/T)*PRADIOM2*PFET3M2-A235*T**n235*np.exp(-1*E235/R/T)*RCH3O*PFET3M2-A233*T**n233*np.exp(-1*E233/R/T)*RMGUAI*PFET3M2+A158*T**n158*np.exp(-1*E158/R/T)*PLIGM2*PRFET3M2+\nA326*T**n326*np.exp(-1*E326/R/T)*LIG*PRFET3M2-A217*T**n217*np.exp(-1*E217/R/T)*RC3H5O2*PFET3M2+A263*T**n263*np.exp(-1*E263/R/T)*KETM2*PRFET3M2+A305*T**n305*np.exp(-1*E305/R/T)*LIG*PRFET3M2-A226*T**n226*np.exp(-1*E226/R/T)*PRKETM2*PFET3M2-A220*T**n220*np.exp(-1*E220/R/T)*PFET3M2*RADIOM2-\nA225*T**n225*np.exp(-1*E225/R/T)*RCH3*PFET3M2+A242*T**n242*np.exp(-1*E242/R/T)*PRFET3M2*ADIOM2-A227*T**n227*np.exp(-1*E227/R/T)*RKET*PFET3M2-A223*T**n223*np.exp(-1*E223/R/T)*RLIGM2B*PFET3M2+A389*T**n389*np.exp(-1*E389/R/T)*PRFET3M2*KET-A234*T**n234*np.exp(-1*E234/R/T)*OH*PFET3M2-\nA218*T**n218*np.exp(-1*E218/R/T)*PRFET3*PFET3M2-A224*T**n224*np.exp(-1*E224/R/T)*RLIGM2A*PFET3M2-A222*T**n222*np.exp(-1*E222/R/T)*PRLIGH*PFET3M2+A200*T**n200*np.exp(-1*E200/R/T)*PRFET3M2*LIGM2+A116*T**n116*np.exp(-1*E116/R/T)*PRFET3M2*LIGH-A237*T**n237*np.exp(-1*E237/R/T)*RADIO*PFET3M2+\nA179*T**n179*np.exp(-1*E179/R/T)*PRFET3M2*LIGM2+A368*T**n368*np.exp(-1*E368/R/T)*PRFET3M2*ADIO+A17*T**n17*np.exp(-1*E17/R/T)*RLIGM2B-A228*T**n228*np.exp(-1*E228/R/T)*PRADIO*PFET3M2-A231*T**n231*np.exp(-1*E231/R/T)*RLIGA*PFET3M2;\n/*ODEforPRFET3M2=>PRFET3M2=36*/\ndydt[PRFET3M2]=+A229*T**n229*np.exp(-1*E229/R/T)*RC3H3O*PFET3M2+A219*T**n219*np.exp(-1*E219/R/T)*RC3H7O2*PFET3M2-A284*T**n284*np.exp(-1*E284/R/T)*PRFET3M2*C10H2+A230*T**n230*np.exp(-1*E230/R/T)*RLIGB*PFET3M2-A137*T**n137*np.exp(-1*E137/R/T)*PLIGH*PRFET3M2\n-A347*T**n347*np.exp(-1*E347/R/T)*PFET3*PRFET3M2+A236*T**n236*np.exp(-1*E236/R/T)*PFET3M2*RPHENOL+A232*T**n232*np.exp(-1*E232/R/T)*PRADIOM2*PFET3M2+A235*T**n235*np.exp(-1*E235/R/T)*RCH3O*PFET3M2+A233*T**n233*np.exp(-1*E233/R/T)*RMGUAI*PFET3M2-A158*T**n158*np.exp(-1*E158/R/T)*PLIGM2*PRFET3M2-\nA326*T**n326*np.exp(-1*E326/R/T)*LIG*PRFET3M2+A217*T**n217*np.exp(-1*E217/R/T)*RC3H5O2*PFET3M2-A263*T**n263*np.exp(-1*E263/R/T)*KETM2*PRFET3M2-2.000000e+00*A65*T**n65*np.exp(-1*E65/R/T)*PRFET3M2*PRFET3M2-A305*T**n305*np.exp(-1*E305/R/T)*LIG*PRFET3M2+\nA226*T**n226*np.exp(-1*E226/R/T)*PRKETM2*PFET3M2+A220*T**n220*np.exp(-1*E220/R/T)*PFET3M2*RADIOM2+A225*T**n225*np.exp(-1*E225/R/T)*RCH3*PFET3M2-A18*T**n18*np.exp(-1*E18/R/T)*PRFET3M2-A242*T**n242*np.exp(-1*E242/R/T)*PRFET3M2*ADIOM2+A227*T**n227*np.exp(-1*E227/R/T)*RKET*PFET3M2+\nA223*T**n223*np.exp(-1*E223/R/T)*RLIGM2B*PFET3M2-A389*T**n389*np.exp(-1*E389/R/T)*PRFET3M2*KET+A234*T**n234*np.exp(-1*E234/R/T)*OH*PFET3M2+A218*T**n218*np.exp(-1*E218/R/T)*PRFET3*PFET3M2+A224*T**n224*np.exp(-1*E224/R/T)*RLIGM2A*PFET3M2+A222*T**n222*np.exp(-1*E222/R/T)*PRLIGH*PFET3M2-\nA200*T**n200*np.exp(-1*E200/R/T)*PRFET3M2*LIGM2-A116*T**n116*np.exp(-1*E116/R/T)*PRFET3M2*LIGH+A237*T**n237*np.exp(-1*E237/R/T)*RADIO*PFET3M2-A179*T**n179*np.exp(-1*E179/R/T)*PRFET3M2*LIGM2-A368*T**n368*np.exp(-1*E368/R/T)*PRFET3M2*ADIO+A228*T**n228*np.exp(-1*E228/R/T)*PRADIO*PFET3M2+\nA231*T**n231*np.exp(-1*E231/R/T)*RLIGA*PFET3M2;\n/*ODEforC3H4O2=>C3H4O2=37*/\ndydt[C3H4O2]=+A18*T**n18*np.exp(-1*E18/R/T)*PRFET3M2+A22*T**n22*np.exp(-1*E22/R/T)*PRFET3;\n/*ODEforRPHENOL=>RPHENOL=38*/\ndydt[RPHENOL]=+A19*T**n19*np.exp(-1*E19/R/T)*RADIO-A278*T**n278*np.exp(-1*E278/R/T)*KETM2*RPHENOL-A173*T**n173*np.exp(-1*E173/R/T)*PLIGM2*RPHENOL-A383*T**n383*np.exp(-1*E383/R/T)*ADIO*RPHENOL-A320*T**n320*np.exp(-1*E320/R/T)*LIG*RPHENOL-\nA236*T**n236*np.exp(-1*E236/R/T)*PFET3M2*RPHENOL-A194*T**n194*np.exp(-1*E194/R/T)*LIGM2*RPHENOL-A299*T**n299*np.exp(-1*E299/R/T)*C10H2*RPHENOL-A80*T**n80*np.exp(-1*E80/R/T)*RPHENOX*RPHENOL-A362*T**n362*np.exp(-1*E362/R/T)*PFET3*RPHENOL-A215*T**n215*np.exp(-1*E215/R/T)*LIGM2*RPHENOL-\nA341*T**n341*np.exp(-1*E341/R/T)*LIG*RPHENOL-A257*T**n257*np.exp(-1*E257/R/T)*ADIOM2*RPHENOL-A131*T**n131*np.exp(-1*E131/R/T)*LIGH*RPHENOL+A21*T**n21*np.exp(-1*E21/R/T)*RLIGB-A404*T**n404*np.exp(-1*E404/R/T)*KET*RPHENOL-A152*T**n152*np.exp(-1*E152/R/T)*PLIGH*RPHENOL;\n/*ODEforRLIGA=>RLIGA=39*/\ndydt[RLIGA]=+A308*T**n308*np.exp(-1*E308/R/T)*LIG*RLIGM2A+A316*T**n316*np.exp(-1*E316/R/T)*LIG*PRADIOM2+A306*T**n306*np.exp(-1*E306/R/T)*LIG*PRLIGH-A189*T**n189*np.exp(-1*E189/R/T)*RLIGA*LIGM2-A210*T**n210*np.exp(-1*E210/R/T)*RLIGA*LIGM2+\nA318*T**n318*np.exp(-1*E318/R/T)*LIG*OH+A310*T**n310*np.exp(-1*E310/R/T)*LIG*PRKETM2-A357*T**n357*np.exp(-1*E357/R/T)*PFET3*RLIGA+A319*T**n319*np.exp(-1*E319/R/T)*LIG*RCH3O+A304*T**n304*np.exp(-1*E304/R/T)*LIG*RADIOM2+A307*T**n307*np.exp(-1*E307/R/T)*RLIGM2B*LIG+\nA317*T**n317*np.exp(-1*E317/R/T)*LIG*RMGUAI+A311*T**n311*np.exp(-1*E311/R/T)*LIG*RKET-A273*T**n273*np.exp(-1*E273/R/T)*KETM2*RLIGA+A305*T**n305*np.exp(-1*E305/R/T)*LIG*PRFET3M2+A312*T**n312*np.exp(-1*E312/R/T)*LIG*PRADIO-A147*T**n147*np.exp(-1*E147/R/T)*PLIGH*RLIGA+\nA303*T**n303*np.exp(-1*E303/R/T)*LIG*RC3H7O2-A336*T**n336*np.exp(-1*E336/R/T)*LIG*RLIGA+A321*T**n321*np.exp(-1*E321/R/T)*RADIO*LIG-A20*T**n20*np.exp(-1*E20/R/T)*RLIGA+A301*T**n301*np.exp(-1*E301/R/T)*LIG*RC3H5O2+A309*T**n309*np.exp(-1*E309/R/T)*LIG*RCH3-\nA168*T**n168*np.exp(-1*E168/R/T)*PLIGM2*RLIGA+A320*T**n320*np.exp(-1*E320/R/T)*LIG*RPHENOL-A378*T**n378*np.exp(-1*E378/R/T)*ADIO*RLIGA-A294*T**n294*np.exp(-1*E294/R/T)*RLIGA*C10H2-A252*T**n252*np.exp(-1*E252/R/T)*RLIGA*ADIOM2-A126*T**n126*np.exp(-1*E126/R/T)*LIGH*RLIGA+\nA313*T**n313*np.exp(-1*E313/R/T)*LIG*RC3H3O-A231*T**n231*np.exp(-1*E231/R/T)*RLIGA*PFET3M2+A302*T**n302*np.exp(-1*E302/R/T)*LIG*PRFET3-2.000000e+00*A76*T**n76*np.exp(-1*E76/R/T)*RLIGA*RLIGA+A314*T**n314*np.exp(-1*E314/R/T)*LIG*RLIGB-\nA399*T**n399*np.exp(-1*E399/R/T)*KET*RLIGA;\n/*ODEforKET=>KET=40*/\ndydt[KET]=-A91*T**n91*np.exp(-1*E91/R/T)*KET-A394*T**n394*np.exp(-1*E394/R/T)*PRKETM2*KET-A45*T**n45*np.exp(-1*E45/R/T)*KET*RPHENOXM2-A41*T**n41*np.exp(-1*E41/R/T)*KET*RPHENOX-A392*T**n392*np.exp(-1*E392/R/T)*RLIGM2A*KET-\nA398*T**n398*np.exp(-1*E398/R/T)*RLIGB*KET-A405*T**n405*np.exp(-1*E405/R/T)*RADIO*KET+A185*T**n185*np.exp(-1*E185/R/T)*RKET*LIGM2+A143*T**n143*np.exp(-1*E143/R/T)*PLIGH*RKET+A311*T**n311*np.exp(-1*E311/R/T)*LIG*RKET+A164*T**n164*np.exp(-1*E164/R/T)*PLIGM2*RKET-\nA400*T**n400*np.exp(-1*E400/R/T)*PRADIOM2*KET+A206*T**n206*np.exp(-1*E206/R/T)*RKET*LIGM2+A227*T**n227*np.exp(-1*E227/R/T)*RKET*PFET3M2-A402*T**n402*np.exp(-1*E402/R/T)*KET*OH-A387*T**n387*np.exp(-1*E387/R/T)*RC3H7O2*KET+A20*T**n20*np.exp(-1*E20/R/T)*RLIGA+\nA269*T**n269*np.exp(-1*E269/R/T)*KETM2*RKET-A389*T**n389*np.exp(-1*E389/R/T)*PRFET3M2*KET+A290*T**n290*np.exp(-1*E290/R/T)*RKET*C10H2-A401*T**n401*np.exp(-1*E401/R/T)*RMGUAI*KET-A397*T**n397*np.exp(-1*E397/R/T)*RC3H3O*KET+A122*T**n122*np.exp(-1*E122/R/T)*LIGH*RKET-\nA388*T**n388*np.exp(-1*E388/R/T)*KET*RADIOM2-A396*T**n396*np.exp(-1*E396/R/T)*KET*PRADIO+A332*T**n332*np.exp(-1*E332/R/T)*LIG*RKET-A393*T**n393*np.exp(-1*E393/R/T)*RCH3*KET-A390*T**n390*np.exp(-1*E390/R/T)*PRLIGH*KET-A385*T**n385*np.exp(-1*E385/R/T)*KET*RC3H5O2-\nA386*T**n386*np.exp(-1*E386/R/T)*KET*PRFET3+A353*T**n353*np.exp(-1*E353/R/T)*PFET3*RKET+A248*T**n248*np.exp(-1*E248/R/T)*RKET*ADIOM2-A391*T**n391*np.exp(-1*E391/R/T)*RLIGM2B*KET-A404*T**n404*np.exp(-1*E404/R/T)*KET*RPHENOL-A403*T**n403*np.exp(-1*E403/R/T)*RCH3O*KET+\nA374*T**n374*np.exp(-1*E374/R/T)*ADIO*RKET-A399*T**n399*np.exp(-1*E399/R/T)*KET*RLIGA;\n/*ODEforRLIGB=>RLIGB=41*/\ndydt[RLIGB]=-A230*T**n230*np.exp(-1*E230/R/T)*RLIGB*PFET3M2-A356*T**n356*np.exp(-1*E356/R/T)*PFET3*RLIGB+A339*T**n339*np.exp(-1*E339/R/T)*LIG*OH+A327*T**n327*np.exp(-1*E327/R/T)*LIG*PRLIGH-A377*T**n377*np.exp(-1*E377/R/T)*RLIGB*ADIO-\n2.000000e+00*A75*T**n75*np.exp(-1*E75/R/T)*RLIGB*RLIGB-A398*T**n398*np.exp(-1*E398/R/T)*RLIGB*KET+A326*T**n326*np.exp(-1*E326/R/T)*LIG*PRFET3M2+A334*T**n334*np.exp(-1*E334/R/T)*LIG*RC3H3O+A336*T**n336*np.exp(-1*E336/R/T)*LIG*RLIGA-\nA73*T**n73*np.exp(-1*E73/R/T)*RLIGB*RPHENOX+A342*T**n342*np.exp(-1*E342/R/T)*RADIO*LIG+A328*T**n328*np.exp(-1*E328/R/T)*RLIGM2B*LIG-A146*T**n146*np.exp(-1*E146/R/T)*PLIGH*RLIGB+A323*T**n323*np.exp(-1*E323/R/T)*LIG*PRFET3-A167*T**n167*np.exp(-1*E167/R/T)*PLIGM2*RLIGB+\nA330*T**n330*np.exp(-1*E330/R/T)*LIG*RCH3+A333*T**n333*np.exp(-1*E333/R/T)*LIG*PRADIO+A325*T**n325*np.exp(-1*E325/R/T)*LIG*RADIOM2+A331*T**n331*np.exp(-1*E331/R/T)*LIG*PRKETM2-A188*T**n188*np.exp(-1*E188/R/T)*RLIGB*LIGM2+A324*T**n324*np.exp(-1*E324/R/T)*LIG*RC3H7O2+\nA332*T**n332*np.exp(-1*E332/R/T)*LIG*RKET+A338*T**n338*np.exp(-1*E338/R/T)*LIG*RMGUAI-A125*T**n125*np.exp(-1*E125/R/T)*RLIGB*LIGH+A329*T**n329*np.exp(-1*E329/R/T)*LIG*RLIGM2A+A341*T**n341*np.exp(-1*E341/R/T)*LIG*RPHENOL+A322*T**n322*np.exp(-1*E322/R/T)*LIG*RC3H5O2-\nA209*T**n209*np.exp(-1*E209/R/T)*RLIGB*LIGM2-A272*T**n272*np.exp(-1*E272/R/T)*KETM2*RLIGB-A21*T**n21*np.exp(-1*E21/R/T)*RLIGB-A293*T**n293*np.exp(-1*E293/R/T)*RLIGB*C10H2+A337*T**n337*np.exp(-1*E337/R/T)*LIG*PRADIOM2+A340*T**n340*np.exp(-1*E340/R/T)*LIG*RCH3O-\nA314*T**n314*np.exp(-1*E314/R/T)*LIG*RLIGB-A251*T**n251*np.exp(-1*E251/R/T)*RLIGB*ADIOM2;\n/*ODEforPFET3=>PFET3=42*/\ndydt[PFET3]=+A260*T**n260*np.exp(-1*E260/R/T)*KETM2*PRFET3-A361*T**n361*np.exp(-1*E361/R/T)*PFET3*RCH3O-A356*T**n356*np.exp(-1*E356/R/T)*PFET3*RLIGB-A347*T**n347*np.exp(-1*E347/R/T)*PFET3*PRFET3M2-A343*T**n343*np.exp(-1*E343/R/T)*PFET3*RC3H5O2-\nA357*T**n357*np.exp(-1*E357/R/T)*PFET3*RLIGA-A362*T**n362*np.exp(-1*E362/R/T)*PFET3*RPHENOL+A155*T**n155*np.exp(-1*E155/R/T)*PLIGM2*PRFET3-A350*T**n350*np.exp(-1*E350/R/T)*RLIGM2A*PFET3-A354*T**n354*np.exp(-1*E354/R/T)*PFET3*PRADIO-A345*T**n345*np.exp(-1*E345/R/T)*RC3H7O2*PFET3-\nA359*T**n359*np.exp(-1*E359/R/T)*PFET3*RMGUAI-A358*T**n358*np.exp(-1*E358/R/T)*PRADIOM2*PFET3-A355*T**n355*np.exp(-1*E355/R/T)*RC3H3O*PFET3-A346*T**n346*np.exp(-1*E346/R/T)*PFET3*RADIOM2+A323*T**n323*np.exp(-1*E323/R/T)*LIG*PRFET3+A176*T**n176*np.exp(-1*E176/R/T)*PRFET3*LIGM2+\nA365*T**n365*np.exp(-1*E365/R/T)*ADIO*PRFET3-A348*T**n348*np.exp(-1*E348/R/T)*PFET3*PRLIGH+A218*T**n218*np.exp(-1*E218/R/T)*PRFET3*PFET3M2+A134*T**n134*np.exp(-1*E134/R/T)*PLIGH*PRFET3+A239*T**n239*np.exp(-1*E239/R/T)*PRFET3*ADIOM2+A197*T**n197*np.exp(-1*E197/R/T)*PRFET3*LIGM2-\nA352*T**n352*np.exp(-1*E352/R/T)*PRKETM2*PFET3+A386*T**n386*np.exp(-1*E386/R/T)*KET*PRFET3+A281*T**n281*np.exp(-1*E281/R/T)*C10H2*PRFET3+A113*T**n113*np.exp(-1*E113/R/T)*LIGH*PRFET3-A353*T**n353*np.exp(-1*E353/R/T)*PFET3*RKET-A349*T**n349*np.exp(-1*E349/R/T)*RLIGM2B*PFET3-\nA351*T**n351*np.exp(-1*E351/R/T)*PFET3*RCH3+A21*T**n21*np.exp(-1*E21/R/T)*RLIGB-A360*T**n360*np.exp(-1*E360/R/T)*PFET3*OH+A302*T**n302*np.exp(-1*E302/R/T)*LIG*PRFET3-A363*T**n363*np.exp(-1*E363/R/T)*RADIO*PFET3;\n/*ODEforPRFET3=>PRFET3=43*/\ndydt[PRFET3]=-2.000000e+00*A78*T**n78*np.exp(-1*E78/R/T)*PRFET3*PRFET3-A260*T**n260*np.exp(-1*E260/R/T)*KETM2*PRFET3+A361*T**n361*np.exp(-1*E361/R/T)*PFET3*RCH3O+A356*T**n356*np.exp(-1*E356/R/T)*PFET3*RLIGB\n+A347*T**n347*np.exp(-1*E347/R/T)*PFET3*PRFET3M2+A343*T**n343*np.exp(-1*E343/R/T)*PFET3*RC3H5O2-A22*T**n22*np.exp(-1*E22/R/T)*PRFET3+A357*T**n357*np.exp(-1*E357/R/T)*PFET3*RLIGA+A362*T**n362*np.exp(-1*E362/R/T)*PFET3*RPHENOL-A155*T**n155*np.exp(-1*E155/R/T)*PLIGM2*PRFET3+\nA350*T**n350*np.exp(-1*E350/R/T)*RLIGM2A*PFET3+A354*T**n354*np.exp(-1*E354/R/T)*PFET3*PRADIO+A345*T**n345*np.exp(-1*E345/R/T)*RC3H7O2*PFET3+A359*T**n359*np.exp(-1*E359/R/T)*PFET3*RMGUAI+A358*T**n358*np.exp(-1*E358/R/T)*PRADIOM2*PFET3+A355*T**n355*np.exp(-1*E355/R/T)*RC3H3O*PFET3+\nA346*T**n346*np.exp(-1*E346/R/T)*PFET3*RADIOM2-A323*T**n323*np.exp(-1*E323/R/T)*LIG*PRFET3-A176*T**n176*np.exp(-1*E176/R/T)*PRFET3*LIGM2-A365*T**n365*np.exp(-1*E365/R/T)*ADIO*PRFET3+A348*T**n348*np.exp(-1*E348/R/T)*PFET3*PRLIGH-A218*T**n218*np.exp(-1*E218/R/T)*PRFET3*PFET3M2-\nA134*T**n134*np.exp(-1*E134/R/T)*PLIGH*PRFET3-A239*T**n239*np.exp(-1*E239/R/T)*PRFET3*ADIOM2-A197*T**n197*np.exp(-1*E197/R/T)*PRFET3*LIGM2+A352*T**n352*np.exp(-1*E352/R/T)*PRKETM2*PFET3-A386*T**n386*np.exp(-1*E386/R/T)*KET*PRFET3-A281*T**n281*np.exp(-1*E281/R/T)*C10H2*PRFET3-\nA113*T**n113*np.exp(-1*E113/R/T)*LIGH*PRFET3+A353*T**n353*np.exp(-1*E353/R/T)*PFET3*RKET+A349*T**n349*np.exp(-1*E349/R/T)*RLIGM2B*PFET3+A351*T**n351*np.exp(-1*E351/R/T)*PFET3*RCH3+A360*T**n360*np.exp(-1*E360/R/T)*PFET3*OH-A302*T**n302*np.exp(-1*E302/R/T)*LIG*PRFET3+\nA363*T**n363*np.exp(-1*E363/R/T)*RADIO*PFET3;\n/*ODEforSYNAPYL=>SYNAPYL=44*/\ndydt[SYNAPYL]=-A87*T**n87*np.exp(-1*E87/R/T)*SYNAPYL-A35*T**n35*np.exp(-1*E35/R/T)*SYNAPYL*RPHENOXM2-A39*T**n39*np.exp(-1*E39/R/T)*SYNAPYL*RPHENOX+A23*T**n23*np.exp(-1*E23/R/T)*RADIOM2;\n/*ODEforRKETM2=>RKETM2=45*/\ndydt[RKETM2]=+A260*T**n260*np.exp(-1*E260/R/T)*KETM2*PRFET3+A278*T**n278*np.exp(-1*E278/R/T)*KETM2*RPHENOL+A269*T**n269*np.exp(-1*E269/R/T)*KETM2*RKET+A279*T**n279*np.exp(-1*E279/R/T)*RADIO*KETM2+A265*T**n265*np.exp(-1*E265/R/T)*RLIGM2B*KETM2+\nA270*T**n270*np.exp(-1*E270/R/T)*KETM2*PRADIO+A264*T**n264*np.exp(-1*E264/R/T)*KETM2*PRLIGH+A276*T**n276*np.exp(-1*E276/R/T)*KETM2*OH+A277*T**n277*np.exp(-1*E277/R/T)*KETM2*RCH3O+A259*T**n259*np.exp(-1*E259/R/T)*KETM2*RC3H5O2+A268*T**n268*np.exp(-1*E268/R/T)*KETM2*PRKETM2+\nA274*T**n274*np.exp(-1*E274/R/T)*KETM2*PRADIOM2+A263*T**n263*np.exp(-1*E263/R/T)*KETM2*PRFET3M2-A24*T**n24*np.exp(-1*E24/R/T)*RKETM2-2.000000e+00*A64*T**n64*np.exp(-1*E64/R/T)*RKETM2*RKETM2+A273*T**n273*np.exp(-1*E273/R/T)*KETM2*RLIGA+\nA262*T**n262*np.exp(-1*E262/R/T)*KETM2*RADIOM2+A267*T**n267*np.exp(-1*E267/R/T)*KETM2*RCH3+A272*T**n272*np.exp(-1*E272/R/T)*KETM2*RLIGB+A271*T**n271*np.exp(-1*E271/R/T)*KETM2*RC3H3O+A275*T**n275*np.exp(-1*E275/R/T)*KETM2*RMGUAI+A261*T**n261*np.exp(-1*E261/R/T)*KETM2*RC3H7O2+\nA266*T**n266*np.exp(-1*E266/R/T)*KETM2*RLIGM2A;\n/*ODEforCOUMARYL=>COUMARYL=46*/\ndydt[COUMARYL]=-A89*T**n89*np.exp(-1*E89/R/T)*COUMARYL+A25*T**n25*np.exp(-1*E25/R/T)*RADIO-A47*T**n47*np.exp(-1*E47/R/T)*COUMARYL*RPHENOXM2-A43*T**n43*np.exp(-1*E43/R/T)*COUMARYL*RPHENOX;\n/*ODEforRKET=>RKET=47*/\ndydt[RKET]=+A394*T**n394*np.exp(-1*E394/R/T)*PRKETM2*KET+A392*T**n392*np.exp(-1*E392/R/T)*RLIGM2A*KET+A398*T**n398*np.exp(-1*E398/R/T)*RLIGB*KET+A405*T**n405*np.exp(-1*E405/R/T)*RADIO*KET-A185*T**n185*np.exp(-1*E185/R/T)*RKET*LIGM2-\nA143*T**n143*np.exp(-1*E143/R/T)*PLIGH*RKET-A311*T**n311*np.exp(-1*E311/R/T)*LIG*RKET-A164*T**n164*np.exp(-1*E164/R/T)*PLIGM2*RKET+A400*T**n400*np.exp(-1*E400/R/T)*PRADIOM2*KET-A206*T**n206*np.exp(-1*E206/R/T)*RKET*LIGM2-A227*T**n227*np.exp(-1*E227/R/T)*RKET*PFET3M2+\nA402*T**n402*np.exp(-1*E402/R/T)*KET*OH+A387*T**n387*np.exp(-1*E387/R/T)*RC3H7O2*KET-A269*T**n269*np.exp(-1*E269/R/T)*KETM2*RKET+A389*T**n389*np.exp(-1*E389/R/T)*PRFET3M2*KET-A26*T**n26*np.exp(-1*E26/R/T)*RKET-2.000000e+00*\nA77*T**n77*np.exp(-1*E77/R/T)*RKET*RKET-A290*T**n290*np.exp(-1*E290/R/T)*RKET*C10H2+A397*T**n397*np.exp(-1*E397/R/T)*RC3H3O*KET+A401*T**n401*np.exp(-1*E401/R/T)*RMGUAI*KET-A122*T**n122*np.exp(-1*E122/R/T)*LIGH*RKET+A388*T**n388*np.exp(-1*E388/R/T)*KET*RADIOM2+\nA396*T**n396*np.exp(-1*E396/R/T)*KET*PRADIO-A332*T**n332*np.exp(-1*E332/R/T)*LIG*RKET+A393*T**n393*np.exp(-1*E393/R/T)*RCH3*KET+A390*T**n390*np.exp(-1*E390/R/T)*PRLIGH*KET+A385*T**n385*np.exp(-1*E385/R/T)*KET*RC3H5O2+A386*T**n386*np.exp(-1*E386/R/T)*KET*PRFET3-\nA353*T**n353*np.exp(-1*E353/R/T)*PFET3*RKET-A248*T**n248*np.exp(-1*E248/R/T)*RKET*ADIOM2+A391*T**n391*np.exp(-1*E391/R/T)*RLIGM2B*KET+A404*T**n404*np.exp(-1*E404/R/T)*KET*RPHENOL+A403*T**n403*np.exp(-1*E403/R/T)*RCH3O*KET-A374*T**n374*np.exp(-1*E374/R/T)*ADIO*RKET+\nA399*T**n399*np.exp(-1*E399/R/T)*KET*RLIGA;\n/*ODEforKETD=>KETD=48*/\ndydt[KETD]=+A26*T**n26*np.exp(-1*E26/R/T)*RKET-A46*T**n46*np.exp(-1*E46/R/T)*KETD*RPHENOXM2-A42*T**n42*np.exp(-1*E42/R/T)*KETD*RPHENOX-A92*T**n92*np.exp(-1*E92/R/T)*KETD;\n/*ODEforRC3H7O2=>RC3H7O2=49*/\ndydt[RC3H7O2]=-A27*T**n27*np.exp(-1*E27/R/T)*RC3H7O2-A219*T**n219*np.exp(-1*E219/R/T)*RC3H7O2*PFET3M2-A366*T**n366*np.exp(-1*E366/R/T)*RC3H7O2*ADIO-2.000000e+00*A66*T**n66*np.exp(-1*E66/R/T)*RC3H7O2*RC3H7O2\n-A282*T**n282*np.exp(-1*E282/R/T)*RC3H7O2*C10H2-A135*T**n135*np.exp(-1*E135/R/T)*RC3H7O2*PLIGH+A32*T**n32*np.exp(-1*E32/R/T)*ADIOM2*RPHENOXM2+A36*T**n36*np.exp(-1*E36/R/T)*ADIOM2*RPHENOX+A40*T**n40*np.exp(-1*E40/R/T)*ADIO*RPHENOX-A240*T**n240*np.exp(-1*E240/R/T)*RC3H7O2*ADIOM2-\nA324*T**n324*np.exp(-1*E324/R/T)*LIG*RC3H7O2+A44*T**n44*np.exp(-1*E44/R/T)*ADIO*RPHENOXM2-A114*T**n114*np.exp(-1*E114/R/T)*RC3H7O2*LIGH-A156*T**n156*np.exp(-1*E156/R/T)*PLIGM2*RC3H7O2-A345*T**n345*np.exp(-1*E345/R/T)*RC3H7O2*PFET3-A198*T**n198*np.exp(-1*E198/R/T)*RC3H7O2*LIGM2-\nA303*T**n303*np.exp(-1*E303/R/T)*LIG*RC3H7O2-A177*T**n177*np.exp(-1*E177/R/T)*RC3H7O2*LIGM2-A261*T**n261*np.exp(-1*E261/R/T)*KETM2*RC3H7O2-A387*T**n387*np.exp(-1*E387/R/T)*RC3H7O2*KET;\n/*ODEforRCH3O=>RCH3O=50*/\ndydt[RCH3O]=+A27*T**n27*np.exp(-1*E27/R/T)*RC3H7O2+A49*T**n49*np.exp(-1*E49/R/T)*C10H2M2*RPHENOXM2-A56*T**n56*np.exp(-1*E56/R/T)*RCH3*RCH3O-A361*T**n361*np.exp(-1*E361/R/T)*PFET3*RCH3O-A256*T**n256*np.exp(-1*E256/R/T)*RCH3O*ADIOM2-\nA193*T**n193*np.exp(-1*E193/R/T)*RCH3O*LIGM2-A52*T**n52*np.exp(-1*E52/R/T)*RCH3O*RPHENOX+A48*T**n48*np.exp(-1*E48/R/T)*C10H2M4*RPHENOXM2-A382*T**n382*np.exp(-1*E382/R/T)*RCH3O*ADIO-A172*T**n172*np.exp(-1*E172/R/T)*PLIGM2*RCH3O+A50*T**n50*np.exp(-1*E50/R/T)*C10H2M4*RPHENOX-\n2.000000e+00*A57*T**n57*np.exp(-1*E57/R/T)*RCH3O*RCH3O-A235*T**n235*np.exp(-1*E235/R/T)*RCH3O*PFET3M2+A28*T**n28*np.exp(-1*E28/R/T)*C10H2M4-A277*T**n277*np.exp(-1*E277/R/T)*KETM2*RCH3O-A59*T**n59*np.exp(-1*E59/R/T)*PRKETM2*RCH3O-\nA319*T**n319*np.exp(-1*E319/R/T)*LIG*RCH3O-A298*T**n298*np.exp(-1*E298/R/T)*RCH3O*C10H2-A58*T**n58*np.exp(-1*E58/R/T)*PRADIOM2*RCH3O+A51*T**n51*np.exp(-1*E51/R/T)*C10H2M2*RPHENOX-A53*T**n53*np.exp(-1*E53/R/T)*RCH3O*RPHENOXM2-A340*T**n340*np.exp(-1*E340/R/T)*LIG*RCH3O-\nA151*T**n151*np.exp(-1*E151/R/T)*PLIGH*RCH3O-A214*T**n214*np.exp(-1*E214/R/T)*RCH3O*LIGM2-A130*T**n130*np.exp(-1*E130/R/T)*RCH3O*LIGH-A403*T**n403*np.exp(-1*E403/R/T)*RCH3O*KET+A29*T**n29*np.exp(-1*E29/R/T)*C10H2M2;\n/*ODEforCH3CHO=>CH3CHO=51*/\ndydt[CH3CHO]=+A27*T**n27*np.exp(-1*E27/R/T)*RC3H7O2+A57*T**n57*np.exp(-1*E57/R/T)*RCH3O*RCH3O;\n/*ODEforC10H2M2=>C10H2M2=52*/\ndydt[C10H2M2]=-5.000000e-01*A49*T**n49*np.exp(-1*E49/R/T)*C10H2M2*RPHENOXM2+5.000000e-01*A50*T**n50*np.exp(-1*E50/R/T)*C10H2M4*RPHENOX-\n5.000000e-01*A51*T**n51*np.exp(-1*E51/R/T)*C10H2M2*RPHENOX+5.000000e-01*A48*T**n48*np.exp(-1*E48/R/T)*C10H2M4*RPHENOXM2-A29*T**n29*np.exp(-1*E29/R/T)*C10H2M2;\n/*ODEforCH2CO=>CH2CO=53*/\ndydt[CH2CO]=+A108*T**n108*np.exp(-1*E108/R/T)*LIGC+A30*T**n30*np.exp(-1*E30/R/T)*PLIGC;\n/*ODEforCO2=>CO2=54*/\ndydt[CO2]=+A107*T**n107*np.exp(-1*E107/R/T)*LIGO+A53*T**n53*np.exp(-1*E53/R/T)*RCH3O*RPHENOXM2+A31*T**n31*np.exp(-1*E31/R/T)*PLIGO+A52*T**n52*np.exp(-1*E52/R/T)*RCH3O*RPHENOX;\n/*ODEforADIOM2=>ADIOM2=55*/\ndydt[ADIOM2]=+A367*T**n367*np.exp(-1*E367/R/T)*ADIO*RADIOM2+A283*T**n283*np.exp(-1*E283/R/T)*C10H2*RADIOM2-A256*T**n256*np.exp(-1*E256/R/T)*RCH3O*ADIOM2-A84*T**n84*np.exp(-1*E84/R/T)*ADIOM2-A32*T**n32*np.exp(-1*E32/R/T)*ADIOM2*RPHENOXM2+\nA115*T**n115*np.exp(-1*E115/R/T)*LIGH*RADIOM2-A249*T**n249*np.exp(-1*E249/R/T)*PRADIO*ADIOM2+A199*T**n199*np.exp(-1*E199/R/T)*LIGM2*RADIOM2-A243*T**n243*np.exp(-1*E243/R/T)*PRLIGH*ADIOM2+A304*T**n304*np.exp(-1*E304/R/T)*LIG*RADIOM2-A244*T**n244*np.exp(-1*E244/R/T)*RLIGM2B*ADIOM2-\nA257*T**n257*np.exp(-1*E257/R/T)*ADIOM2*RPHENOL+A262*T**n262*np.exp(-1*E262/R/T)*KETM2*RADIOM2+A220*T**n220*np.exp(-1*E220/R/T)*PFET3M2*RADIOM2-A242*T**n242*np.exp(-1*E242/R/T)*PRFET3M2*ADIOM2-A238*T**n238*np.exp(-1*E238/R/T)*RC3H5O2*ADIOM2-A255*T**n255*np.exp(-1*E255/R/T)*OH*ADIOM2-\nA247*T**n247*np.exp(-1*E247/R/T)*PRKETM2*ADIOM2+A346*T**n346*np.exp(-1*E346/R/T)*PFET3*RADIOM2-A245*T**n245*np.exp(-1*E245/R/T)*RLIGM2A*ADIOM2+A178*T**n178*np.exp(-1*E178/R/T)*LIGM2*RADIOM2-A246*T**n246*np.exp(-1*E246/R/T)*RCH3*ADIOM2+A136*T**n136*np.exp(-1*E136/R/T)*PLIGH*RADIOM2-\nA36*T**n36*np.exp(-1*E36/R/T)*ADIOM2*RPHENOX-A253*T**n253*np.exp(-1*E253/R/T)*PRADIOM2*ADIOM2+A388*T**n388*np.exp(-1*E388/R/T)*KET*RADIOM2+A325*T**n325*np.exp(-1*E325/R/T)*LIG*RADIOM2-A240*T**n240*np.exp(-1*E240/R/T)*RC3H7O2*ADIOM2+A157*T**n157*np.exp(-1*E157/R/T)*PLIGM2*RADIOM2-\nA239*T**n239*np.exp(-1*E239/R/T)*PRFET3*ADIOM2-A252*T**n252*np.exp(-1*E252/R/T)*RLIGA*ADIOM2-A254*T**n254*np.exp(-1*E254/R/T)*RMGUAI*ADIOM2-A250*T**n250*np.exp(-1*E250/R/T)*RC3H3O*ADIOM2-A248*T**n248*np.exp(-1*E248/R/T)*RKET*ADIOM2-A258*T**n258*np.exp(-1*E258/R/T)*RADIO*ADIOM2-\nA251*T**n251*np.exp(-1*E251/R/T)*RLIGB*ADIOM2;\n/*ODEforPH2=>PH2=56*/\ndydt[PH2]=+2.000000e+00*A78*T**n78*np.exp(-1*E78/R/T)*PRFET3*PRFET3+3.000000e+00*A74*T**n74*np.exp(-1*E74/R/T)*RADIO*RADIO+\n2.000000e+00*A60*T**n60*np.exp(-1*E60/R/T)*RADIOM2*RADIOM2+A68*T**n68*np.exp(-1*E68/R/T)*RC3H3O*RC3H3O+2.000000e+00*A63*T**n63*np.exp(-1*E63/R/T)*RMGUAI*RMGUAI+A45*T**n45*np.exp(-1*E45/R/T)*KET*RPHENOXM2+\nA32*T**n32*np.exp(-1*E32/R/T)*ADIOM2*RPHENOXM2-A96*T**n96*np.exp(-1*E96/R/T)*PH2+2.000000e+00*A41*T**n41*np.exp(-1*E41/R/T)*KET*RPHENOX+2.000000e+00*A40*T**n40*np.exp(-1*E40/R/T)*ADIO*RPHENOX+\n2.000000e+00*A80*T**n80*np.exp(-1*E80/R/T)*RPHENOX*RPHENOL+5.000000e+00*A75*T**n75*np.exp(-1*E75/R/T)*RLIGB*RLIGB+A44*T**n44*np.exp(-1*E44/R/T)*ADIO*RPHENOXM2+\n2.000000e+00*A43*T**n43*np.exp(-1*E43/R/T)*COUMARYL*RPHENOX+A33*T**n33*np.exp(-1*E33/R/T)*KETM2*RPHENOXM2+A65*T**n65*np.exp(-1*E65/R/T)*PRFET3M2*PRFET3M2+2.000000e+00*A64*T**n64*np.exp(-1*E64/R/T)*RKETM2*RKETM2+\n3.000000e+00*A61*T**n61*np.exp(-1*E61/R/T)*RLIGM2B*RLIGM2B+2.000000e+00*A47*T**n47*np.exp(-1*E47/R/T)*COUMARYL*RPHENOXM2+5.000000e-01*\nA54*T**n54*np.exp(-1*E54/R/T)*RCH3*RPHENOXM2+A110*T**n110*np.exp(-1*E110/R/T)*PRADIOM2*PRADIOM2+3.000000e+00*A62*T**n62*np.exp(-1*E62/R/T)*RLIGM2A*RLIGM2A+2.000000e+00*A38*T**n38*np.exp(-1*E38/R/T)*RPHENOX*KETDM2\n+3.000000e+00*A73*T**n73*np.exp(-1*E73/R/T)*RLIGB*RPHENOX+2.000000e+00*A39*T**n39*np.exp(-1*E39/R/T)*SYNAPYL*RPHENOX+2.000000e+00*\nA77*T**n77*np.exp(-1*E77/R/T)*RKET*RKET+A66*T**n66*np.exp(-1*E66/R/T)*RC3H7O2*RC3H7O2+2.000000e+00*A42*T**n42*np.exp(-1*E42/R/T)*KETD*RPHENOX+A81*T**n81*np.exp(-1*E81/R/T)*RC3H3O*RPHENOX+A36*T**n36*np.exp(-1*E36/R/T)*ADIOM2*RPHENOX+\n2.000000e+00*A37*T**n37*np.exp(-1*E37/R/T)*KETM2*RPHENOX+A46*T**n46*np.exp(-1*E46/R/T)*KETD*RPHENOXM2+1.500000e+00*A35*T**n35*np.exp(-1*E35/R/T)*SYNAPYL*RPHENOXM2+\n2.000000e+00*A109*T**n109*np.exp(-1*E109/R/T)*PRADIO*PRADIO+A34*T**n34*np.exp(-1*E34/R/T)*RPHENOXM2*KETDM2+5.000000e+00*A76*T**n76*np.exp(-1*E76/R/T)*RLIGA*RLIGA+\n1.500000e+00*A55*T**n55*np.exp(-1*E55/R/T)*RCH3*RPHENOX+A67*T**n67*np.exp(-1*E67/R/T)*RC3H5O2*RC3H5O2;\n/*ODEforPCOH=>PCOH=57*/\ndydt[PCOH]=+A74*T**n74*np.exp(-1*E74/R/T)*RADIO*RADIO+2.000000e+00*A60*T**n60*np.exp(-1*E60/R/T)*RADIOM2*RADIOM2+A49*T**n49*np.exp(-1*E49/R/T)*C10H2M2*RPHENOXM2+3.000000e+00\n*A77*T**n77*np.exp(-1*E77/R/T)*RKET*RKET+A39*T**n39*np.exp(-1*E39/R/T)*SYNAPYL*RPHENOX+2.000000e+00*A68*T**n68*np.exp(-1*E68/R/T)*RC3H3O*RC3H3O+2.000000e+00*A63*T**n63*np.exp(-1*E63/R/T)*RMGUAI*RMGUAI\n+A48*T**n48*np.exp(-1*E48/R/T)*C10H2M4*RPHENOXM2+A45*T**n45*np.exp(-1*E45/R/T)*KET*RPHENOXM2+2.000000e+00*A32*T**n32*np.exp(-1*E32/R/T)*ADIOM2*RPHENOXM2+A81*T**n81*np.exp(-1*E81/R/T)*RC3H3O*RPHENOX+A36*T**n36*np.exp(-1*E36/R/T)*ADIOM2*RPHENOX-\nA95*T**n95*np.exp(-1*E95/R/T)*PCOH+A46*T**n46*np.exp(-1*E46/R/T)*KETD*RPHENOXM2+2.000000e+00*A33*T**n33*np.exp(-1*E33/R/T)*KETM2*RPHENOXM2+A44*T**n44*np.exp(-1*E44/R/T)*ADIO*RPHENOXM2+2.000000e+00*\nA35*T**n35*np.exp(-1*E35/R/T)*SYNAPYL*RPHENOXM2+A109*T**n109*np.exp(-1*E109/R/T)*PRADIO*PRADIO+2.000000e+00*A65*T**n65*np.exp(-1*E65/R/T)*PRFET3M2*PRFET3M2+2.000000e+00*A34*T**n34*np.exp(-1*E34/R/T)*RPHENOXM2*KETDM2\n+4.000000e+00*A64*T**n64*np.exp(-1*E64/R/T)*RKETM2*RKETM2+4.000000e+00*A61*T**n61*np.exp(-1*E61/R/T)*RLIGM2B*RLIGM2B+A47*T**n47*np.exp(-1*E47/R/T)*COUMARYL*RPHENOXM2+\n4.000000e+00*A79*T**n79*np.exp(-1*E79/R/T)*RLIGH*RLIGH+A54*T**n54*np.exp(-1*E54/R/T)*RCH3*RPHENOXM2+2.000000e+00*A110*T**n110*np.exp(-1*E110/R/T)*PRADIOM2*PRADIOM2+\n4.000000e+00*A62*T**n62*np.exp(-1*E62/R/T)*RLIGM2A*RLIGM2A+A38*T**n38*np.exp(-1*E38/R/T)*RPHENOX*KETDM2+A55*T**n55*np.exp(-1*E55/R/T)*RCH3*RPHENOX+2.000000e+00*A67*T**n67*np.exp(-1*E67/R/T)*RC3H5O2*RC3H5O2;\n/*ODEforRC3H5O2=>RC3H5O2=58*/\ndydt[RC3H5O2]=-A301*T**n301*np.exp(-1*E301/R/T)*LIG*RC3H5O2-A133*T**n133*np.exp(-1*E133/R/T)*PLIGH*RC3H5O2-A364*T**n364*np.exp(-1*E364/R/T)*ADIO*RC3H5O2+A45*T**n45*np.exp(-1*E45/R/T)*KET*RPHENOXM2-A343*T**n343*np.exp(-1*E343/R/T)*PFET3*RC3H5O2+\nA41*T**n41*np.exp(-1*E41/R/T)*KET*RPHENOX-A280*T**n280*np.exp(-1*E280/R/T)*C10H2*RC3H5O2+A37*T**n37*np.exp(-1*E37/R/T)*KETM2*RPHENOX+A33*T**n33*np.exp(-1*E33/R/T)*KETM2*RPHENOXM2-A217*T**n217*np.exp(-1*E217/R/T)*RC3H5O2*PFET3M2-A259*T**n259*np.exp(-1*E259/R/T)*KETM2*RC3H5O2-\nA175*T**n175*np.exp(-1*E175/R/T)*RC3H5O2*LIGM2-A385*T**n385*np.exp(-1*E385/R/T)*KET*RC3H5O2-A112*T**n112*np.exp(-1*E112/R/T)*LIGH*RC3H5O2-A196*T**n196*np.exp(-1*E196/R/T)*RC3H5O2*LIGM2-A322*T**n322*np.exp(-1*E322/R/T)*LIG*RC3H5O2-A154*T**n154*np.exp(-1*E154/R/T)*PLIGM2*RC3H5O2-\nA238*T**n238*np.exp(-1*E238/R/T)*RC3H5O2*ADIOM2-2.000000e+00*A67*T**n67*np.exp(-1*E67/R/T)*RC3H5O2*RC3H5O2;\n/*ODEforRC3H3O=>RC3H3O=59*/\ndydt[RC3H3O]=-A229*T**n229*np.exp(-1*E229/R/T)*RC3H3O*PFET3M2-A292*T**n292*np.exp(-1*E292/R/T)*RC3H3O*C10H2+A39*T**n39*np.exp(-1*E39/R/T)*SYNAPYL*RPHENOX-2.000000e+00*A68*T**n68*np.exp(-1*E68/R/T)*RC3H3O*RC3H3O\n-A397*T**n397*np.exp(-1*E397/R/T)*RC3H3O*KET+A42*T**n42*np.exp(-1*E42/R/T)*KETD*RPHENOX-A187*T**n187*np.exp(-1*E187/R/T)*RC3H3O*LIGM2-A81*T**n81*np.exp(-1*E81/R/T)*RC3H3O*RPHENOX+A46*T**n46*np.exp(-1*E46/R/T)*KETD*RPHENOXM2+A43*T**n43*np.exp(-1*E43/R/T)*COUMARYL*RPHENOX+\nA35*T**n35*np.exp(-1*E35/R/T)*SYNAPYL*RPHENOXM2-A208*T**n208*np.exp(-1*E208/R/T)*RC3H3O*LIGM2-A124*T**n124*np.exp(-1*E124/R/T)*RC3H3O*LIGH+A34*T**n34*np.exp(-1*E34/R/T)*RPHENOXM2*KETDM2-A145*T**n145*np.exp(-1*E145/R/T)*RC3H3O*PLIGH+A47*T**n47*np.exp(-1*E47/R/T)*COUMARYL*RPHENOXM2-\nA250*T**n250*np.exp(-1*E250/R/T)*RC3H3O*ADIOM2-A334*T**n334*np.exp(-1*E334/R/T)*LIG*RC3H3O-A313*T**n313*np.exp(-1*E313/R/T)*LIG*RC3H3O-A376*T**n376*np.exp(-1*E376/R/T)*RC3H3O*ADIO-A271*T**n271*np.exp(-1*E271/R/T)*KETM2*RC3H3O-A166*T**n166*np.exp(-1*E166/R/T)*PLIGM2*RC3H3O+\nA38*T**n38*np.exp(-1*E38/R/T)*RPHENOX*KETDM2-A355*T**n355*np.exp(-1*E355/R/T)*RC3H3O*PFET3;\n/*ODEforPCOS=>PCOS=60*/\ndydt[PCOS]=+2.000000e+00*A78*T**n78*np.exp(-1*E78/R/T)*PRFET3*PRFET3+A74*T**n74*np.exp(-1*E74/R/T)*RADIO*RADIO+A77*T**n77*np.exp(-1*E77/R/T)*RKET*RKET+A39*T**n39*np.exp(-1*E39/R/T)*SYNAPYL*RPHENOX+\n2.000000e+00*A42*T**n42*np.exp(-1*E42/R/T)*KETD*RPHENOX+A45*T**n45*np.exp(-1*E45/R/T)*KET*RPHENOXM2+A83*T**n83*np.exp(-1*E83/R/T)*CHAR*RPHENOXM2+A81*T**n81*np.exp(-1*E81/R/T)*RC3H3O*RPHENOX+A36*T**n36*np.exp(-1*E36/R/T)*ADIOM2*RPHENOX+\n2.000000e+00*A41*T**n41*np.exp(-1*E41/R/T)*KET*RPHENOX+A50*T**n50*np.exp(-1*E50/R/T)*C10H2M4*RPHENOX+2.000000e+00*A40*T**n40*np.exp(-1*E40/R/T)*ADIO*RPHENOX+\n2.000000e+00*A80*T**n80*np.exp(-1*E80/R/T)*RPHENOX*RPHENOL+2.000000e+00*A37*T**n37*np.exp(-1*E37/R/T)*KETM2*RPHENOX-A94*T**n94*np.exp(-1*E94/R/T)*PCOS+A46*T**n46*np.exp(-1*E46/R/T)*KETD*RPHENOXM2+\n4.000000e+00*A75*T**n75*np.exp(-1*E75/R/T)*RLIGB*RLIGB+2.000000e+00*A43*T**n43*np.exp(-1*E43/R/T)*COUMARYL*RPHENOX+A44*T**n44*np.exp(-1*E44/R/T)*ADIO*RPHENOXM2+A109*T**n109*np.exp(-1*E109/R/T)*PRADIO*PRADIO\n+A47*T**n47*np.exp(-1*E47/R/T)*COUMARYL*RPHENOXM2+A51*T**n51*np.exp(-1*E51/R/T)*C10H2M2*RPHENOX+4.000000e+00*A76*T**n76*np.exp(-1*E76/R/T)*RLIGA*RLIGA+A38*T**n38*np.exp(-1*E38/R/T)*RPHENOX*KETDM2+A82*T**n82*np.exp(-1*E82/R/T)*CHAR*RPHENOX+\n3.000000e+00*A73*T**n73*np.exp(-1*E73/R/T)*RLIGB*RPHENOX;\n/*ODEforADIO=>ADIO=61*/\ndydt[ADIO]=-A367*T**n367*np.exp(-1*E367/R/T)*ADIO*RADIOM2+A216*T**n216*np.exp(-1*E216/R/T)*RADIO*LIGM2+A153*T**n153*np.exp(-1*E153/R/T)*RADIO*PLIGH+A132*T**n132*np.exp(-1*E132/R/T)*RADIO*LIGH+A279*T**n279*np.exp(-1*E279/R/T)*RADIO*KETM2-\nA372*T**n372*np.exp(-1*E372/R/T)*RCH3*ADIO-A364*T**n364*np.exp(-1*E364/R/T)*ADIO*RC3H5O2-A377*T**n377*np.exp(-1*E377/R/T)*RLIGB*ADIO-A40*T**n40*np.exp(-1*E40/R/T)*ADIO*RPHENOX-A373*T**n373*np.exp(-1*E373/R/T)*PRKETM2*ADIO-A44*T**n44*np.exp(-1*E44/R/T)*ADIO*RPHENOXM2+\nA405*T**n405*np.exp(-1*E405/R/T)*RADIO*KET-A381*T**n381*np.exp(-1*E381/R/T)*ADIO*OH+A174*T**n174*np.exp(-1*E174/R/T)*RADIO*PLIGM2-A371*T**n371*np.exp(-1*E371/R/T)*RLIGM2A*ADIO-A380*T**n380*np.exp(-1*E380/R/T)*RMGUAI*ADIO-A90*T**n90*np.exp(-1*E90/R/T)*ADIO+\nA321*T**n321*np.exp(-1*E321/R/T)*RADIO*LIG+A342*T**n342*np.exp(-1*E342/R/T)*RADIO*LIG-A379*T**n379*np.exp(-1*E379/R/T)*PRADIOM2*ADIO+A195*T**n195*np.exp(-1*E195/R/T)*RADIO*LIGM2-A365*T**n365*np.exp(-1*E365/R/T)*ADIO*PRFET3-A366*T**n366*np.exp(-1*E366/R/T)*RC3H7O2*ADIO-\nA383*T**n383*np.exp(-1*E383/R/T)*ADIO*RPHENOL-A382*T**n382*np.exp(-1*E382/R/T)*RCH3O*ADIO-A378*T**n378*np.exp(-1*E378/R/T)*ADIO*RLIGA-A370*T**n370*np.exp(-1*E370/R/T)*RLIGM2B*ADIO+A237*T**n237*np.exp(-1*E237/R/T)*RADIO*PFET3M2-A375*T**n375*np.exp(-1*E375/R/T)*ADIO*PRADIO-\nA368*T**n368*np.exp(-1*E368/R/T)*PRFET3M2*ADIO-A376*T**n376*np.exp(-1*E376/R/T)*RC3H3O*ADIO+A258*T**n258*np.exp(-1*E258/R/T)*RADIO*ADIOM2+A363*T**n363*np.exp(-1*E363/R/T)*RADIO*PFET3+A300*T**n300*np.exp(-1*E300/R/T)*RADIO*C10H2-A369*T**n369*np.exp(-1*E369/R/T)*PRLIGH*ADIO-\nA374*T**n374*np.exp(-1*E374/R/T)*ADIO*RKET;\n/*ODEforRCH3=>RCH3=62*/\ndydt[RCH3]=-A309*T**n309*np.exp(-1*E309/R/T)*LIG*RCH3-A56*T**n56*np.exp(-1*E56/R/T)*RCH3*RCH3O-A246*T**n246*np.exp(-1*E246/R/T)*RCH3*ADIOM2-2.000000e+00*A70*T**n70*np.exp(-1*E70/R/T)*RCH3*RCH3-\nA372*T**n372*np.exp(-1*E372/R/T)*RCH3*ADIO-A330*T**n330*np.exp(-1*E330/R/T)*LIG*RCH3+A52*T**n52*np.exp(-1*E52/R/T)*RCH3O*RPHENOX-A204*T**n204*np.exp(-1*E204/R/T)*RCH3*LIGM2-A120*T**n120*np.exp(-1*E120/R/T)*RCH3*LIGH-A141*T**n141*np.exp(-1*E141/R/T)*PLIGH*RCH3-\nA71*T**n71*np.exp(-1*E71/R/T)*PRADIOM2*RCH3-A393*T**n393*np.exp(-1*E393/R/T)*RCH3*KET-A69*T**n69*np.exp(-1*E69/R/T)*RCH3*OH-A183*T**n183*np.exp(-1*E183/R/T)*RCH3*LIGM2-A72*T**n72*np.exp(-1*E72/R/T)*PRKETM2*RCH3-A54*T**n54*np.exp(-1*E54/R/T)*RCH3*RPHENOXM2-\nA351*T**n351*np.exp(-1*E351/R/T)*PFET3*RCH3+A53*T**n53*np.exp(-1*E53/R/T)*RCH3O*RPHENOXM2-A267*T**n267*np.exp(-1*E267/R/T)*KETM2*RCH3-A162*T**n162*np.exp(-1*E162/R/T)*PLIGM2*RCH3-A225*T**n225*np.exp(-1*E225/R/T)*RCH3*PFET3M2-A288*T**n288*np.exp(-1*E288/R/T)*RCH3*C10H2-\nA55*T**n55*np.exp(-1*E55/R/T)*RCH3*RPHENOX+A101*T**n101*np.exp(-1*E101/R/T)*PCH3;\n/*ODEforPCH3=>PCH3=63*/\ndydt[PCH3]=+A72*T**n72*np.exp(-1*E72/R/T)*PRKETM2*RCH3+A71*T**n71*np.exp(-1*E71/R/T)*PRADIOM2*RCH3+A55*T**n55*np.exp(-1*E55/R/T)*RCH3*RPHENOX-A101*T**n101*np.exp(-1*E101/R/T)*PCH3+2.000000e+00*\nA79*T**n79*np.exp(-1*E79/R/T)*RLIGH*RLIGH+A54*T**n54*np.exp(-1*E54/R/T)*RCH3*RPHENOXM2;\n/*ODEforETOH=>ETOH=64*/\ndydt[ETOH]=+A56*T**n56*np.exp(-1*E56/R/T)*RCH3*RCH3O;\n/*ODEforH2O=>H2O=65*/\ndydt[H2O]=+A234*T**n234*np.exp(-1*E234/R/T)*OH*PFET3M2+2.000000e+00*A66*T**n66*np.exp(-1*E66/R/T)*RC3H7O2*RC3H7O2+A339*T**n339*np.exp(-1*E339/R/T)*LIG*OH+A129*T**n129*np.exp(-1*E129/R/T)*LIGH*OH+\nA192*T**n192*np.exp(-1*E192/R/T)*LIGM2*OH+A297*T**n297*np.exp(-1*E297/R/T)*C10H2*OH+A318*T**n318*np.exp(-1*E318/R/T)*LIG*OH+A57*T**n57*np.exp(-1*E57/R/T)*RCH3O*RCH3O+A276*T**n276*np.exp(-1*E276/R/T)*KETM2*OH+A171*T**n171*np.exp(-1*E171/R/T)*PLIGM2*OH+\nA213*T**n213*np.exp(-1*E213/R/T)*LIGM2*OH+A381*T**n381*np.exp(-1*E381/R/T)*ADIO*OH+4.000000e+00*A79*T**n79*np.exp(-1*E79/R/T)*RLIGH*RLIGH+2.000000e+00*A62*T**n62*np.exp(-1*E62/R/T)*RLIGM2A*RLIGM2A\n+A360*T**n360*np.exp(-1*E360/R/T)*PFET3*OH+2.000000e+00*A76*T**n76*np.exp(-1*E76/R/T)*RLIGA*RLIGA+A402*T**n402*np.exp(-1*E402/R/T)*KET*OH+A73*T**n73*np.exp(-1*E73/R/T)*RLIGB*RPHENOX+A255*T**n255*np.exp(-1*E255/R/T)*OH*ADIOM2\n+2.000000e+00*A67*T**n67*np.exp(-1*E67/R/T)*RC3H5O2*RC3H5O2+A150*T**n150*np.exp(-1*E150/R/T)*PLIGH*OH;\n/*ODEforPCH2OH=>PCH2OH=66*/\ndydt[PCH2OH]=+2.000000e+00*A74*T**n74*np.exp(-1*E74/R/T)*RADIO*RADIO+2.000000e+00*A60*T**n60*np.exp(-1*E60/R/T)*RADIOM2*RADIOM2-\nA98*T**n98*np.exp(-1*E98/R/T)*PCH2OH+2.000000e+00*A109*T**n109*np.exp(-1*E109/R/T)*PRADIO*PRADIO+2.000000e+00*A77*T**n77*np.exp(-1*E77/R/T)*RKET*RKET+A58*T**n58*np.exp(-1*E58/R/T)*PRADIOM2*RCH3O\n+2.000000e+00*A64*T**n64*np.exp(-1*E64/R/T)*RKETM2*RKETM2+2.000000e+00*A61*T**n61*np.exp(-1*E61/R/T)*RLIGM2B*RLIGM2B+2.000000e+00*\nA110*T**n110*np.exp(-1*E110/R/T)*PRADIOM2*PRADIOM2+A59*T**n59*np.exp(-1*E59/R/T)*PRKETM2*RCH3O+2.000000e+00*A75*T**n75*np.exp(-1*E75/R/T)*RLIGB*RLIGB;\n/*ODEforPCHP2=>PCHP2=67*/\ndydt[PCHP2]=+2.000000e+00*A78*T**n78*np.exp(-1*E78/R/T)*PRFET3*PRFET3+2.000000e+00*A74*T**n74*np.exp(-1*E74/R/T)*RADIO*RADIO+\n2.000000e+00*A60*T**n60*np.exp(-1*E60/R/T)*RADIOM2*RADIOM2+2.000000e+00*A109*T**n109*np.exp(-1*E109/R/T)*PRADIO*PRADIO-A103*T**n103*np.exp(-1*E103/R/T)*PCHP2+\n2.000000e+00*A65*T**n65*np.exp(-1*E65/R/T)*PRFET3M2*PRFET3M2+2.000000e+00*A77*T**n77*np.exp(-1*E77/R/T)*RKET*RKET+2.000000e+00*\nA64*T**n64*np.exp(-1*E64/R/T)*RKETM2*RKETM2+A102*T**n102*np.exp(-1*E102/R/T)*PCH2P+2.000000e+00*A79*T**n79*np.exp(-1*E79/R/T)*RLIGH*RLIGH+2.000000e+00*\nA110*T**n110*np.exp(-1*E110/R/T)*PRADIOM2*PRADIOM2;\n/*ODEforPCHOHP=>PCHOHP=68*/\ndydt[PCHOHP]=-A99*T**n99*np.exp(-1*E99/R/T)*PCHOHP+2.000000e+00*A78*T**n78*np.exp(-1*E78/R/T)*PRFET3*PRFET3+2.000000e+00*\nA74*T**n74*np.exp(-1*E74/R/T)*RADIO*RADIO+2.000000e+00*A60*T**n60*np.exp(-1*E60/R/T)*RADIOM2*RADIOM2+2.000000e+00*A109*T**n109*np.exp(-1*E109/R/T)*PRADIO*PRADIO+\n2.000000e+00*A65*T**n65*np.exp(-1*E65/R/T)*PRFET3M2*PRFET3M2+2.000000e+00*A66*T**n66*np.exp(-1*E66/R/T)*RC3H7O2*RC3H7O2+2.000000e+00*\nA61*T**n61*np.exp(-1*E61/R/T)*RLIGM2B*RLIGM2B+2.000000e+00*A110*T**n110*np.exp(-1*E110/R/T)*PRADIOM2*PRADIOM2+2.000000e+00*A75*T**n75*np.exp(-1*E75/R/T)*RLIGB*RLIGB;\n/*ODEforCHAR=>CHAR=69*/\ndydt[CHAR]=+5.000000e-01*A284*T**n284*np.exp(-1*E284/R/T)*PRFET3M2*C10H2+5.000000e-01*A292*T**n292*np.exp(-1*E292/R/T)*RC3H3O*C10H2+\n1.000000e-01*A103*T**n103*np.exp(-1*E103/R/T)*PCHP2+2.000000e-01*A97*T**n97*np.exp(-1*E97/R/T)*PC2H2+5.000000e-01*\nA290*T**n290*np.exp(-1*E290/R/T)*RKET*C10H2+5.000000e-01*A283*T**n283*np.exp(-1*E283/R/T)*C10H2*RADIOM2+5.000000e-01*A282*T**n282*np.exp(-1*E282/R/T)*RC3H7O2*C10H2+\n5.000000e-01*A287*T**n287*np.exp(-1*E287/R/T)*RLIGM2A*C10H2+5.000000e-01*A297*T**n297*np.exp(-1*E297/R/T)*C10H2*OH+1.000000e-01*\nA100*T**n100*np.exp(-1*E100/R/T)*PCOHP2+5.000000e-01*A289*T**n289*np.exp(-1*E289/R/T)*PRKETM2*C10H2+5.000000e-01*A286*T**n286*np.exp(-1*E286/R/T)*RLIGM2B*C10H2+\n5.000000e-01*A280*T**n280*np.exp(-1*E280/R/T)*C10H2*RC3H5O2+5.000000e-01*A299*T**n299*np.exp(-1*E299/R/T)*C10H2*RPHENOL+5.000000e-01*\nA295*T**n295*np.exp(-1*E295/R/T)*PRADIOM2*C10H2+5.000000e-01*A294*T**n294*np.exp(-1*E294/R/T)*RLIGA*C10H2+2.000000e-01*A75*T**n75*np.exp(-1*E75/R/T)*RLIGB*RLIGB+\n5.000000e-01*A296*T**n296*np.exp(-1*E296/R/T)*RMGUAI*C10H2+5.000000e-01*A298*T**n298*np.exp(-1*E298/R/T)*RCH3O*C10H2+5.000000e-01*\nA285*T**n285*np.exp(-1*E285/R/T)*PRLIGH*C10H2+5.000000e-01*A281*T**n281*np.exp(-1*E281/R/T)*C10H2*PRFET3+2.000000e-01*A61*T**n61*np.exp(-1*E61/R/T)*RLIGM2B*RLIGM2B+\n5.000000e-01*A293*T**n293*np.exp(-1*E293/R/T)*RLIGB*C10H2+5.000000e-01*A288*T**n288*np.exp(-1*E288/R/T)*RCH3*C10H2+5.000000e-01*\nA291*T**n291*np.exp(-1*E291/R/T)*C10H2*PRADIO+5.000000e-01*A300*T**n300*np.exp(-1*E300/R/T)*RADIO*C10H2;\n/*ODEforPCOHP2=>PCOHP2=70*/\ndydt[PCOHP2]=-A100*T**n100*np.exp(-1*E100/R/T)*PCOHP2+2.000000e+00*A62*T**n62*np.exp(-1*E62/R/T)*RLIGM2A*RLIGM2A+2.000000e+00*\nA76*T**n76*np.exp(-1*E76/R/T)*RLIGA*RLIGA+A73*T**n73*np.exp(-1*E73/R/T)*RLIGB*RPHENOX+2.000000e+00*A79*T**n79*np.exp(-1*E79/R/T)*RLIGH*RLIGH;\n/*ODEforPC2H2=>PC2H2=71*/\ndydt[PC2H2]=+5.000000e-01*A99*T**n99*np.exp(-1*E99/R/T)*PCHOHP+A81*T**n81*np.exp(-1*E81/R/T)*RC3H3O*RPHENOX+2.000000e+00*A62*T**n62*np.exp(-1*E62/R/T)*RLIGM2A*RLIGM2A\n-A97*T**n97*np.exp(-1*E97/R/T)*PC2H2+2.000000e+00*A76*T**n76*np.exp(-1*E76/R/T)*RLIGA*RLIGA+2.000000e+00*A68*T**n68*np.exp(-1*E68/R/T)*RC3H3O*RC3H3O+\n2.000000e+00*A66*T**n66*np.exp(-1*E66/R/T)*RC3H7O2*RC3H7O2+A73*T**n73*np.exp(-1*E73/R/T)*RLIGB*RPHENOX+2.000000e+00*A79*T**n79*np.exp(-1*E79/R/T)*RLIGH*RLIGH+\n2.000000e+00*A67*T**n67*np.exp(-1*E67/R/T)*RC3H5O2*RC3H5O2;\n/*ODEforPCHO=>PCHO=72*/\ndydt[PCHO]=+2.000000e+00*A78*T**n78*np.exp(-1*E78/R/T)*PRFET3*PRFET3+2.000000e+00*A65*T**n65*np.exp(-1*E65/R/T)*PRFET3M2*PRFET3M2-\nA111*T**n111*np.exp(-1*E111/R/T)*PCHO;\n/*ODEforCH3OH=>CH3OH=73*/\ndydt[CH3OH]=+A319*T**n319*np.exp(-1*E319/R/T)*LIG*RCH3O+A69*T**n69*np.exp(-1*E69/R/T)*RCH3*OH+A298*T**n298*np.exp(-1*E298/R/T)*RCH3O*C10H2+A361*T**n361*np.exp(-1*E361/R/T)*PFET3*RCH3O+A256*T**n256*np.exp(-1*E256/R/T)*RCH3O*ADIOM2+\nA193*T**n193*np.exp(-1*E193/R/T)*RCH3O*LIGM2+A382*T**n382*np.exp(-1*E382/R/T)*RCH3O*ADIO+A172*T**n172*np.exp(-1*E172/R/T)*PLIGM2*RCH3O+A340*T**n340*np.exp(-1*E340/R/T)*LIG*RCH3O+A151*T**n151*np.exp(-1*E151/R/T)*PLIGH*RCH3O+A235*T**n235*np.exp(-1*E235/R/T)*RCH3O*PFET3M2+\nA130*T**n130*np.exp(-1*E130/R/T)*RCH3O*LIGH+A214*T**n214*np.exp(-1*E214/R/T)*RCH3O*LIGM2+A403*T**n403*np.exp(-1*E403/R/T)*RCH3O*KET+A277*T**n277*np.exp(-1*E277/R/T)*KETM2*RCH3O;\n/*ODEforC2H6=>C2H6=74*/\ndydt[C2H6]=+A70*T**n70*np.exp(-1*E70/R/T)*RCH3*RCH3;\n/*ODEforPCH2P=>PCH2P=75*/\ndydt[PCH2P]=+A98*T**n98*np.exp(-1*E98/R/T)*PCH2OH-A102*T**n102*np.exp(-1*E102/R/T)*PCH2P+2.000000e+00*A79*T**n79*np.exp(-1*E79/R/T)*RLIGH*RLIGH;\n/*ODEforVADIOM2=>VADIOM2=76*/\ndydt[VADIOM2]=+A84*T**n84*np.exp(-1*E84/R/T)*ADIOM2;\n/*ODEforVKETM2=>VKETM2=77*/\ndydt[VKETM2]=+A85*T**n85*np.exp(-1*E85/R/T)*KETM2;\n/*ODEforVKETDM2=>VKETDM2=78*/\ndydt[VKETDM2]=+A86*T**n86*np.exp(-1*E86/R/T)*KETDM2;\n/*ODEforVSYNAPYL=>VSYNAPYL=79*/\ndydt[VSYNAPYL]=+A87*T**n87*np.exp(-1*E87/R/T)*SYNAPYL;\n/*ODEforMGUAI=>MGUAI=80*/\ndydt[MGUAI]=+A191*T**n191*np.exp(-1*E191/R/T)*RMGUAI*LIGM2+A170*T**n170*np.exp(-1*E170/R/T)*PLIGM2*RMGUAI+A296*T**n296*np.exp(-1*E296/R/T)*RMGUAI*C10H2+A401*T**n401*np.exp(-1*E401/R/T)*RMGUAI*KET+A317*T**n317*np.exp(-1*E317/R/T)*LIG*RMGUAI+\nA254*T**n254*np.exp(-1*E254/R/T)*RMGUAI*ADIOM2-A88*T**n88*np.exp(-1*E88/R/T)*MGUAI+A128*T**n128*np.exp(-1*E128/R/T)*RMGUAI*LIGH+A359*T**n359*np.exp(-1*E359/R/T)*PFET3*RMGUAI+A380*T**n380*np.exp(-1*E380/R/T)*RMGUAI*ADIO+A275*T**n275*np.exp(-1*E275/R/T)*KETM2*RMGUAI+\nA233*T**n233*np.exp(-1*E233/R/T)*RMGUAI*PFET3M2+A149*T**n149*np.exp(-1*E149/R/T)*RMGUAI*PLIGH+A212*T**n212*np.exp(-1*E212/R/T)*RMGUAI*LIGM2+A338*T**n338*np.exp(-1*E338/R/T)*LIG*RMGUAI;\n/*ODEforVMGUAI=>VMGUAI=81*/\ndydt[VMGUAI]=+A88*T**n88*np.exp(-1*E88/R/T)*MGUAI;\n/*ODEforVCOUMARYL=>VCOUMARYL=82*/\ndydt[VCOUMARYL]=+A89*T**n89*np.exp(-1*E89/R/T)*COUMARYL;\n/*ODEforVADIO=>VADIO=83*/\ndydt[VADIO]=+A90*T**n90*np.exp(-1*E90/R/T)*ADIO;\n/*ODEforVKET=>VKET=84*/\ndydt[VKET]=+A91*T**n91*np.exp(-1*E91/R/T)*KET;\n/*ODEforVKETD=>VKETD=85*/\ndydt[VKETD]=+A92*T**n92*np.exp(-1*E92/R/T)*KETD;\n/*ODEforPHENOL=>PHENOL=86*/\ndydt[PHENOL]=+A362*T**n362*np.exp(-1*E362/R/T)*PFET3*RPHENOL+A215*T**n215*np.exp(-1*E215/R/T)*LIGM2*RPHENOL+A278*T**n278*np.exp(-1*E278/R/T)*KETM2*RPHENOL+A173*T**n173*np.exp(-1*E173/R/T)*PLIGM2*RPHENOL+A341*T**n341*np.exp(-1*E341/R/T)*LIG*RPHENOL+\nA383*T**n383*np.exp(-1*E383/R/T)*ADIO*RPHENOL+A257*T**n257*np.exp(-1*E257/R/T)*ADIOM2*RPHENOL+A320*T**n320*np.exp(-1*E320/R/T)*LIG*RPHENOL-A93*T**n93*np.exp(-1*E93/R/T)*PHENOL+A131*T**n131*np.exp(-1*E131/R/T)*LIGH*RPHENOL+A194*T**n194*np.exp(-1*E194/R/T)*LIGM2*RPHENOL+\nA236*T**n236*np.exp(-1*E236/R/T)*PFET3M2*RPHENOL+A404*T**n404*np.exp(-1*E404/R/T)*KET*RPHENOL+A152*T**n152*np.exp(-1*E152/R/T)*PLIGH*RPHENOL+A299*T**n299*np.exp(-1*E299/R/T)*C10H2*RPHENOL;\n/*ODEforVPHENOL=>VPHENOL=87*/\ndydt[VPHENOL]=+A93*T**n93*np.exp(-1*E93/R/T)*PHENOL;\n/*ODEforLIGC=>LIGC=88*/\ndydt[LIGC]=-A108*T**n108*np.exp(-1*E108/R/T)*LIGC+A105*T**n105*np.exp(-1*E105/R/T)*PLIGC;\n/*ODEforLIGO=>LIGO=89*/\ndydt[LIGO]=-A107*T**n107*np.exp(-1*E107/R/T)*LIGO+A106*T**n106*np.exp(-1*E106/R/T)*PLIGO;\n/*ODEforC3H8O2=>C3H8O2=90*/\ndydt[C3H8O2]=+A114*T**n114*np.exp(-1*E114/R/T)*RC3H7O2*LIGH+A219*T**n219*np.exp(-1*E219/R/T)*RC3H7O2*PFET3M2+A156*T**n156*np.exp(-1*E156/R/T)*PLIGM2*RC3H7O2+A366*T**n366*np.exp(-1*E366/R/T)*RC3H7O2*ADIO+A345*T**n345*np.exp(-1*E345/R/T)*RC3H7O2*PFET3+\nA282*T**n282*np.exp(-1*E282/R/T)*RC3H7O2*C10H2+A135*T**n135*np.exp(-1*E135/R/T)*RC3H7O2*PLIGH+A303*T**n303*np.exp(-1*E303/R/T)*LIG*RC3H7O2+A198*T**n198*np.exp(-1*E198/R/T)*RC3H7O2*LIGM2+A177*T**n177*np.exp(-1*E177/R/T)*RC3H7O2*LIGM2+A261*T**n261*np.exp(-1*E261/R/T)*KETM2*RC3H7O2+\nA387*T**n387*np.exp(-1*E387/R/T)*RC3H7O2*KET+A240*T**n240*np.exp(-1*E240/R/T)*RC3H7O2*ADIOM2+A324*T**n324*np.exp(-1*E324/R/T)*LIG*RC3H7O2;\n/*ODEforCH4=>CH4=91*/\ndydt[CH4]=+A393*T**n393*np.exp(-1*E393/R/T)*RCH3*KET+A309*T**n309*np.exp(-1*E309/R/T)*LIG*RCH3+A183*T**n183*np.exp(-1*E183/R/T)*RCH3*LIGM2+A246*T**n246*np.exp(-1*E246/R/T)*RCH3*ADIOM2+A372*T**n372*np.exp(-1*E372/R/T)*RCH3*ADIO+\nA330*T**n330*np.exp(-1*E330/R/T)*LIG*RCH3+A351*T**n351*np.exp(-1*E351/R/T)*PFET3*RCH3+A162*T**n162*np.exp(-1*E162/R/T)*PLIGM2*RCH3+A267*T**n267*np.exp(-1*E267/R/T)*KETM2*RCH3+A225*T**n225*np.exp(-1*E225/R/T)*RCH3*PFET3M2+A288*T**n288*np.exp(-1*E288/R/T)*RCH3*C10H2+\nA120*T**n120*np.exp(-1*E120/R/T)*RCH3*LIGH+A204*T**n204*np.exp(-1*E204/R/T)*RCH3*LIGM2+A141*T**n141*np.exp(-1*E141/R/T)*PLIGH*RCH3;\n/*ODEforC3H4O=>C3H4O=92*/\ndydt[C3H4O]=+A229*T**n229*np.exp(-1*E229/R/T)*RC3H3O*PFET3M2+A208*T**n208*np.exp(-1*E208/R/T)*RC3H3O*LIGM2+A292*T**n292*np.exp(-1*E292/R/T)*RC3H3O*C10H2+A124*T**n124*np.exp(-1*E124/R/T)*RC3H3O*LIGH+A145*T**n145*np.exp(-1*E145/R/T)*RC3H3O*PLIGH+\nA397*T**n397*np.exp(-1*E397/R/T)*RC3H3O*KET+A250*T**n250*np.exp(-1*E250/R/T)*RC3H3O*ADIOM2+A334*T**n334*np.exp(-1*E334/R/T)*LIG*RC3H3O+A313*T**n313*np.exp(-1*E313/R/T)*LIG*RC3H3O+A187*T**n187*np.exp(-1*E187/R/T)*RC3H3O*LIGM2+A376*T**n376*np.exp(-1*E376/R/T)*RC3H3O*ADIO+\nA271*T**n271*np.exp(-1*E271/R/T)*KETM2*RC3H3O+A166*T**n166*np.exp(-1*E166/R/T)*PLIGM2*RC3H3O+A355*T**n355*np.exp(-1*E355/R/T)*RC3H3O*PFET3;\n'

In [14]:

    

# search the str that start with 'dydt[ADIO]=' and end with ';\n' 
# re.S means allow search for multi lines
re.findall(r'dydt[ADIO]=(.*);\n',reaction_file,re.S)


    

    
Out[14]:





[]

In [ ]:

    




    

In [ ]: