network (Lck, ZAP-70 constant, no CD8), ZAP-docking levels = 2 , Lck-MM intermediates: True
== reduced network: species=19 complexes=21 linkage=2 rank=17 deficiency=2 reactions=36 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):
   J1+     Complex -> Node3                         k = 1/100
   J1-     Node3 -> Complex                         k = 1
   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+     Node6 -> Complex_ITAM1_PP                k = 1000
   J2+     Node3 -> Complex_ITAM1_P                 k = 1000
   J3+     Complex_ITAM1_P -> Node6                 k = 1
   J3-     Node6 -> Complex_ITAM1_P                 k = 1
   J6+     Complex_ZAP -> Node10                    k = 100
   J6-     Node10 -> Complex_ZAP                    k = 1
   J7+     Node10 -> Complex_ZAP_ITAM2_P            k = 1000
   J8+     Complex_ZAP_ITAM2_P -> Node12            k = 1
   J8-     Node12 -> Complex_ZAP_ITAM2_P            k = 1
   J9+     Node12 -> Complex_ZAP_ITAM2_PP           k = 1000
   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
   J241+   Complex_2ZAP -> Node136                  k = 1
   J241-   Node136 -> Complex_2ZAP                  k = 1
   J243+   Node136 -> Complex_2ZAP_ITAM3_P          k = 1000
   J242+   Complex_2ZAP_ITAM3_P -> Node138          k = 1
   J242-   Node138 -> Complex_2ZAP_ITAM3_P          k = 1
   J244+   Node138 -> Complex_2ZAP_ITAM3_PP         k = 1000
   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.963013 at L=2.884e-02 ; subsequent min phi=0.061360 at L=1.380e+02
L1 = 6/625  L3 = 410   phi(L1) = 0.954571769690  phi(L3) = 0.0926701336653
R_tot = 11456783810427636498997799321512053421471070184609042665842968867631039/24085004697151811179984203864147357338161622920897738030139022472500  (~475.681)
L_tot = 21873107225838526440543244021435624414115180845753136869150835074448/48170009394303622359968407728294714676323245841795476060278044945  (~454.081)
exact steady state at L = 6/625 : f(x)=0 verified exactly; totals R_tot, L_tot match; all 19 coordinates positive
exact steady state at L = 410 : f(x)=0 verified exactly; totals R_tot, L_tot match; all 19 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.0868570169049, 0.0868570169049]
steady state at L1: reduced char. poly degree 17; Routh-Hurwitz first column all positive: True
steady state at L3: reduced char. poly degree 17; Routh-Hurwitz first column all positive: True
middle state (approximate): eigenvalues with largest real parts: 0.004946+0i, 3.555e-15+0i, -2.29e-14+0i, -0.144+0i
