network (Lck, ZAP-70 constant, no CD8), ZAP-docking levels = 2 , Lck-MM intermediates: False
== reduced network: species=13 complexes=15 linkage=2 rank=11 deficiency=2 reactions=24 weakly_rev=False
   linkage-class deficiencies [0, 0]  sum 0  terminal SLCs per LC [1, 1]
   Deficiency One Theorem hypotheses satisfied: False
rate constants (irreversible reactions):
   J0+     TCR + pMHC1 -> Complex                   k = 100
   J0-     Complex -> TCR + pMHC1                   k = 1/100
   J185+   TCR_ZAP -> TCR                           k = 1/10
   J10+    Complex_ZAP -> TCR_ZAP + pMHC1           k = 1
   J10-    TCR_ZAP + pMHC1 -> Complex_ZAP           k = 1
   J5+     Complex_ITAM1_PP -> Complex_ZAP          k = 1/10
   J5-     Complex_ZAP -> Complex_ITAM1_PP          k = 1/100
   J4+     Complex_ITAM1_P -> Complex_ITAM1_PP      k = 1
   J2+     Complex -> Complex_ITAM1_P               k = 1/100
   J7+     Complex_ZAP -> Complex_ZAP_ITAM2_P       k = 100
   J9+     Complex_ZAP_ITAM2_P -> Complex_ZAP_ITAM2_PP k = 1
   J15+    Complex_ZAP_ITAM2_PP -> TCR_ZAP + pMHC1  k = 1/100
   J14+    Complex_ZAP_ITAM2_P -> TCR_ZAP + pMHC1   k = 1
   J248+   TCR_2ZAP -> TCR_ZAP                      k = 1/1000
   J12+    Complex_ITAM1_P -> TCR + pMHC1           k = 1
   J13+    Complex_ITAM1_PP -> TCR + pMHC1          k = 10
   J240+   Complex_ZAP_ITAM2_PP -> Complex_2ZAP     k = 1
   J240-   Complex_2ZAP -> Complex_ZAP_ITAM2_PP     k = 1/10
   J245+   Complex_2ZAP -> TCR_2ZAP + pMHC1         k = 100
   J245-   TCR_2ZAP + pMHC1 -> Complex_2ZAP         k = 1/100
   J243+   Complex_2ZAP -> Complex_2ZAP_ITAM3_P     k = 1
   J244+   Complex_2ZAP_ITAM3_P -> Complex_2ZAP_ITAM3_PP k = 1
   J247+   Complex_2ZAP_ITAM3_PP -> TCR_2ZAP + pMHC1 k = 1
   J246+   Complex_2ZAP_ITAM3_P -> TCR_2ZAP + pMHC1 k = 1
numerical: local max phi=0.962951 at L=2.884e-02 ; subsequent min phi=0.061285 at L=1.380e+02
L1 = 6/625  L3 = 410   phi(L1) = 0.954512798096  phi(L3) = 0.0926212653253
R_tot = 24046314446258627834704160606502759/50550731957549605604883273663125  (~475.687)
L_tot = 9090297098225022304132298635407/20020091864376081427676544025  (~454.059)
exact steady state at L = 6/625 : f(x)=0 verified exactly; totals R_tot, L_tot match; all 13 coordinates positive
exact steady state at L = 410 : f(x)=0 verified exactly; totals R_tot, L_tot match; all 13 coordinates positive
g(L)=L+R_tot*phi(L)-L_tot : g(99/10000)>0 and g(11/125)<0 exactly -> a third steady state with L in between
third root L2 in [0.0867243584148, 0.0867243584149]
steady state at L1: reduced char. poly degree 11; Routh-Hurwitz first column all positive: True
steady state at L3: reduced char. poly degree 11; Routh-Hurwitz first column all positive: True
middle state (approximate): eigenvalues with largest real parts: 0.00494+0i, -1.312e-17+6.285e-16i, -1.312e-17-6.285e-16i, -0.1439+0i
