=== N(2,coll,const) variant A
  states = 12
  reactions = 24
  phi1 = 0.9545127980958669
  phi3 = 0.09262126532531366
  Rtot = 24046314446258627834704160606502759/50550731957549605604883273663125  (~475.686771)
  Ltot = 9090297098225022304132298635407/20020091864376081427676544025  (~454.058711)
  deg_wB = 3
  deg_wF = 2
  deg_p = 4
  squarefree = True
  roots_in_(0,Ltot) = 3
  roots_in_(0,inf) = 3
  p(L1) = 0  (~0.000000)
  p(L3) = 0  (~0.000000)
  sign p(LO) = True
  sign p(HI) = False
  roots_in_(LO,HI) = 1
  L2 = 0.0867243584148171
  p_prime(L2)<0 = True
  routh_degree_L1 = 11
  routh_all_positive_L1 = True
  lemma2.2_L1 = True
  rational_state_L1 = True
  max_re_nonzero_L1 = -0.0008470089887357591
  frac_C0_L1 = 0.9492448296955976
  routh_degree_L3 = 11
  routh_all_positive_L3 = True
  lemma2.2_L3 = True
  rational_state_L3 = True
  max_re_nonzero_L3 = -4.93921928147163e-05
  frac_T2_L3 = 0.9073740326322448
  L2_positive_eigs = [(0.004940170739752411+0j)]
  Rtot equals paper value: True
  Ltot equals paper value: True
=== N(2,kept,const) variant A
  states = 18
  reactions = 36
  phi1 = 0.9545717696898708
  phi3 = 0.09267013366532152
  Rtot = <long>  (~475.681195)
  Ltot = <long>  (~454.081440)
  deg_wB = 3
  deg_wF = 2
  deg_p = 4
  squarefree = True
  roots_in_(0,Ltot) = 3
  roots_in_(0,inf) = 3
  p(L1) = 0  (~0.000000)
  p(L3) = 0  (~0.000000)
  sign p(LO) = True
  sign p(HI) = False
  roots_in_(LO,HI) = 1
  L2 = 0.0868570169049325
  p_prime(L2)<0 = True
  routh_degree_L1 = 17
  routh_all_positive_L1 = True
  lemma2.2_L1 = True
  rational_state_L1 = True
  max_re_nonzero_L1 = -0.0008471903906227638
  frac_C0_L1 = 0.9492923920022183
  routh_degree_L3 = 17
  routh_all_positive_L3 = True
  lemma2.2_L3 = True
  rational_state_L3 = True
  max_re_nonzero_L3 = -4.9318157560542885e-05
  frac_T2_L3 = 0.9073251621378812
  L2_positive_eigs = [(0.004945556749642606+0j)]
=== N(3,coll,const) variant A
  states = 17
  reactions = 35
  phi1 = 0.9545127980090932
  phi3 = 0.09265444338471424
  Rtot = 308985806878679027119754571114189984611091/649532279582020475676233757039525805625  (~475.705083)
  Ltot = 11797485703149968343038054650385971589943/25981291183280819027049350281581032225  (~454.076190)
  deg_wB = 4
  deg_wF = 3
  deg_p = 5
  squarefree = True
  roots_in_(0,Ltot) = 3
  roots_in_(0,inf) = 3
  p(L1) = 0  (~0.000000)
  p(L3) = 0  (~0.000000)
  sign p(LO) = True
  sign p(HI) = False
  roots_in_(LO,HI) = 1
  L2 = 0.08672432534046563
  p_prime(L2)<0 = True
  routh_degree_L1 = 16
  routh_all_positive_L1 = True
  lemma2.2_L1 = True
  rational_state_L1 = True
  max_re_nonzero_L1 = -0.0008470089822414541
  frac_C0_L1 = 0.9492448296084345
  routh_degree_L3 = 16
  routh_all_positive_L3 = True
  lemma2.2_L3 = True
  rational_state_L3 = True
  max_re_nonzero_L3 = -4.939217941866793e-05
  frac_T2_L3 = 0.9073407655407558
  L2_positive_eigs = [(0.004940183492658633+0j)]
=== N(3,kept,const) variant A
  states = 26
  reactions = 53
  phi1 = 0.9545717696034085
  phi3 = 0.0927032713308495
  Rtot = <long>  (~475.699484)
  Ltot = <long>  (~454.098898)
  deg_wB = 4
  deg_wF = 3
  deg_p = 5
  squarefree = True
  roots_in_(0,Ltot) = 3
  roots_in_(0,inf) = 3
  p(L1) = 0  (~0.000000)
  p(L3) = 0  (~0.000000)
  sign p(LO) = True
  sign p(HI) = False
  roots_in_(LO,HI) = 1
  L2 = 0.0868569837584505
  p_prime(L2)<0 = True
  routh_degree_L1 = 25
  routh_all_positive_L1 = True
  lemma2.2_L1 = True
  rational_state_L1 = True
  max_re_nonzero_L1 = -0.0008471903835435446
  frac_C0_L1 = 0.9492923919152895
  routh_degree_L3 = 25
  routh_all_positive_L3 = True
  lemma2.2_L3 = True
  rational_state_L3 = True
  max_re_nonzero_L3 = -4.931814377219248e-05
  frac_T2_L3 = 0.9072919356219793
  L2_positive_eigs = [(0.004945569544443103+0j)]
=== N(3,kept,const) variant B
  states = 26
  reactions = 53
  phi1 = 0.9545717696138832
  phi3 = 0.0927070160470396
  Rtot = <long>  (~475.701551)
  Ltot = <long>  (~454.100871)
  deg_wB = 4
  deg_wF = 3
  deg_p = 5
  squarefree = True
  roots_in_(0,Ltot) = 3
  roots_in_(0,inf) = 3
  p(L1) = 0  (~0.000000)
  p(L3) = 0  (~0.000000)
  sign p(LO) = True
  sign p(HI) = False
  roots_in_(LO,HI) = 1
  L2 = 0.08685698021147209
  p_prime(L2)<0 = True
  routh_degree_L1 = 25
  routh_all_positive_L1 = True
  lemma2.2_L1 = True
  rational_state_L1 = True
  max_re_nonzero_L1 = -0.0008471903828329835
  frac_C0_L1 = 0.9492923918864612
  routh_degree_L3 = 25
  routh_all_positive_L3 = True
  lemma2.2_L3 = True
  rational_state_L3 = True
  max_re_nonzero_L3 = -4.931814120606042e-05
  frac_T2_L3 = 0.907288213066769
  L2_positive_eigs = [(0.004945570844401627+0j)]
