(* Generated by JWS Online *)

(* This is an experimental feature of JWS Online. Please report any mistakes.*)

(* Note that the following notable SBML entities or features are not supported in notebook outputyet: *)
    (* Events *)
    (* Constraints *)
    (* Units and UnitDefinitions *)
    (* AlgebraicRules *)
    (* conversionFactors *)

variables = {
AA\[LetterSpace]aw[t], AA\[LetterSpace]b[t], EO\[LetterSpace]a\[LetterSpace]aw[t], EO\[LetterSpace]a\[LetterSpace]b[t], EO\[LetterSpace]aw[t], EO\[LetterSpace]b[t], EO\[LetterSpace]bm[t], EO\[LetterSpace]i\[LetterSpace]aw[t], EO\[LetterSpace]i\[LetterSpace]b[t], HETE\[LetterSpace]aw[t], HETE\[LetterSpace]b[t], HPETE\[LetterSpace]aw[t], HPETE\[LetterSpace]b[t], Hn\[LetterSpace]aw[t], Hn\[LetterSpace]b[t], IL\[LetterSpace]aw[t], IL\[LetterSpace]b[t], IL\[LetterSpace]bm[t], LTA4\[LetterSpace]aw[t], LTA4\[LetterSpace]b[t], LTC4\[LetterSpace]aw[t], LTC4\[LetterSpace]aw\[LetterSpace]out[t], LTC4\[LetterSpace]b[t], LTC4\[LetterSpace]b\[LetterSpace]out[t], LTD4\[LetterSpace]aw[t], LTD4\[LetterSpace]b[t], LTE4\[LetterSpace]aw[t], LTE4\[LetterSpace]b[t], ML\[LetterSpace]blood[t], ML\[LetterSpace]intes[t], ZF\[LetterSpace]airways[t], ZF\[LetterSpace]blood[t], ZF\[LetterSpace]intes[t]
};

initialValues = {
AA\[LetterSpace]aw[0] == 0.1304887, AA\[LetterSpace]b[0] == 0.2890944, EO\[LetterSpace]a\[LetterSpace]aw[0] == 0.001358713, EO\[LetterSpace]a\[LetterSpace]b[0] == 0.0005383185, EO\[LetterSpace]aw[0] == 1.176558*^-06, EO\[LetterSpace]b[0] == 4.68114*^-07, EO\[LetterSpace]bm[0] == 1.637391*^-06, EO\[LetterSpace]i\[LetterSpace]aw[0] == 1.358698*^-07, EO\[LetterSpace]i\[LetterSpace]b[0] == 5.479176*^-08, HETE\[LetterSpace]aw[0] == 1.960926, HETE\[LetterSpace]b[0] == 5.037703, HPETE\[LetterSpace]aw[0] == 1.780604, HPETE\[LetterSpace]b[0] == 0.2965837, Hn\[LetterSpace]aw[0] == 14162.15, Hn\[LetterSpace]b[0] == 14995.88, IL\[LetterSpace]aw[0] == 0.6627439, IL\[LetterSpace]b[0] == 0.5994857, IL\[LetterSpace]bm[0] == 0.4023394, LTA4\[LetterSpace]aw[0] == 41.69257, LTA4\[LetterSpace]b[0] == 1.081167, LTC4\[LetterSpace]aw[0] == 6.806687, LTC4\[LetterSpace]aw\[LetterSpace]out[0] == 359.9744, LTC4\[LetterSpace]b[0] == 0.8869873, LTC4\[LetterSpace]b\[LetterSpace]out[0] == 2168.571, LTD4\[LetterSpace]aw[0] == 205.7602, LTD4\[LetterSpace]b[0] == 1308.488, LTE4\[LetterSpace]aw[0] == 322.6366, LTE4\[LetterSpace]b[0] == 2053.137, ML\[LetterSpace]blood[0] == 0.0, ML\[LetterSpace]intes[0] == 0.0, ZF\[LetterSpace]airways[0] == 0.0, ZF\[LetterSpace]blood[0] == 0.0, ZF\[LetterSpace]intes[0] == 0.0
};

rates = {
v1, v10in, v10out, v11, v12, v13, v14, v15, v16, v17, v18, v19, v2, v20, v21, v22, v23, v24, v25, v26, v27, v28, v29, v3, v30, v31, v32, v33, v34, v35, v36, v37, v38, v39, v4, v40in, v40out, v41, v42, v43, v44, v45, v46, v47, v48, v49, v5, v50, v51, v52, v53, v54, v55, v56, v57, v58, v59, v6, v60, v61, v62, v63, v64, v7, v8, v9
};

rateEquations = {
v1 -> Default*((PL*PLA2\[LetterSpace]Ca*Vmax\[LetterSpace]PLA2)/(Km\[LetterSpace]PLA2\[LetterSpace]APC + PL) - (V\[LetterSpace]CoA*AA\[LetterSpace]b[t])/(Km\[LetterSpace]CoA\[LetterSpace]AA + AA\[LetterSpace]b[t])), v10in -> Default*Kd50*LTC4\[LetterSpace]b[t], v10out -> 1000000*Default*Kd50*V\[LetterSpace]LTC\[LetterSpace]CB*LTC4\[LetterSpace]b[t], v11 -> fup\[LetterSpace]LT*k\[LetterSpace]ggt*Vd\[LetterSpace]LTC*LTC4\[LetterSpace]b\[LetterSpace]out[t], v12 -> fup\[LetterSpace]LT*k\[LetterSpace]dp*Vd\[LetterSpace]LTD*LTD4\[LetterSpace]b[t], v13 -> fup\[LetterSpace]LT*k\[LetterSpace]ltc\[LetterSpace]ltd\[LetterSpace]el*Vd\[LetterSpace]LTC*LTC4\[LetterSpace]b\[LetterSpace]out[t], v14 -> fup\[LetterSpace]LT*k\[LetterSpace]ltc\[LetterSpace]ltd\[LetterSpace]el*Vd\[LetterSpace]LTD*LTD4\[LetterSpace]b[t], v15 -> fup\[LetterSpace]LT*(k\[LetterSpace]acet + k\[LetterSpace]lte\[LetterSpace]el)*Vd\[LetterSpace]LTE*LTE4\[LetterSpace]b[t], v16 -> (ca*ka*OL\[LetterSpace]b^h\[LetterSpace]act*V\[LetterSpace]B*EO\[LetterSpace]b[t])/(EC50\[LetterSpace]act^h\[LetterSpace]act + OL\[LetterSpace]b^h\[LetterSpace]act), v17 -> ca*k\[LetterSpace]EO\[LetterSpace]m*V\[LetterSpace]B*EO\[LetterSpace]i\[LetterSpace]b[t], v18 -> ca*kia*V\[LetterSpace]B*EO\[LetterSpace]a\[LetterSpace]b[t], v19 -> k\[LetterSpace]EO\[LetterSpace]d*V\[LetterSpace]B*EO\[LetterSpace]b[t], v2 -> (Default*FLO3t\[LetterSpace]b*k\[LetterSpace]lo*AA\[LetterSpace]b[t])/K\[LetterSpace]AA, v20 -> k\[LetterSpace]EO\[LetterSpace]a\[LetterSpace]d*V\[LetterSpace]B*EO\[LetterSpace]a\[LetterSpace]b[t], v21 -> (k\[LetterSpace]EO\[LetterSpace]t\[LetterSpace]baw*Rec\[LetterSpace]occup\[LetterSpace]migr^h\[LetterSpace]migr*V\[LetterSpace]B*EO\[LetterSpace]i\[LetterSpace]b[t])/(EC50\[LetterSpace]migr^h\[LetterSpace]migr + Rec\[LetterSpace]occup\[LetterSpace]migr^h\[LetterSpace]migr), v22 -> (k\[LetterSpace]EO\[LetterSpace]t\[LetterSpace]baw*Rec\[LetterSpace]occup\[LetterSpace]migr^h\[LetterSpace]migr*V\[LetterSpace]B*EO\[LetterSpace]a\[LetterSpace]b[t])/(EC50\[LetterSpace]migr^h\[LetterSpace]migr + Rec\[LetterSpace]occup\[LetterSpace]migr^h\[LetterSpace]migr), v23 -> (k\[LetterSpace]EO\[LetterSpace]t\[LetterSpace]baw*Rec\[LetterSpace]occup\[LetterSpace]migr^h\[LetterSpace]migr*V\[LetterSpace]B*EO\[LetterSpace]b[t])/(EC50\[LetterSpace]migr^h\[LetterSpace]migr + Rec\[LetterSpace]occup\[LetterSpace]migr^h\[LetterSpace]migr), v24 -> k\[LetterSpace]Hn\[LetterSpace]p*V\[LetterSpace]B*(EO\[LetterSpace]a\[LetterSpace]b[t] + EO\[LetterSpace]b[t] + EO\[LetterSpace]i\[LetterSpace]b[t]), v25 -> fup\[LetterSpace]Hn*k\[LetterSpace]Hn\[LetterSpace]d*Vd\[LetterSpace]Hn*Hn\[LetterSpace]b[t], v26 -> k\[LetterSpace]IL\[LetterSpace]p*V\[LetterSpace]B*EO\[LetterSpace]a\[LetterSpace]b[t], v27 -> k\[LetterSpace]IL\[LetterSpace]d*Vd\[LetterSpace]IL5*IL\[LetterSpace]b[t], v28 -> k\[LetterSpace]IL\[LetterSpace]t\[LetterSpace]bbm*(IL\[LetterSpace]b[t] - IL\[LetterSpace]bm[t]) - J\[LetterSpace]BM\[LetterSpace]lymfl*IL\[LetterSpace]bm[t], v29 -> V\[LetterSpace]BM*(k1\[LetterSpace]min + (k1*IL\[LetterSpace]bm[t]^h\[LetterSpace]matur)/(Km\[LetterSpace]1^h\[LetterSpace]matur + IL\[LetterSpace]bm[t]^h\[LetterSpace]matur)), v3 -> Default*(1 + AA\[LetterSpace]b[t]/Ki\[LetterSpace]AA)*(FLO5HP\[LetterSpace]b*k\[LetterSpace]3 - FLO3t\[LetterSpace]b*k3*HPETE\[LetterSpace]b[t]), v30 -> ca*k\[LetterSpace]EO\[LetterSpace]t\[LetterSpace]bmb*V\[LetterSpace]BM*EO\[LetterSpace]bm[t], v31 -> Default*((PL*PLA2\[LetterSpace]Ca*Vmax\[LetterSpace]PLA2)/(Km\[LetterSpace]PLA2\[LetterSpace]APC + PL) - (V\[LetterSpace]CoA*AA\[LetterSpace]aw[t])/(Km\[LetterSpace]CoA\[LetterSpace]AA + AA\[LetterSpace]aw[t])), v32 -> (Default*FLO3t\[LetterSpace]aw*k\[LetterSpace]lo*AA\[LetterSpace]aw[t])/K\[LetterSpace]AA, v33 -> Default*(1 + AA\[LetterSpace]aw[t]/Ki\[LetterSpace]AA)*(FLO5HP\[LetterSpace]aw*k\[LetterSpace]3 - FLO3t\[LetterSpace]aw*k3*HPETE\[LetterSpace]aw[t]), v34 -> Default*r1*(-((FLO3\[LetterSpace]aw*HETE\[LetterSpace]aw[t])/Ke\[LetterSpace]ox) + FLO2\[LetterSpace]aw*HPETE\[LetterSpace]aw[t]), v35 -> (B\[LetterSpace]aw*Default*GPx)/A\[LetterSpace]aw, v36 -> Default*FLO5HP\[LetterSpace]aw*k\[LetterSpace]lta\[LetterSpace]syn, v37 -> (Default*nom\[LetterSpace]LTCs\[LetterSpace]aw)/den\[LetterSpace]LTCs\[LetterSpace]aw, v38 -> (B\[LetterSpace]hedh\[LetterSpace]aw*Default*HEDH5)/A\[LetterSpace]hedh\[LetterSpace]aw, v39 -> Default*Kd12*LTA4\[LetterSpace]aw[t], v4 -> Default*r1*(-((FLO3\[LetterSpace]b*HETE\[LetterSpace]b[t])/Ke\[LetterSpace]ox) + FLO2\[LetterSpace]b*HPETE\[LetterSpace]b[t]), v40in -> Default*Kd50*LTC4\[LetterSpace]aw[t], v40out -> 1000000*Default*Kd50*V\[LetterSpace]LTC\[LetterSpace]CAW*LTC4\[LetterSpace]aw[t], v41 -> k\[LetterSpace]ggt*Vd\[LetterSpace]AW\[LetterSpace]LTC*LTC4\[LetterSpace]aw\[LetterSpace]out[t], v42 -> k\[LetterSpace]dp*Vd\[LetterSpace]AW\[LetterSpace]LTD*LTD4\[LetterSpace]aw[t], v43 -> Q\[LetterSpace]AW\[LetterSpace]blf*R\[LetterSpace]LTE\[LetterSpace]B*((R\[LetterSpace]LTE\[LetterSpace]AW*LTE4\[LetterSpace]aw[t])/Kp\[LetterSpace]LTE\[LetterSpace]AW - LTE4\[LetterSpace]b[t]), v44 -> Q\[LetterSpace]AW\[LetterSpace]blf*R\[LetterSpace]LTD\[LetterSpace]B*((R\[LetterSpace]LTD\[LetterSpace]AW*LTD4\[LetterSpace]aw[t])/Kp\[LetterSpace]LTD\[LetterSpace]AW - LTD4\[LetterSpace]b[t]), v45 -> Q\[LetterSpace]AW\[LetterSpace]blf*R\[LetterSpace]LTC\[LetterSpace]B*((R\[LetterSpace]LTC\[LetterSpace]AW*LTC4\[LetterSpace]aw\[LetterSpace]out[t])/Kp\[LetterSpace]LTC\[LetterSpace]AW - LTC4\[LetterSpace]b\[LetterSpace]out[t]), v46 -> (ca*ka*OL\[LetterSpace]aw^h\[LetterSpace]act*V\[LetterSpace]AW*EO\[LetterSpace]aw[t])/(EC50\[LetterSpace]act^h\[LetterSpace]act + OL\[LetterSpace]aw^h\[LetterSpace]act), v47 -> ca*k\[LetterSpace]EO\[LetterSpace]m*V\[LetterSpace]AW*EO\[LetterSpace]i\[LetterSpace]aw[t], v48 -> ca*kia*V\[LetterSpace]AW*EO\[LetterSpace]a\[LetterSpace]aw[t], v49 -> k\[LetterSpace]EO\[LetterSpace]d*V\[LetterSpace]AW*EO\[LetterSpace]aw[t], v5 -> (B\[LetterSpace]b*Default*GPx)/A\[LetterSpace]b, v50 -> k\[LetterSpace]EO\[LetterSpace]a\[LetterSpace]d*V\[LetterSpace]AW*EO\[LetterSpace]a\[LetterSpace]aw[t], v51 -> k\[LetterSpace]Hn\[LetterSpace]p*V\[LetterSpace]AW*(EO\[LetterSpace]aw[t] + EO\[LetterSpace]a\[LetterSpace]aw[t] + EO\[LetterSpace]i\[LetterSpace]aw[t]), v52 -> Q\[LetterSpace]AW\[LetterSpace]blf*R\[LetterSpace]Hn\[LetterSpace]B*((R\[LetterSpace]Hn\[LetterSpace]AW*Hn\[LetterSpace]aw[t])/Kp\[LetterSpace]Hn\[LetterSpace]AW - Hn\[LetterSpace]b[t]), v53 -> k\[LetterSpace]IL\[LetterSpace]p*V\[LetterSpace]AW*EO\[LetterSpace]a\[LetterSpace]aw[t], v54 -> J\[LetterSpace]AW\[LetterSpace]lymfl*IL\[LetterSpace]aw[t] + k\[LetterSpace]IL\[LetterSpace]t\[LetterSpace]awb*(IL\[LetterSpace]aw[t] - IL\[LetterSpace]b[t]), v55 -> (k\[LetterSpace]acet + k\[LetterSpace]lte\[LetterSpace]el)*Vd\[LetterSpace]AW\[LetterSpace]LTE*LTE4\[LetterSpace]aw[t], v56 -> k\[LetterSpace]ltc\[LetterSpace]ltd\[LetterSpace]el*Vd\[LetterSpace]AW\[LetterSpace]LTD*LTD4\[LetterSpace]aw[t], v57 -> k\[LetterSpace]ltc\[LetterSpace]ltd\[LetterSpace]el*Vd\[LetterSpace]AW\[LetterSpace]LTC*LTC4\[LetterSpace]aw\[LetterSpace]out[t], v58 -> k\[LetterSpace]IL\[LetterSpace]d*V\[LetterSpace]AW*IL\[LetterSpace]aw[t], v59 -> k\[LetterSpace]Hn\[LetterSpace]d*Vd\[LetterSpace]AW\[LetterSpace]Hn*Hn\[LetterSpace]aw[t], v6 -> Default*FLO5HP\[LetterSpace]b*k\[LetterSpace]lta\[LetterSpace]syn, v60 -> Default*k\[LetterSpace]abs\[LetterSpace]zf*((1000*DOSE\[LetterSpace]zf*(1 - a + a*ft\[LetterSpace]zf)*F\[LetterSpace]zf*oral)/M\[LetterSpace]ZF + ZF\[LetterSpace]intes[t]), v61 -> Q\[LetterSpace]AW\[LetterSpace]blf*R\[LetterSpace]ZF\[LetterSpace]B*(-((R\[LetterSpace]ZF\[LetterSpace]AW*ZF\[LetterSpace]airways[t])/Kp\[LetterSpace]ZF\[LetterSpace]AW) + ZF\[LetterSpace]blood[t]), v62 -> k\[LetterSpace]elim\[LetterSpace]zf*Vd\[LetterSpace]ZF*ZF\[LetterSpace]blood[t], v63 -> Default*k\[LetterSpace]abs\[LetterSpace]ml*((1000000000*DOSE\[LetterSpace]ml*(1 - a + a*ft\[LetterSpace]ml)*F\[LetterSpace]ml*oral)/M\[LetterSpace]ML + ML\[LetterSpace]intes[t]), v64 -> k\[LetterSpace]elim\[LetterSpace]ml*Vd\[LetterSpace]ML*ML\[LetterSpace]blood[t], v7 -> (Default*nom\[LetterSpace]LTCs\[LetterSpace]b)/den\[LetterSpace]LTCs\[LetterSpace]b, v8 -> (B\[LetterSpace]hedh\[LetterSpace]b*Default*HEDH5)/A\[LetterSpace]hedh\[LetterSpace]b, v9 -> Default*Kd12*LTA4\[LetterSpace]b[t]
};

parameters = {
Ca -> 1.0, Ca\[LetterSpace]FEV\[LetterSpace]ex -> 10000.0, Cao\[LetterSpace]FEV -> 10000.0, DOSE\[LetterSpace]ml -> 0.0, DOSE\[LetterSpace]zf -> 0.0, EC50\[LetterSpace]ML\[LetterSpace]FEV -> 500.0, EC50\[LetterSpace]act -> 0.75, EC50\[LetterSpace]migr -> 0.115, Et\[LetterSpace]LTCs -> 1.0, FLOa -> 1.5, F\[LetterSpace]ml -> 0.660688, F\[LetterSpace]zf -> 0.082, GPx -> 1.6, GSH\[LetterSpace]aw -> 1000.0, GSH\[LetterSpace]b -> 5000.0, GS\[LetterSpace]pool\[LetterSpace]aw -> 10000.0, GS\[LetterSpace]pool\[LetterSpace]b -> 10000.0, HEDH5 -> 0.5, J\[LetterSpace]AW\[LetterSpace]lymfl -> 0.00115, J\[LetterSpace]BM\[LetterSpace]lymfl -> 0.00049, K\[LetterSpace]AA -> 10.74959, K\[LetterSpace]Ca2 -> 14.36738, K\[LetterSpace]Ca3 -> 7116.527, K\[LetterSpace]Ca\[LetterSpace]FEV -> 150000.0, K\[LetterSpace]GSH -> 744.9176, K\[LetterSpace]LTA -> 1.75951, K\[LetterSpace]LTA\[LetterSpace]GSH -> 1696.6, K\[LetterSpace]LTC -> 0.1951215, K\[LetterSpace]PLA2\[LetterSpace]Ca -> 0.1, Kd12 -> 0.007, Kd50 -> 0.43, KdZ -> 20.0, Kd\[LetterSpace]Hn\[LetterSpace]FEV -> 6300000.0, Kd\[LetterSpace]IL\[LetterSpace]migr -> 50.0, Kd\[LetterSpace]LT -> 1000.0, Kd\[LetterSpace]LTE\[LetterSpace]migr -> 50000.0, Kd\[LetterSpace]LTR1\[LetterSpace]FEV -> 1000.0, Kd\[LetterSpace]LTR2\[LetterSpace]FEV -> 10000.0, Kd\[LetterSpace]LT\[LetterSpace]2 -> 10000.0, Kd\[LetterSpace]gpx\[LetterSpace]GSSG -> 0.07154222, Kd\[LetterSpace]gpx\[LetterSpace]HETE5 -> 6.043446, Kd\[LetterSpace]hedh\[LetterSpace]HETE5 -> 0.331696, Kd\[LetterSpace]hedh\[LetterSpace]NADP -> 2.895899, Kd\[LetterSpace]hedh\[LetterSpace]NADPH -> 2.685581, Kd\[LetterSpace]hedh\[LetterSpace]oxoETE5 -> 1.667019, Ke\[LetterSpace]ox -> 99.99979, Ke\[LetterSpace]red -> 5.761955*^-07, Ki\[LetterSpace]AA -> 551.8748, Ki\[LetterSpace]HETE -> 0.5408177, Ki\[LetterSpace]ML\[LetterSpace]EOa -> 500.0, Km\[LetterSpace]1 -> 2.0, Km\[LetterSpace]CoA\[LetterSpace]AA -> 0.005, Km\[LetterSpace]PLA2\[LetterSpace]APC -> 20.0, Km\[LetterSpace]gpx\[LetterSpace]GSH -> 600.0, Km\[LetterSpace]gpx\[LetterSpace]HPETE5 -> 5.974381, Kp\[LetterSpace]Hn\[LetterSpace]AW -> 3950.0, Kp\[LetterSpace]LTC\[LetterSpace]AW -> 0.22, Kp\[LetterSpace]LTD\[LetterSpace]AW -> 0.22, Kp\[LetterSpace]LTE\[LetterSpace]AW -> 0.22, Kp\[LetterSpace]ZF\[LetterSpace]AW -> 0.204, Kp\[LetterSpace]ZF\[LetterSpace]IW -> 0.631, LOH\[LetterSpace]aw -> 0.0, LOH\[LetterSpace]b -> 0.0, LOOH\[LetterSpace]aw -> 20.0, LOOH\[LetterSpace]b -> 1.0, MAX\[LetterSpace]FEV -> 4.94, ML\[LetterSpace]airways -> 0.0, ML\[LetterSpace]ex -> 0.0, M\[LetterSpace]Hn -> 111.2, M\[LetterSpace]IL5 -> 45000.0, M\[LetterSpace]LTC -> 625.8, M\[LetterSpace]LTD -> 496.7, M\[LetterSpace]LTE -> 439.6, M\[LetterSpace]ML -> 586.18, M\[LetterSpace]ZF -> 236.0, NADP\[LetterSpace]aw -> 2900.0, NADP\[LetterSpace]b -> 2000.0, NP\[LetterSpace]pool\[LetterSpace]aw -> 3000.0, NP\[LetterSpace]pool\[LetterSpace]b -> 3000.0, N\[LetterSpace]A\[LetterSpace]pmole -> 602000000000.0, OL\[LetterSpace]b\[LetterSpace]ex -> 0.0, PL -> 110.0, Q\[LetterSpace]AW\[LetterSpace]blf -> 5.23, R1\[LetterSpace]portion\[LetterSpace]EOa -> 0.39, R1\[LetterSpace]portion\[LetterSpace]FEV -> 0.016, R2\[LetterSpace]portion\[LetterSpace]FEV -> 0.037, R\[LetterSpace]FEV -> 3.007, R\[LetterSpace]Hn\[LetterSpace]AW -> 5130.0, R\[LetterSpace]Hn\[LetterSpace]B -> 141.0, R\[LetterSpace]LTC\[LetterSpace]AW -> 1.4, R\[LetterSpace]LTC\[LetterSpace]B -> 0.538, R\[LetterSpace]LTD\[LetterSpace]AW -> 1.4, R\[LetterSpace]LTD\[LetterSpace]B -> 0.538, R\[LetterSpace]LTE\[LetterSpace]AW -> 1.4, R\[LetterSpace]LTE\[LetterSpace]B -> 0.538, R\[LetterSpace]ZF\[LetterSpace]AW -> 2.96, R\[LetterSpace]ZF\[LetterSpace]B -> 0.533, R\[LetterSpace]in\[LetterSpace]relax\[LetterSpace]FEV -> 2.073, T -> 1440.0, TSN\[LetterSpace]0 -> 15.5, V\[LetterSpace]CoA -> 350.0, Vmax\[LetterSpace]PLA2 -> 450.0, a -> 1.0, al -> 1.0, ca -> 10.0, diam\[LetterSpace]EO -> 0.00012, fup\[LetterSpace]Hn -> 0.77, fup\[LetterSpace]LT -> 0.16, fup\[LetterSpace]ML -> 0.004, fup\[LetterSpace]ZF -> 0.069, h\[LetterSpace]act -> 3.0, h\[LetterSpace]matur -> 1.0, h\[LetterSpace]migr -> 3.0, k1 -> 1*^-06, k1\[LetterSpace]min -> 1.6*^-07, k3 -> 34.0, k\[LetterSpace]3 -> 263640.0, k\[LetterSpace]EO\[LetterSpace]a\[LetterSpace]d -> 0.00015, k\[LetterSpace]EO\[LetterSpace]d -> 0.0003, k\[LetterSpace]EO\[LetterSpace]m -> 10.0, k\[LetterSpace]EO\[LetterSpace]t\[LetterSpace]baw -> 0.04, k\[LetterSpace]EO\[LetterSpace]t\[LetterSpace]bmb -> 0.02, k\[LetterSpace]Hn\[LetterSpace]d -> 0.033, k\[LetterSpace]Hn\[LetterSpace]p -> 18000000000.0, k\[LetterSpace]IL\[LetterSpace]d -> 0.0046, k\[LetterSpace]IL\[LetterSpace]p -> 16.0, k\[LetterSpace]IL\[LetterSpace]t\[LetterSpace]awb -> 0.05, k\[LetterSpace]IL\[LetterSpace]t\[LetterSpace]bbm -> 0.001, k\[LetterSpace]LTCs\[LetterSpace]back -> 1003.709, k\[LetterSpace]LTCs\[LetterSpace]fow -> 1068016.0, k\[LetterSpace]abs\[LetterSpace]ml -> 0.012, k\[LetterSpace]abs\[LetterSpace]zf -> 0.018, k\[LetterSpace]acet -> 0.002703885, k\[LetterSpace]dp -> 0.067, k\[LetterSpace]elim\[LetterSpace]ml -> 0.00225, k\[LetterSpace]elim\[LetterSpace]zf -> 0.004, k\[LetterSpace]fev\[LetterSpace]eff -> 3000000.0, k\[LetterSpace]ggt -> 0.1, k\[LetterSpace]gpx\[LetterSpace]cat -> 0.4884995, k\[LetterSpace]hedh\[LetterSpace]1 -> 88.33768, k\[LetterSpace]hedh\[LetterSpace]2 -> 1724.404, k\[LetterSpace]hedh\[LetterSpace]3 -> 31.49749, k\[LetterSpace]hedh\[LetterSpace]4 -> 8.077531, k\[LetterSpace]lo -> 4642.68, k\[LetterSpace]lta\[LetterSpace]syn -> 54420.0, k\[LetterSpace]ltc\[LetterSpace]ltd\[LetterSpace]el -> 0.1, k\[LetterSpace]lte\[LetterSpace]el -> 0.04, k\[LetterSpace]ox -> 0.0002658, k\[LetterSpace]ox2 -> 67.2, k\[LetterSpace]red -> 0.0002538, k\[LetterSpace]red2 -> 4.428*^-05, ka -> 500.0, ki -> 25000.0, kia -> 0.001, n\[LetterSpace]FEV -> 1.0, naEO\[LetterSpace]LTCsyn -> 0.0, npi -> 3.14, oral -> 1.0, oxoETE\[LetterSpace]aw -> 0.0, oxoETE\[LetterSpace]b -> 0.0, portion\[LetterSpace]migr -> 0.3, zf\[LetterSpace]inh -> 0.0, Default -> 1.0, V\[LetterSpace]AW -> 0.209999993443489, V\[LetterSpace]B -> 2.80999994277954, V\[LetterSpace]BM -> 0.824999988079071, Vd\[LetterSpace]AW\[LetterSpace]Hn -> 5640.0, Vd\[LetterSpace]AW\[LetterSpace]LTC -> 1.53999996185303, Vd\[LetterSpace]AW\[LetterSpace]LTD -> 1.53999996185303, Vd\[LetterSpace]AW\[LetterSpace]LTE -> 1.53999996185303, Vd\[LetterSpace]AW\[LetterSpace]ZF -> 3.25999999046326, Vd\[LetterSpace]Hn -> 78100.0, Vd\[LetterSpace]IL5 -> 10.1999998092651, Vd\[LetterSpace]LTC -> 9.47999954223633, Vd\[LetterSpace]LTD -> 9.47999954223633, Vd\[LetterSpace]LTE -> 9.47999954223633, Vd\[LetterSpace]ML -> 10.0, Vd\[LetterSpace]ZF -> 9.44999980926514
};

assignments = {
ft\[LetterSpace]ml -> Ceiling[t/T], ZF\[LetterSpace]blood\[LetterSpace]conc -> fup\[LetterSpace]ZF*ZF\[LetterSpace]blood[t], ZF\[LetterSpace]airways\[LetterSpace]conc -> ZF\[LetterSpace]airways[t], ML\[LetterSpace]blood\[LetterSpace]conc -> fup\[LetterSpace]ML*ML\[LetterSpace]blood[t], ML\[LetterSpace]airways\[LetterSpace]conc -> fup\[LetterSpace]ML*ML\[LetterSpace]blood[t], w\[LetterSpace]EO -> (diam\[LetterSpace]EO^3*npi)/6, V\[LetterSpace]aCB -> N\[LetterSpace]A\[LetterSpace]pmole*V\[LetterSpace]B*w\[LetterSpace]EO*EO\[LetterSpace]a\[LetterSpace]b[t], V\[LetterSpace]aCAW -> N\[LetterSpace]A\[LetterSpace]pmole*V\[LetterSpace]AW*w\[LetterSpace]EO*EO\[LetterSpace]a\[LetterSpace]aw[t], GSSG\[LetterSpace]b -> (-GSH\[LetterSpace]b + GS\[LetterSpace]pool\[LetterSpace]b)/2, PLA2\[LetterSpace]D -> 1 + Ca/K\[LetterSpace]PLA2\[LetterSpace]Ca, Ki\[LetterSpace]AA\[LetterSpace]AA -> Ki\[LetterSpace]AA, Ki\[LetterSpace]HPETE\[LetterSpace]AA -> Ki\[LetterSpace]AA, OOH\[LetterSpace]b -> LOOH\[LetterSpace]b + HPETE\[LetterSpace]b[t], OH\[LetterSpace]b -> LOH\[LetterSpace]b + HETE\[LetterSpace]b[t], r1 -> k\[LetterSpace]ox + (Ca*k\[LetterSpace]ox2)/K\[LetterSpace]Ca2, r2 -> k\[LetterSpace]red + (Ca*k\[LetterSpace]red2)/K\[LetterSpace]Ca2, RELFLO5\[LetterSpace]b -> ((k\[LetterSpace]lo*AA\[LetterSpace]b[t])/K\[LetterSpace]AA + k3*(1 + AA\[LetterSpace]b[t]/Ki\[LetterSpace]AA)*HPETE\[LetterSpace]b[t])/(k\[LetterSpace]lta\[LetterSpace]syn + k\[LetterSpace]3*(1 + AA\[LetterSpace]b[t]/Ki\[LetterSpace]HPETE\[LetterSpace]AA)), FLO3\[LetterSpace]b -> FLOa/DFLOa\[LetterSpace]b, FLO2\[LetterSpace]b -> FLO3\[LetterSpace]b*REDOX\[LetterSpace]b, FLO3t\[LetterSpace]b -> FLO3\[LetterSpace]b*(1 + Ca/K\[LetterSpace]Ca3), FLO2t\[LetterSpace]b -> FLO2\[LetterSpace]b*(1 + Ca/K\[LetterSpace]Ca2), FLO5HP\[LetterSpace]b -> FLO3t\[LetterSpace]b*RELFLO5\[LetterSpace]b, C\[LetterSpace]b -> (GSH\[LetterSpace]b^2*HPETE\[LetterSpace]b[t])/(Km\[LetterSpace]gpx\[LetterSpace]GSH^2*Km\[LetterSpace]gpx\[LetterSpace]HPETE5), delta\[LetterSpace]LTCs\[LetterSpace]b -> 1 + GSH\[LetterSpace]b/K\[LetterSpace]GSH + LTA4\[LetterSpace]b[t]/K\[LetterSpace]LTA + (GSH\[LetterSpace]b*LTA4\[LetterSpace]b[t])/(K\[LetterSpace]LTA*K\[LetterSpace]LTA\[LetterSpace]GSH) + LTC4\[LetterSpace]b[t]/K\[LetterSpace]LTC, C\[LetterSpace]hedh\[LetterSpace]b -> ((k\[LetterSpace]hedh\[LetterSpace]2*NADP\[LetterSpace]b)/Kd\[LetterSpace]hedh\[LetterSpace]NADP + (k\[LetterSpace]hedh\[LetterSpace]3*oxoETE\[LetterSpace]b)/Kd\[LetterSpace]hedh\[LetterSpace]oxoETE5)*(1 + NADPH\[LetterSpace]b/Kd\[LetterSpace]hedh\[LetterSpace]NADPH + HETE\[LetterSpace]b[t]/Kd\[LetterSpace]hedh\[LetterSpace]HETE5), GSSG\[LetterSpace]aw -> (-GSH\[LetterSpace]aw + GS\[LetterSpace]pool\[LetterSpace]aw)/2, NADPH\[LetterSpace]aw -> -NADP\[LetterSpace]aw + NP\[LetterSpace]pool\[LetterSpace]aw, OOH\[LetterSpace]aw -> LOOH\[LetterSpace]aw + HPETE\[LetterSpace]aw[t], OH\[LetterSpace]aw -> LOH\[LetterSpace]aw + HETE\[LetterSpace]aw[t], RELFLO5\[LetterSpace]aw -> ((k\[LetterSpace]lo*AA\[LetterSpace]aw[t])/K\[LetterSpace]AA + k3*(1 + AA\[LetterSpace]aw[t]/Ki\[LetterSpace]AA)*HPETE\[LetterSpace]aw[t])/(k\[LetterSpace]lta\[LetterSpace]syn + k\[LetterSpace]3*(1 + AA\[LetterSpace]aw[t]/Ki\[LetterSpace]HPETE\[LetterSpace]AA)), FLO3\[LetterSpace]aw -> FLOa/DFLOa\[LetterSpace]aw, FLO2\[LetterSpace]aw -> FLO3\[LetterSpace]aw*REDOX\[LetterSpace]aw, FLO3t\[LetterSpace]aw -> FLO3\[LetterSpace]aw*(1 + Ca/K\[LetterSpace]Ca3), FLO2t\[LetterSpace]aw -> FLO2\[LetterSpace]aw*(1 + Ca/K\[LetterSpace]Ca2), FLO5HP\[LetterSpace]aw -> FLO3t\[LetterSpace]aw*RELFLO5\[LetterSpace]aw, C\[LetterSpace]aw -> (GSH\[LetterSpace]aw^2*HPETE\[LetterSpace]aw[t])/(Km\[LetterSpace]gpx\[LetterSpace]GSH^2*Km\[LetterSpace]gpx\[LetterSpace]HPETE5), B\[LetterSpace]aw -> (GSH\[LetterSpace]aw^2*k\[LetterSpace]gpx\[LetterSpace]cat*HPETE\[LetterSpace]aw[t])/(Km\[LetterSpace]gpx\[LetterSpace]GSH^2*Km\[LetterSpace]gpx\[LetterSpace]HPETE5), den\[LetterSpace]LTCs\[LetterSpace]aw -> delta\[LetterSpace]LTCs\[LetterSpace]aw*K\[LetterSpace]GSH*K\[LetterSpace]LTA\[LetterSpace]GSH*K\[LetterSpace]LTC, C\[LetterSpace]hedh\[LetterSpace]aw -> ((k\[LetterSpace]hedh\[LetterSpace]2*NADP\[LetterSpace]aw)/Kd\[LetterSpace]hedh\[LetterSpace]NADP + (k\[LetterSpace]hedh\[LetterSpace]3*oxoETE\[LetterSpace]aw)/Kd\[LetterSpace]hedh\[LetterSpace]oxoETE5)*(1 + NADPH\[LetterSpace]aw/Kd\[LetterSpace]hedh\[LetterSpace]NADPH + HETE\[LetterSpace]aw[t]/Kd\[LetterSpace]hedh\[LetterSpace]HETE5), B\[LetterSpace]hedh\[LetterSpace]aw -> -((k\[LetterSpace]hedh\[LetterSpace]3*k\[LetterSpace]hedh\[LetterSpace]4*NADPH\[LetterSpace]aw*oxoETE\[LetterSpace]aw)/(Kd\[LetterSpace]hedh\[LetterSpace]NADPH*Kd\[LetterSpace]hedh\[LetterSpace]oxoETE5)) + (k\[LetterSpace]hedh\[LetterSpace]1*k\[LetterSpace]hedh\[LetterSpace]2*NADP\[LetterSpace]aw*HETE\[LetterSpace]aw[t])/(Kd\[LetterSpace]hedh\[LetterSpace]HETE5*Kd\[LetterSpace]hedh\[LetterSpace]NADP), Rec\[LetterSpace]occup\[LetterSpace]migr -> ((1 - portion\[LetterSpace]migr)*IL\[LetterSpace]b[t])/(Kd\[LetterSpace]IL\[LetterSpace]migr + IL\[LetterSpace]b[t]) + (fup\[LetterSpace]LT*portion\[LetterSpace]migr*LTE4\[LetterSpace]b[t])/(Kd\[LetterSpace]LTE\[LetterSpace]migr + fup\[LetterSpace]LT*LTE4\[LetterSpace]b[t]), r\[LetterSpace]in\[LetterSpace]FEV -> -((R\[LetterSpace]FEV - R\[LetterSpace]in\[LetterSpace]relax\[LetterSpace]FEV)*Sqrt[R\[LetterSpace]FEV/r\[LetterSpace]out\[LetterSpace]FEV]) + r\[LetterSpace]out\[LetterSpace]FEV, FEV1 -> (MAX\[LetterSpace]FEV*r\[LetterSpace]in\[LetterSpace]FEV^2)/R\[LetterSpace]in\[LetterSpace]relax\[LetterSpace]FEV^2, time\[LetterSpace]hour -> t/60, time\[LetterSpace]day -> 0.0006944444444444445*t, N\[LetterSpace]EO\[LetterSpace]bm -> 6.02*^8*EO\[LetterSpace]bm[t], N\[LetterSpace]EO\[LetterSpace]b -> 6.02*^8*EO\[LetterSpace]b[t], N\[LetterSpace]EO\[LetterSpace]i\[LetterSpace]b -> 6.02*^8*EO\[LetterSpace]i\[LetterSpace]b[t], EO\[LetterSpace]b\[LetterSpace]tot -> EO\[LetterSpace]a\[LetterSpace]b[t] + EO\[LetterSpace]b[t] + EO\[LetterSpace]i\[LetterSpace]b[t], N\[LetterSpace]EO\[LetterSpace]b\[LetterSpace]tot -> 6.02*^8*(EO\[LetterSpace]a\[LetterSpace]b[t] + EO\[LetterSpace]b[t] + EO\[LetterSpace]i\[LetterSpace]b[t]), N\[LetterSpace]EO\[LetterSpace]aw -> 6.02*^8*EO\[LetterSpace]aw[t], N\[LetterSpace]EO\[LetterSpace]i\[LetterSpace]aw -> 6.02*^8*EO\[LetterSpace]i\[LetterSpace]aw[t], N\[LetterSpace]EO\[LetterSpace]a\[LetterSpace]aw -> 6.02*^8*EO\[LetterSpace]a\[LetterSpace]aw[t], EO\[LetterSpace]aw\[LetterSpace]tot -> EO\[LetterSpace]aw[t] + EO\[LetterSpace]a\[LetterSpace]aw[t] + EO\[LetterSpace]i\[LetterSpace]aw[t], N\[LetterSpace]EO\[LetterSpace]aw\[LetterSpace]tot -> 6.02*^8*(EO\[LetterSpace]aw[t] + EO\[LetterSpace]a\[LetterSpace]aw[t] + EO\[LetterSpace]i\[LetterSpace]aw[t]), EO\[LetterSpace]b\[LetterSpace]tot\[LetterSpace]per\[LetterSpace]ss -> 187617.26078799248*(-0.000533 + EO\[LetterSpace]b\[LetterSpace]tot), N\[LetterSpace]EO\[LetterSpace]aw\[LetterSpace]tot\[LetterSpace]perc -> 0.00012213964225298783*N\[LetterSpace]EO\[LetterSpace]aw\[LetterSpace]tot, N\[LetterSpace]EO\[LetterSpace]b\[LetterSpace]tot\[LetterSpace]perc -> 0.000023441047289734392*N\[LetterSpace]EO\[LetterSpace]b\[LetterSpace]tot, LTC4\[LetterSpace]b\[LetterSpace]pM -> 1000000*LTC4\[LetterSpace]b[t], N\[LetterSpace]EO\[LetterSpace]a -> (N\[LetterSpace]EO\[LetterSpace]a\[LetterSpace]aw*V\[LetterSpace]AW + N\[LetterSpace]EO\[LetterSpace]a\[LetterSpace]b*V\[LetterSpace]B)/(V\[LetterSpace]AW + V\[LetterSpace]B), N\[LetterSpace]EO -> (N\[LetterSpace]EO\[LetterSpace]aw\[LetterSpace]tot*V\[LetterSpace]AW + N\[LetterSpace]EO\[LetterSpace]b\[LetterSpace]tot*V\[LetterSpace]B)/(V\[LetterSpace]AW + V\[LetterSpace]B), N\[LetterSpace]EO\[LetterSpace]aw\[LetterSpace]perc -> 0.00012213964225298783*N\[LetterSpace]EO\[LetterSpace]aw\[LetterSpace]tot, Hn\[LetterSpace]aw\[LetterSpace]perc -> 0.00706114955514758*Hn\[LetterSpace]aw[t], N\[LetterSpace]EO\[LetterSpace]perc -> 0.0002787394288071625*N\[LetterSpace]EO, OL\[LetterSpace]ASM -> ((1 - R1\[LetterSpace]portion\[LetterSpace]FEV - R2\[LetterSpace]portion\[LetterSpace]FEV)*Hn\[LetterSpace]aw[t])/(Kd\[LetterSpace]Hn\[LetterSpace]FEV*(1 + Hn\[LetterSpace]aw[t]/Kd\[LetterSpace]Hn\[LetterSpace]FEV)) + (R1\[LetterSpace]portion\[LetterSpace]FEV*LTD4\[LetterSpace]aw[t])/(Kd\[LetterSpace]LTR1\[LetterSpace]FEV*(1 + ML\[LetterSpace]airways\[LetterSpace]conc/EC50\[LetterSpace]ML\[LetterSpace]FEV + LTD4\[LetterSpace]aw[t]/Kd\[LetterSpace]LTR1\[LetterSpace]FEV)) + (R2\[LetterSpace]portion\[LetterSpace]FEV*(LTC4\[LetterSpace]aw\[LetterSpace]out[t] + LTD4\[LetterSpace]aw[t]))/(Kd\[LetterSpace]LTR2\[LetterSpace]FEV*(1 + (LTC4\[LetterSpace]aw\[LetterSpace]out[t] + LTD4\[LetterSpace]aw[t])/Kd\[LetterSpace]LTR2\[LetterSpace]FEV)), Ca\[LetterSpace]FEV\[LetterSpace]LTR2 -> (k\[LetterSpace]fev\[LetterSpace]eff*(((1 - R1\[LetterSpace]portion\[LetterSpace]FEV - R2\[LetterSpace]portion\[LetterSpace]FEV)*Hn\[LetterSpace]aw[t])/(Kd\[LetterSpace]Hn\[LetterSpace]FEV*(1 + Hn\[LetterSpace]aw[t]/Kd\[LetterSpace]Hn\[LetterSpace]FEV)) + (R2\[LetterSpace]portion\[LetterSpace]FEV*(LTC4\[LetterSpace]aw\[LetterSpace]out[t] + LTD4\[LetterSpace]aw[t]))/(Kd\[LetterSpace]LTR2\[LetterSpace]FEV*(1 + (LTC4\[LetterSpace]aw\[LetterSpace]out[t] + LTD4\[LetterSpace]aw[t])/Kd\[LetterSpace]LTR2\[LetterSpace]FEV))))/(-Cao\[LetterSpace]FEV + Ca\[LetterSpace]FEV), Ca\[LetterSpace]FEV\[LetterSpace]LTR1 -> (k\[LetterSpace]fev\[LetterSpace]eff*R1\[LetterSpace]portion\[LetterSpace]FEV*LTD4\[LetterSpace]aw[t])/((-Cao\[LetterSpace]FEV + Ca\[LetterSpace]FEV)*Kd\[LetterSpace]LTR1\[LetterSpace]FEV*(1 + ML\[LetterSpace]airways\[LetterSpace]conc/EC50\[LetterSpace]ML\[LetterSpace]FEV + LTD4\[LetterSpace]aw[t]/Kd\[LetterSpace]LTR1\[LetterSpace]FEV)), LTD4\[LetterSpace]aw\[LetterSpace]pers -> 0.4860031104199067*LTD4\[LetterSpace]aw[t], ML\[LetterSpace]uM -> ML\[LetterSpace]blood[t]/1000000, LTD4\[LetterSpace]b\[LetterSpace]free -> fup\[LetterSpace]LT*LTD4\[LetterSpace]b[t], LTs\[LetterSpace]aw\[LetterSpace]pg -> ((M\[LetterSpace]LTC*LTC4\[LetterSpace]aw\[LetterSpace]out[t])/fup\[LetterSpace]LT + (M\[LetterSpace]LTD*LTD4\[LetterSpace]aw[t])/fup\[LetterSpace]LT + (M\[LetterSpace]LTE*LTE4\[LetterSpace]aw[t])/fup\[LetterSpace]LT)/1000, N\[LetterSpace]EO\[LetterSpace]a\[LetterSpace]b -> 6.02*^8*EO\[LetterSpace]a\[LetterSpace]b[t], FEV1\[LetterSpace]percent -> 28.34467120181406*(-3.528 + FEV1), TSN\[LetterSpace]Hn -> (Ca\[LetterSpace]FEV^n\[LetterSpace]FEV*(R\[LetterSpace]FEV - (R\[LetterSpace]FEV*(R\[LetterSpace]FEV - R\[LetterSpace]in\[LetterSpace]relax\[LetterSpace]FEV)^2)^(1/3))*TSN\[LetterSpace]0)/(Ca\[LetterSpace]FEV^n\[LetterSpace]FEV + K\[LetterSpace]Ca\[LetterSpace]FEV^n\[LetterSpace]FEV), TSN -> (Ca\[LetterSpace]FEV\[LetterSpace]ex^n\[LetterSpace]FEV*(R\[LetterSpace]FEV - (R\[LetterSpace]FEV*(R\[LetterSpace]FEV - R\[LetterSpace]in\[LetterSpace]relax\[LetterSpace]FEV)^2)^(1/3))*TSN\[LetterSpace]0)/(Ca\[LetterSpace]FEV\[LetterSpace]ex^n\[LetterSpace]FEV + K\[LetterSpace]Ca\[LetterSpace]FEV^n\[LetterSpace]FEV), r\[LetterSpace]out\[LetterSpace]FEV -> R\[LetterSpace]FEV - (Ca\[LetterSpace]FEV^n\[LetterSpace]FEV*(R\[LetterSpace]FEV - (R\[LetterSpace]FEV*(R\[LetterSpace]FEV - R\[LetterSpace]in\[LetterSpace]relax\[LetterSpace]FEV)^2)^(1/3)))/(Ca\[LetterSpace]FEV^n\[LetterSpace]FEV + K\[LetterSpace]Ca\[LetterSpace]FEV^n\[LetterSpace]FEV), Ca\[LetterSpace]FEV -> Cao\[LetterSpace]FEV + k\[LetterSpace]fev\[LetterSpace]eff*(((1 - R1\[LetterSpace]portion\[LetterSpace]FEV - R2\[LetterSpace]portion\[LetterSpace]FEV)*Hn\[LetterSpace]aw[t])/(Kd\[LetterSpace]Hn\[LetterSpace]FEV*(1 + Hn\[LetterSpace]aw[t]/Kd\[LetterSpace]Hn\[LetterSpace]FEV)) + (R1\[LetterSpace]portion\[LetterSpace]FEV*LTD4\[LetterSpace]aw[t])/(Kd\[LetterSpace]LTR1\[LetterSpace]FEV*(1 + ML\[LetterSpace]airways\[LetterSpace]conc/EC50\[LetterSpace]ML\[LetterSpace]FEV + LTD4\[LetterSpace]aw[t]/Kd\[LetterSpace]LTR1\[LetterSpace]FEV)) + (R2\[LetterSpace]portion\[LetterSpace]FEV*(LTC4\[LetterSpace]aw\[LetterSpace]out[t] + LTD4\[LetterSpace]aw[t]))/(Kd\[LetterSpace]LTR2\[LetterSpace]FEV*(1 + (LTC4\[LetterSpace]aw\[LetterSpace]out[t] + LTD4\[LetterSpace]aw[t])/Kd\[LetterSpace]LTR2\[LetterSpace]FEV))), OL\[LetterSpace]aw -> (R1\[LetterSpace]portion\[LetterSpace]EOa*LTD4\[LetterSpace]aw[t])/(Kd\[LetterSpace]LT*(1 + ML\[LetterSpace]airways\[LetterSpace]conc/Ki\[LetterSpace]ML\[LetterSpace]EOa + LTD4\[LetterSpace]aw[t]/Kd\[LetterSpace]LT)) + ((1 - R1\[LetterSpace]portion\[LetterSpace]EOa)*(LTC4\[LetterSpace]aw\[LetterSpace]out[t] + LTD4\[LetterSpace]aw[t]))/(Kd\[LetterSpace]LT\[LetterSpace]2*(1 + (LTC4\[LetterSpace]aw\[LetterSpace]out[t] + LTD4\[LetterSpace]aw[t])/Kd\[LetterSpace]LT\[LetterSpace]2)), OL\[LetterSpace]b -> (fup\[LetterSpace]LT*R1\[LetterSpace]portion\[LetterSpace]EOa*LTD4\[LetterSpace]b[t])/(Kd\[LetterSpace]LT*(1 + ML\[LetterSpace]blood\[LetterSpace]conc/Ki\[LetterSpace]ML\[LetterSpace]EOa + (fup\[LetterSpace]LT*LTD4\[LetterSpace]b[t])/Kd\[LetterSpace]LT)) + (fup\[LetterSpace]LT*(1 - R1\[LetterSpace]portion\[LetterSpace]EOa)*(LTC4\[LetterSpace]b\[LetterSpace]out[t] + LTD4\[LetterSpace]b[t]))/(Kd\[LetterSpace]LT\[LetterSpace]2*(1 + (fup\[LetterSpace]LT*(LTC4\[LetterSpace]b\[LetterSpace]out[t] + LTD4\[LetterSpace]b[t]))/Kd\[LetterSpace]LT\[LetterSpace]2)), A\[LetterSpace]hedh\[LetterSpace]aw -> C\[LetterSpace]hedh\[LetterSpace]aw + (1 + NADP\[LetterSpace]aw/Kd\[LetterSpace]hedh\[LetterSpace]NADP + oxoETE\[LetterSpace]aw/Kd\[LetterSpace]hedh\[LetterSpace]oxoETE5)*((k\[LetterSpace]hedh\[LetterSpace]4*NADPH\[LetterSpace]aw)/Kd\[LetterSpace]hedh\[LetterSpace]NADPH + (k\[LetterSpace]hedh\[LetterSpace]1*HETE\[LetterSpace]aw[t])/Kd\[LetterSpace]hedh\[LetterSpace]HETE5), nom\[LetterSpace]LTCs\[LetterSpace]aw -> Et\[LetterSpace]LTCs*(GSH\[LetterSpace]aw*K\[LetterSpace]LTC*k\[LetterSpace]LTCs\[LetterSpace]fow*LTA4\[LetterSpace]aw[t] - K\[LetterSpace]GSH*K\[LetterSpace]LTA\[LetterSpace]GSH*k\[LetterSpace]LTCs\[LetterSpace]back*LTC4\[LetterSpace]aw[t]), delta\[LetterSpace]LTCs\[LetterSpace]aw -> 1 + GSH\[LetterSpace]aw/K\[LetterSpace]GSH + LTA4\[LetterSpace]aw[t]/K\[LetterSpace]LTA + (GSH\[LetterSpace]aw*LTA4\[LetterSpace]aw[t])/(K\[LetterSpace]LTA*K\[LetterSpace]LTA\[LetterSpace]GSH) + LTC4\[LetterSpace]aw[t]/K\[LetterSpace]LTC, A\[LetterSpace]aw -> C\[LetterSpace]aw + (GSH\[LetterSpace]aw^2*(1 + GSSG\[LetterSpace]aw/Kd\[LetterSpace]gpx\[LetterSpace]GSSG))/Km\[LetterSpace]gpx\[LetterSpace]GSH^2 + ((1 + HETE\[LetterSpace]aw[t]/Kd\[LetterSpace]gpx\[LetterSpace]HETE5)*HPETE\[LetterSpace]aw[t])/Km\[LetterSpace]gpx\[LetterSpace]HPETE5, DFLOa\[LetterSpace]aw -> (1 + Ca/K\[LetterSpace]Ca3)*(1 + ((1 + Ca/K\[LetterSpace]Ca2)*REDOX\[LetterSpace]aw)/(1 + Ca/K\[LetterSpace]Ca3) + (al*ZF\[LetterSpace]airways\[LetterSpace]conc)/KdZ + (AA\[LetterSpace]aw[t]*(1 + AA\[LetterSpace]aw[t]/Ki\[LetterSpace]AA\[LetterSpace]AA))/K\[LetterSpace]AA + RELFLO5\[LetterSpace]aw*(1 + AA\[LetterSpace]aw[t]/Ki\[LetterSpace]HPETE\[LetterSpace]AA) + HETE\[LetterSpace]aw[t]/Ki\[LetterSpace]HETE + (AA\[LetterSpace]aw[t]*(1 + HETE\[LetterSpace]aw[t]/K\[LetterSpace]AA\[LetterSpace]HETE))/Ki\[LetterSpace]AA), REDOX\[LetterSpace]aw -> ((OH\[LetterSpace]aw*r1)/Ke\[LetterSpace]ox + r2/Ke\[LetterSpace]red + (al*ki*(1 + Ca/K\[LetterSpace]Ca3)*ZF\[LetterSpace]airways\[LetterSpace]conc)/KdZ)/(OOH\[LetterSpace]aw*r1 + r2), B\[LetterSpace]hedh\[LetterSpace]b -> -((k\[LetterSpace]hedh\[LetterSpace]3*k\[LetterSpace]hedh\[LetterSpace]4*NADPH\[LetterSpace]b*oxoETE\[LetterSpace]b)/(Kd\[LetterSpace]hedh\[LetterSpace]NADPH*Kd\[LetterSpace]hedh\[LetterSpace]oxoETE5)) + (k\[LetterSpace]hedh\[LetterSpace]1*k\[LetterSpace]hedh\[LetterSpace]2*NADP\[LetterSpace]b*HETE\[LetterSpace]b[t])/(Kd\[LetterSpace]hedh\[LetterSpace]HETE5*Kd\[LetterSpace]hedh\[LetterSpace]NADP), A\[LetterSpace]hedh\[LetterSpace]b -> C\[LetterSpace]hedh\[LetterSpace]b + (1 + NADP\[LetterSpace]b/Kd\[LetterSpace]hedh\[LetterSpace]NADP + oxoETE\[LetterSpace]b/Kd\[LetterSpace]hedh\[LetterSpace]oxoETE5)*((k\[LetterSpace]hedh\[LetterSpace]4*NADPH\[LetterSpace]b)/Kd\[LetterSpace]hedh\[LetterSpace]NADPH + (k\[LetterSpace]hedh\[LetterSpace]1*HETE\[LetterSpace]b[t])/Kd\[LetterSpace]hedh\[LetterSpace]HETE5), den\[LetterSpace]LTCs\[LetterSpace]b -> delta\[LetterSpace]LTCs\[LetterSpace]b*K\[LetterSpace]GSH*K\[LetterSpace]LTA\[LetterSpace]GSH*K\[LetterSpace]LTC, nom\[LetterSpace]LTCs\[LetterSpace]b -> Et\[LetterSpace]LTCs*(GSH\[LetterSpace]b*K\[LetterSpace]LTC*k\[LetterSpace]LTCs\[LetterSpace]fow*LTA4\[LetterSpace]b[t] - K\[LetterSpace]GSH*K\[LetterSpace]LTA\[LetterSpace]GSH*k\[LetterSpace]LTCs\[LetterSpace]back*LTC4\[LetterSpace]b[t]), B\[LetterSpace]b -> (GSH\[LetterSpace]b^2*k\[LetterSpace]gpx\[LetterSpace]cat*HPETE\[LetterSpace]b[t])/(Km\[LetterSpace]gpx\[LetterSpace]GSH^2*Km\[LetterSpace]gpx\[LetterSpace]HPETE5), A\[LetterSpace]b -> C\[LetterSpace]b + (GSH\[LetterSpace]b^2*(1 + GSSG\[LetterSpace]b/Kd\[LetterSpace]gpx\[LetterSpace]GSSG))/Km\[LetterSpace]gpx\[LetterSpace]GSH^2 + ((1 + HETE\[LetterSpace]b[t]/Kd\[LetterSpace]gpx\[LetterSpace]HETE5)*HPETE\[LetterSpace]b[t])/Km\[LetterSpace]gpx\[LetterSpace]HPETE5, DFLOa\[LetterSpace]b -> (1 + Ca/K\[LetterSpace]Ca3)*(1 + ((1 + Ca/K\[LetterSpace]Ca2)*REDOX\[LetterSpace]b)/(1 + Ca/K\[LetterSpace]Ca3) + (al*ZF\[LetterSpace]blood\[LetterSpace]conc)/KdZ + (AA\[LetterSpace]b[t]*(1 + AA\[LetterSpace]b[t]/Ki\[LetterSpace]AA\[LetterSpace]AA))/K\[LetterSpace]AA + RELFLO5\[LetterSpace]b*(1 + AA\[LetterSpace]b[t]/Ki\[LetterSpace]HPETE\[LetterSpace]AA) + HETE\[LetterSpace]b[t]/Ki\[LetterSpace]HETE + (AA\[LetterSpace]b[t]*(1 + HETE\[LetterSpace]b[t]/K\[LetterSpace]AA\[LetterSpace]HETE))/Ki\[LetterSpace]AA), REDOX\[LetterSpace]b -> ((OH\[LetterSpace]b*r1)/Ke\[LetterSpace]ox + r2/Ke\[LetterSpace]red + (al*ki*(1 + Ca/K\[LetterSpace]Ca3)*ZF\[LetterSpace]blood\[LetterSpace]conc)/KdZ)/(OOH\[LetterSpace]b*r1 + r2), K\[LetterSpace]AA\[LetterSpace]HETE -> Ki\[LetterSpace]HETE, PLA2\[LetterSpace]Ca -> Ca/(K\[LetterSpace]PLA2\[LetterSpace]Ca*PLA2\[LetterSpace]D), NADPH\[LetterSpace]b -> -NADP\[LetterSpace]b + NP\[LetterSpace]pool\[LetterSpace]b, V\[LetterSpace]LTC\[LetterSpace]CAW -> N\[LetterSpace]A\[LetterSpace]pmole*V\[LetterSpace]AW*w\[LetterSpace]EO*(EO\[LetterSpace]a\[LetterSpace]aw[t] + naEO\[LetterSpace]LTCsyn*(EO\[LetterSpace]aw[t] + EO\[LetterSpace]i\[LetterSpace]aw[t])), V\[LetterSpace]CAW -> N\[LetterSpace]A\[LetterSpace]pmole*V\[LetterSpace]AW*w\[LetterSpace]EO*(EO\[LetterSpace]aw[t] + EO\[LetterSpace]a\[LetterSpace]aw[t] + EO\[LetterSpace]i\[LetterSpace]aw[t]), V\[LetterSpace]LTC\[LetterSpace]CB -> N\[LetterSpace]A\[LetterSpace]pmole*V\[LetterSpace]B*w\[LetterSpace]EO*(EO\[LetterSpace]a\[LetterSpace]b[t] + naEO\[LetterSpace]LTCsyn*(EO\[LetterSpace]b[t] + EO\[LetterSpace]i\[LetterSpace]b[t])), V\[LetterSpace]CB -> N\[LetterSpace]A\[LetterSpace]pmole*V\[LetterSpace]B*w\[LetterSpace]EO*(EO\[LetterSpace]a\[LetterSpace]b[t] + EO\[LetterSpace]b[t] + EO\[LetterSpace]i\[LetterSpace]b[t]), ft\[LetterSpace]zf -> Ceiling[(-960. + t)/T] + Ceiling[(-660. + t)/T] + Ceiling[(-360. + t)/T] + Ceiling[t/T]
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

units = {
{"time" -> "", "metabolite" -> "", "extent" -> ""}
};

(* Time evolution *)
odes = {
AA\[LetterSpace]aw'[t] == 1.0*v31 -1.0*v32, AA\[LetterSpace]b'[t] == 1.0*v1 -1.0*v2, EO\[LetterSpace]a\[LetterSpace]aw'[t] == 1.0*v22 +1.0*v47 -1.0*v48 -1.0*v50, EO\[LetterSpace]a\[LetterSpace]b'[t] == 1.0*v17 -1.0*v18 -1.0*v20 -1.0*v22, EO\[LetterSpace]aw'[t] == 1.0*v23 +1.0*v48 -1.0*v46 -1.0*v49, EO\[LetterSpace]b'[t] == 1.0*v18 +1.0*v30 -1.0*v16 -1.0*v19 -1.0*v23, EO\[LetterSpace]bm'[t] == 1.0*v29 -1.0*v30, EO\[LetterSpace]i\[LetterSpace]aw'[t] == 1.0*v21 +1.0*v46 -1.0*v47, EO\[LetterSpace]i\[LetterSpace]b'[t] == 1.0*v16 -1.0*v17 -1.0*v21, HETE\[LetterSpace]aw'[t] == 1.0*v34 +1.0*v35 -1.0*v38, HETE\[LetterSpace]b'[t] == 1.0*v4 +1.0*v5 -1.0*v8, HPETE\[LetterSpace]aw'[t] == 1.0*v33 -1.0*v34 -1.0*v35, HPETE\[LetterSpace]b'[t] == 1.0*v3 -1.0*v4 -1.0*v5, Hn\[LetterSpace]aw'[t] == 1.0*v51 -1.0*v52 -1.0*v59, Hn\[LetterSpace]b'[t] == 1.0*v24 +1.0*v52 -1.0*v25, IL\[LetterSpace]aw'[t] == 1.0*v53 -1.0*v54 -1.0*v58, IL\[LetterSpace]b'[t] == 1.0*v26 +1.0*v54 -1.0*v27 -1.0*v28, IL\[LetterSpace]bm'[t] == 1.0*v28 , LTA4\[LetterSpace]aw'[t] == 1.0*v36 -1.0*v37 -1.0*v39, LTA4\[LetterSpace]b'[t] == 1.0*v6 -1.0*v7 -1.0*v9, LTC4\[LetterSpace]aw'[t] == 1.0*v37 -1.0*v40in, LTC4\[LetterSpace]aw\[LetterSpace]out'[t] == 1.0*v40out -1.0*v41 -1.0*v45 -1.0*v57, LTC4\[LetterSpace]b'[t] == 1.0*v7 -1.0*v10in, LTC4\[LetterSpace]b\[LetterSpace]out'[t] == 1.0*v10out +1.0*v45 -1.0*v11 -1.0*v13, LTD4\[LetterSpace]aw'[t] == 1.0*v41 -1.0*v42 -1.0*v44 -1.0*v56, LTD4\[LetterSpace]b'[t] == 1.0*v11 +1.0*v44 -1.0*v12 -1.0*v14, LTE4\[LetterSpace]aw'[t] == 1.0*v42 -1.0*v43 -1.0*v55, LTE4\[LetterSpace]b'[t] == 1.0*v12 +1.0*v43 -1.0*v15, ML\[LetterSpace]blood'[t] == 1.0*v63 -1.0*v64, ML\[LetterSpace]intes'[t] ==  -1.0*v63, ZF\[LetterSpace]airways'[t] == 1.0*v61 , ZF\[LetterSpace]blood'[t] == 1.0*v60 -1.0*v61 -1.0*v62, ZF\[LetterSpace]intes'[t] ==  -1.0*v60
};
timeCourse = NDSolve[Join[odes, initialValues]//.rateEquations//.assignments//.parameters, variables, {t, 0, 100}];

(* Steady-state solution initialized with result of time evolution *)
findRootEquations = odes /.D[_[t],t]->0;
findRootVariables = Partition[Flatten[{#, #/.timeCourse/.t->100} &/@variables],2];
steadyStateVariables = FindRoot[findRootEquations//.rateEquations//.assignments//.parameters, findRootVariables, MaxIterations->100]
fluxes = #//.assignments//.parameters/.steadyStateVariables&/@rateEquations

(* Plot the time evolution of the variables *)
plotTable=Table[Plot[variables[[i]]/.parameters/.timeCourse,{t,0,100},PlotLegends->variables[[i]],PlotRange->Full],{i,Length[variables]}]