(* 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 = {
A\[LetterSpace]b[t], A\[LetterSpace]m[t], A\[LetterSpace]p[t], A\[LetterSpace]r[t], FAAH\[LetterSpace]b[t], FAAH\[LetterSpace]m[t], FAAH\[LetterSpace]r[t], FAAHinh\[LetterSpace]b[t], FAAHinh\[LetterSpace]m[t], FAAHinh\[LetterSpace]r[t], L\[LetterSpace]b[t], L\[LetterSpace]m[t], L\[LetterSpace]p[t], L\[LetterSpace]r[t], NAPE\[LetterSpace]b[t], NAPE\[LetterSpace]r[t], NLPE\[LetterSpace]b[t], NLPE\[LetterSpace]r[t], NOPE\[LetterSpace]b[t], NOPE\[LetterSpace]r[t], NPPE\[LetterSpace]b[t], NPPE\[LetterSpace]r[t], NSPE\[LetterSpace]b[t], NSPE\[LetterSpace]r[t], O\[LetterSpace]b[t], O\[LetterSpace]m[t], O\[LetterSpace]p[t], O\[LetterSpace]r[t], PFM\[LetterSpace]gut[t], PFM\[LetterSpace]p[t], PFM\[LetterSpace]r[t], P\[LetterSpace]b[t], P\[LetterSpace]m[t], P\[LetterSpace]p[t], P\[LetterSpace]r[t], S\[LetterSpace]b[t], S\[LetterSpace]m[t], S\[LetterSpace]p[t], S\[LetterSpace]r[t]
};

initialValues = {
A\[LetterSpace]b[0] == 0.7493309, A\[LetterSpace]m[0] == 0.97761, A\[LetterSpace]p[0] == 0.8740574, A\[LetterSpace]r[0] == 0.5419204, FAAH\[LetterSpace]b[0] == 15.366, FAAH\[LetterSpace]m[0] == 10.686, FAAH\[LetterSpace]r[0] == 2.165868, FAAHinh\[LetterSpace]b[0] == 0.0, FAAHinh\[LetterSpace]m[0] == 0.0, FAAHinh\[LetterSpace]r[0] == 0.0, L\[LetterSpace]b[0] == 2.319571, L\[LetterSpace]m[0] == 0.0, L\[LetterSpace]p[0] == 1.916254, L\[LetterSpace]r[0] == 1.705466, NAPE\[LetterSpace]b[0] == 3.879041*^-05, NAPE\[LetterSpace]r[0] == 4.241633*^-06, NLPE\[LetterSpace]b[0] == 7.550331*^-05, NLPE\[LetterSpace]r[0] == 8.256095*^-06, NOPE\[LetterSpace]b[0] == 0.0008814287, NOPE\[LetterSpace]r[0] == 9.638198*^-05, NPPE\[LetterSpace]b[0] == 0.0001732296, NPPE\[LetterSpace]r[0] == 1.894222*^-05, NSPE\[LetterSpace]b[0] == 0.0001272629, NSPE\[LetterSpace]r[0] == 1.391587*^-05, O\[LetterSpace]b[0] == 20.77858, O\[LetterSpace]m[0] == 16.3219, O\[LetterSpace]p[0] == 5.085073, O\[LetterSpace]r[0] == 14.23822, PFM\[LetterSpace]gut[0] == 0.0, PFM\[LetterSpace]p[0] == 0.0, PFM\[LetterSpace]r[0] == 0.0, P\[LetterSpace]b[0] == 6.541209, P\[LetterSpace]m[0] == 5.809415, P\[LetterSpace]p[0] == 4.849307, P\[LetterSpace]r[0] == 4.121915, S\[LetterSpace]b[0] == 3.427807, S\[LetterSpace]m[0] == 2.968774, S\[LetterSpace]p[0] == 0.273772, S\[LetterSpace]r[0] == 2.515968
};

rates = {
absorp, dist, elim, vA\[LetterSpace]UE\[LetterSpace]b, vA\[LetterSpace]UE\[LetterSpace]r, vA\[LetterSpace]b\[LetterSpace]m, vA\[LetterSpace]degr\[LetterSpace]b, vA\[LetterSpace]degr\[LetterSpace]m, vA\[LetterSpace]degr\[LetterSpace]r, vA\[LetterSpace]m\[LetterSpace]p, vA\[LetterSpace]r\[LetterSpace]p, vA\[LetterSpace]syn\[LetterSpace]b, vA\[LetterSpace]syn\[LetterSpace]r, vFAAH\[LetterSpace]degr\[LetterSpace]b, vFAAH\[LetterSpace]degr\[LetterSpace]m, vFAAH\[LetterSpace]degr\[LetterSpace]r, vFAAH\[LetterSpace]inh\[LetterSpace]b, vFAAH\[LetterSpace]inh\[LetterSpace]degr\[LetterSpace]b, vFAAH\[LetterSpace]inh\[LetterSpace]degr\[LetterSpace]m, vFAAH\[LetterSpace]inh\[LetterSpace]degr\[LetterSpace]r, vFAAH\[LetterSpace]inh\[LetterSpace]m, vFAAH\[LetterSpace]inh\[LetterSpace]r, vFAAH\[LetterSpace]syn\[LetterSpace]b, vFAAH\[LetterSpace]syn\[LetterSpace]m, vFAAH\[LetterSpace]syn\[LetterSpace]r, vL\[LetterSpace]UE\[LetterSpace]b, vL\[LetterSpace]UE\[LetterSpace]r, vL\[LetterSpace]b\[LetterSpace]m, vL\[LetterSpace]degr\[LetterSpace]b, vL\[LetterSpace]degr\[LetterSpace]m, vL\[LetterSpace]degr\[LetterSpace]r, vL\[LetterSpace]m\[LetterSpace]p, vL\[LetterSpace]r\[LetterSpace]p, vL\[LetterSpace]syn\[LetterSpace]b, vL\[LetterSpace]syn\[LetterSpace]r, vNAPE\[LetterSpace]syn\[LetterSpace]b, vNAPE\[LetterSpace]syn\[LetterSpace]r, vNLPE\[LetterSpace]syn\[LetterSpace]b, vNLPE\[LetterSpace]syn\[LetterSpace]r, vNOPE\[LetterSpace]syn\[LetterSpace]b, vNOPE\[LetterSpace]syn\[LetterSpace]r, vNPPE\[LetterSpace]syn\[LetterSpace]b, vNPPE\[LetterSpace]syn\[LetterSpace]r, vNSPE\[LetterSpace]syn\[LetterSpace]b, vNSPE\[LetterSpace]syn\[LetterSpace]r, vO\[LetterSpace]UE\[LetterSpace]b, vO\[LetterSpace]UE\[LetterSpace]r, vO\[LetterSpace]b\[LetterSpace]m, vO\[LetterSpace]degr\[LetterSpace]b, vO\[LetterSpace]degr\[LetterSpace]m, vO\[LetterSpace]degr\[LetterSpace]r, vO\[LetterSpace]r\[LetterSpace]p, vO\[LetterSpace]syn\[LetterSpace]b, vO\[LetterSpace]syn\[LetterSpace]r, vP\[LetterSpace]UE\[LetterSpace]b, vP\[LetterSpace]UE\[LetterSpace]r, vP\[LetterSpace]b\[LetterSpace]m, vP\[LetterSpace]degr\[LetterSpace]b, vP\[LetterSpace]degr\[LetterSpace]m, vP\[LetterSpace]degr\[LetterSpace]r, vP\[LetterSpace]m\[LetterSpace]p, vP\[LetterSpace]r\[LetterSpace]p, vP\[LetterSpace]syn\[LetterSpace]b, vP\[LetterSpace]syn\[LetterSpace]r, vS\[LetterSpace]UE\[LetterSpace]b, vS\[LetterSpace]UE\[LetterSpace]r, vS\[LetterSpace]b\[LetterSpace]m, vS\[LetterSpace]degr\[LetterSpace]b, vS\[LetterSpace]degr\[LetterSpace]m, vS\[LetterSpace]degr\[LetterSpace]r, vS\[LetterSpace]m\[LetterSpace]p, vS\[LetterSpace]r\[LetterSpace]p, vS\[LetterSpace]syn\[LetterSpace]b, vS\[LetterSpace]syn\[LetterSpace]r, vo\[LetterSpace]m\[LetterSpace]p
};

rateEquations = {
absorp -> kabs\[LetterSpace]PFM*MD, dist -> kout\[LetterSpace]PFM*PFM\[LetterSpace]p[t] - kin\[LetterSpace]PFM*PFM\[LetterSpace]r[t], elim -> klinear\[LetterSpace]PFM*PFM\[LetterSpace]p[t] + (Vm\[LetterSpace]PFM*PFM\[LetterSpace]p[t])/(Vss\[LetterSpace]PFM*(Km\[LetterSpace]PFM + PFM\[LetterSpace]p[t]/Vss\[LetterSpace]PFM)), vA\[LetterSpace]UE\[LetterSpace]b -> BRAIN*b\[LetterSpace]FAAH\[LetterSpace]Brain*kcl\[LetterSpace]A*A\[LetterSpace]b[t], vA\[LetterSpace]UE\[LetterSpace]r -> c\[LetterSpace]NAAA\[LetterSpace]ROB*kcl\[LetterSpace]A*A\[LetterSpace]r[t], vA\[LetterSpace]b\[LetterSpace]m -> (ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]A*MEC*(A\[LetterSpace]b[t] - A\[LetterSpace]m[t]))/(Km\[LetterSpace]p\[LetterSpace]m\[LetterSpace]A + A\[LetterSpace]b[t] + A\[LetterSpace]m[t]), vA\[LetterSpace]degr\[LetterSpace]b -> (a\[LetterSpace]FAAH\[LetterSpace]A*BRAIN*kcat\[LetterSpace]FAAH*A\[LetterSpace]b[t]*FAAH\[LetterSpace]b[t])/(FAAH\[LetterSpace]D\[LetterSpace]b*Km\[LetterSpace]FAAH\[LetterSpace]A), vA\[LetterSpace]degr\[LetterSpace]m -> (a\[LetterSpace]FAAH\[LetterSpace]A*kcat\[LetterSpace]FAAH*MEC*A\[LetterSpace]m[t]*FAAH\[LetterSpace]m[t])/(FAAH\[LetterSpace]D\[LetterSpace]m*Km\[LetterSpace]FAAH\[LetterSpace]A), vA\[LetterSpace]degr\[LetterSpace]r -> (a\[LetterSpace]FAAH\[LetterSpace]A*kcat\[LetterSpace]FAAH*ROB*A\[LetterSpace]r[t]*FAAH\[LetterSpace]r[t])/(FAAH\[LetterSpace]D\[LetterSpace]r*Km\[LetterSpace]FAAH\[LetterSpace]A), vA\[LetterSpace]m\[LetterSpace]p -> (ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]A*MEC*(A\[LetterSpace]m[t] - Ktr\[LetterSpace]p\[LetterSpace]m\[LetterSpace]A*A\[LetterSpace]p[t]))/(Km\[LetterSpace]p\[LetterSpace]m\[LetterSpace]A + A\[LetterSpace]m[t] + A\[LetterSpace]p[t]), vA\[LetterSpace]r\[LetterSpace]p -> (ktr\[LetterSpace]r\[LetterSpace]p*PLASMA*(-(Ktr\[LetterSpace]p\[LetterSpace]r\[LetterSpace]A*A\[LetterSpace]p[t]) + A\[LetterSpace]r[t]))/(Km\[LetterSpace]p\[LetterSpace]m\[LetterSpace]A + A\[LetterSpace]p[t] + A\[LetterSpace]r[t]), vA\[LetterSpace]syn\[LetterSpace]b -> (BRAIN*k\[LetterSpace]NA\[LetterSpace]PE*PLD\[LetterSpace]b*NAPE\[LetterSpace]b[t])/(den\[LetterSpace]b*Km\[LetterSpace]NA\[LetterSpace]PE), vA\[LetterSpace]syn\[LetterSpace]r -> (k\[LetterSpace]NA\[LetterSpace]PE*PLD\[LetterSpace]r*ROB*NAPE\[LetterSpace]r[t])/(den\[LetterSpace]r*Km\[LetterSpace]NA\[LetterSpace]PE), vFAAH\[LetterSpace]degr\[LetterSpace]b -> BRAIN*k\[LetterSpace]deg\[LetterSpace]FAAH*FAAH\[LetterSpace]b[t], vFAAH\[LetterSpace]degr\[LetterSpace]m -> k\[LetterSpace]deg\[LetterSpace]FAAH*MEC*FAAH\[LetterSpace]m[t], vFAAH\[LetterSpace]degr\[LetterSpace]r -> k\[LetterSpace]deg\[LetterSpace]FAAH*ROB*FAAH\[LetterSpace]r[t], vFAAH\[LetterSpace]inh\[LetterSpace]b -> BRAIN*k\[LetterSpace]inh*PF\[LetterSpace]b*FAAH\[LetterSpace]b[t], vFAAH\[LetterSpace]inh\[LetterSpace]degr\[LetterSpace]b -> BRAIN*k\[LetterSpace]deg\[LetterSpace]FAAH*FAAHinh\[LetterSpace]b[t], vFAAH\[LetterSpace]inh\[LetterSpace]degr\[LetterSpace]m -> k\[LetterSpace]deg\[LetterSpace]FAAH*MEC*FAAHinh\[LetterSpace]m[t], vFAAH\[LetterSpace]inh\[LetterSpace]degr\[LetterSpace]r -> k\[LetterSpace]deg\[LetterSpace]FAAH*ROB*FAAHinh\[LetterSpace]r[t], vFAAH\[LetterSpace]inh\[LetterSpace]m -> k\[LetterSpace]inh*MEC*PF\[LetterSpace]m*FAAH\[LetterSpace]m[t], vFAAH\[LetterSpace]inh\[LetterSpace]r -> k\[LetterSpace]inh*PF\[LetterSpace]r*ROB*FAAH\[LetterSpace]r[t], vFAAH\[LetterSpace]syn\[LetterSpace]b -> BRAIN*b\[LetterSpace]FAAH\[LetterSpace]Brain*FAAH\[LetterSpace]t*k\[LetterSpace]deg\[LetterSpace]FAAH, vFAAH\[LetterSpace]syn\[LetterSpace]m -> b\[LetterSpace]FAAH\[LetterSpace]MEC*FAAH\[LetterSpace]t*k\[LetterSpace]deg\[LetterSpace]FAAH*MEC, vFAAH\[LetterSpace]syn\[LetterSpace]r -> c\[LetterSpace]FAAH\[LetterSpace]ROB*FAAH\[LetterSpace]t*k\[LetterSpace]deg\[LetterSpace]FAAH, vL\[LetterSpace]UE\[LetterSpace]b -> BRAIN*b\[LetterSpace]FAAH\[LetterSpace]Brain*kcl\[LetterSpace]L*L\[LetterSpace]b[t], vL\[LetterSpace]UE\[LetterSpace]r -> c\[LetterSpace]NAAA\[LetterSpace]ROB*kcl\[LetterSpace]L*L\[LetterSpace]r[t], vL\[LetterSpace]b\[LetterSpace]m -> ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]L*MEC*(L\[LetterSpace]b[t] - L\[LetterSpace]m[t]), vL\[LetterSpace]degr\[LetterSpace]b -> (a\[LetterSpace]FAAH\[LetterSpace]L*BRAIN*kcat\[LetterSpace]FAAH*FAAH\[LetterSpace]b[t]*L\[LetterSpace]b[t])/(FAAH\[LetterSpace]D\[LetterSpace]b*Km\[LetterSpace]FAAH\[LetterSpace]L), vL\[LetterSpace]degr\[LetterSpace]m -> (a\[LetterSpace]FAAH\[LetterSpace]L*kcat\[LetterSpace]FAAH*MEC*FAAH\[LetterSpace]m[t]*L\[LetterSpace]m[t])/(FAAH\[LetterSpace]D\[LetterSpace]m*Km\[LetterSpace]FAAH\[LetterSpace]L), vL\[LetterSpace]degr\[LetterSpace]r -> (a\[LetterSpace]FAAH\[LetterSpace]L*kcat\[LetterSpace]FAAH*ROB*FAAH\[LetterSpace]r[t]*L\[LetterSpace]r[t])/(FAAH\[LetterSpace]D\[LetterSpace]r*Km\[LetterSpace]FAAH\[LetterSpace]L), vL\[LetterSpace]m\[LetterSpace]p -> ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]L*MEC*(L\[LetterSpace]m[t] - Ktr\[LetterSpace]p\[LetterSpace]m\[LetterSpace]L*L\[LetterSpace]p[t]), vL\[LetterSpace]r\[LetterSpace]p -> ktr\[LetterSpace]r\[LetterSpace]p*PLASMA*(-(Ktr\[LetterSpace]p\[LetterSpace]r\[LetterSpace]L*L\[LetterSpace]p[t]) + L\[LetterSpace]r[t]), vL\[LetterSpace]syn\[LetterSpace]b -> (BRAIN*k\[LetterSpace]NL\[LetterSpace]PE*PLD\[LetterSpace]b*NLPE\[LetterSpace]b[t])/(den\[LetterSpace]b*Km\[LetterSpace]NL\[LetterSpace]PE), vL\[LetterSpace]syn\[LetterSpace]r -> (k\[LetterSpace]NL\[LetterSpace]PE*PLD\[LetterSpace]r*ROB*NLPE\[LetterSpace]r[t])/(den\[LetterSpace]r*Km\[LetterSpace]NL\[LetterSpace]PE), vNAPE\[LetterSpace]syn\[LetterSpace]b -> a\[LetterSpace]NAT\[LetterSpace]A*BRAIN*b\[LetterSpace]NAT\[LetterSpace]Brain*p\[LetterSpace]A*Vmax\[LetterSpace]NAT, vNAPE\[LetterSpace]syn\[LetterSpace]r -> a\[LetterSpace]NAT\[LetterSpace]A*c\[LetterSpace]NAT\[LetterSpace]ROB*p\[LetterSpace]A*Vmax\[LetterSpace]NAT, vNLPE\[LetterSpace]syn\[LetterSpace]b -> a\[LetterSpace]NAT\[LetterSpace]L*BRAIN*b\[LetterSpace]NAT\[LetterSpace]Brain*p\[LetterSpace]L*Vmax\[LetterSpace]NAT, vNLPE\[LetterSpace]syn\[LetterSpace]r -> a\[LetterSpace]NAT\[LetterSpace]L*c\[LetterSpace]NAT\[LetterSpace]ROB*p\[LetterSpace]L*Vmax\[LetterSpace]NAT, vNOPE\[LetterSpace]syn\[LetterSpace]b -> a\[LetterSpace]NAT\[LetterSpace]O*BRAIN*b\[LetterSpace]NAT\[LetterSpace]Brain*p\[LetterSpace]O*Vmax\[LetterSpace]NAT, vNOPE\[LetterSpace]syn\[LetterSpace]r -> a\[LetterSpace]NAT\[LetterSpace]O*c\[LetterSpace]NAT\[LetterSpace]ROB*p\[LetterSpace]O*Vmax\[LetterSpace]NAT, vNPPE\[LetterSpace]syn\[LetterSpace]b -> a\[LetterSpace]NAT\[LetterSpace]P*BRAIN*b\[LetterSpace]NAT\[LetterSpace]Brain*p\[LetterSpace]P*Vmax\[LetterSpace]NAT, vNPPE\[LetterSpace]syn\[LetterSpace]r -> a\[LetterSpace]NAT\[LetterSpace]P*c\[LetterSpace]NAT\[LetterSpace]ROB*p\[LetterSpace]P*Vmax\[LetterSpace]NAT, vNSPE\[LetterSpace]syn\[LetterSpace]b -> a\[LetterSpace]NAT\[LetterSpace]S*BRAIN*b\[LetterSpace]NAT\[LetterSpace]Brain*p\[LetterSpace]S*Vmax\[LetterSpace]NAT, vNSPE\[LetterSpace]syn\[LetterSpace]r -> a\[LetterSpace]NAT\[LetterSpace]S*c\[LetterSpace]NAT\[LetterSpace]ROB*p\[LetterSpace]S*Vmax\[LetterSpace]NAT, vO\[LetterSpace]UE\[LetterSpace]b -> BRAIN*b\[LetterSpace]FAAH\[LetterSpace]Brain*kcl\[LetterSpace]O*O\[LetterSpace]b[t], vO\[LetterSpace]UE\[LetterSpace]r -> c\[LetterSpace]NAAA\[LetterSpace]ROB*kcl\[LetterSpace]O*O\[LetterSpace]r[t], vO\[LetterSpace]b\[LetterSpace]m -> ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]O*MEC*(O\[LetterSpace]b[t] - O\[LetterSpace]m[t]), vO\[LetterSpace]degr\[LetterSpace]b -> (a\[LetterSpace]FAAH\[LetterSpace]O*BRAIN*kcat\[LetterSpace]FAAH*FAAH\[LetterSpace]b[t]*O\[LetterSpace]b[t])/(FAAH\[LetterSpace]D\[LetterSpace]b*Km\[LetterSpace]FAAH\[LetterSpace]O), vO\[LetterSpace]degr\[LetterSpace]m -> (a\[LetterSpace]FAAH\[LetterSpace]O*kcat\[LetterSpace]FAAH*MEC*FAAH\[LetterSpace]m[t]*O\[LetterSpace]m[t])/(FAAH\[LetterSpace]D\[LetterSpace]m*Km\[LetterSpace]FAAH\[LetterSpace]O), vO\[LetterSpace]degr\[LetterSpace]r -> (a\[LetterSpace]FAAH\[LetterSpace]O*kcat\[LetterSpace]FAAH*ROB*FAAH\[LetterSpace]r[t]*O\[LetterSpace]r[t])/(FAAH\[LetterSpace]D\[LetterSpace]r*Km\[LetterSpace]FAAH\[LetterSpace]O), vO\[LetterSpace]r\[LetterSpace]p -> ktr\[LetterSpace]r\[LetterSpace]p*PLASMA*(-(Ktr\[LetterSpace]p\[LetterSpace]r\[LetterSpace]O*O\[LetterSpace]p[t]) + O\[LetterSpace]r[t]), vO\[LetterSpace]syn\[LetterSpace]b -> (BRAIN*k\[LetterSpace]NO\[LetterSpace]PE*PLD\[LetterSpace]b*NOPE\[LetterSpace]b[t])/(den\[LetterSpace]b*Km\[LetterSpace]NO\[LetterSpace]PE), vO\[LetterSpace]syn\[LetterSpace]r -> (k\[LetterSpace]NO\[LetterSpace]PE*PLD\[LetterSpace]r*ROB*NOPE\[LetterSpace]r[t])/(den\[LetterSpace]r*Km\[LetterSpace]NO\[LetterSpace]PE), vP\[LetterSpace]UE\[LetterSpace]b -> BRAIN*b\[LetterSpace]FAAH\[LetterSpace]Brain*kcl\[LetterSpace]P*P\[LetterSpace]b[t], vP\[LetterSpace]UE\[LetterSpace]r -> c\[LetterSpace]NAAA\[LetterSpace]ROB*kcl\[LetterSpace]P*P\[LetterSpace]r[t], vP\[LetterSpace]b\[LetterSpace]m -> ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]P*MEC*(P\[LetterSpace]b[t] - P\[LetterSpace]m[t]), vP\[LetterSpace]degr\[LetterSpace]b -> (a\[LetterSpace]FAAH\[LetterSpace]P*BRAIN*kcat\[LetterSpace]FAAH*FAAH\[LetterSpace]b[t]*P\[LetterSpace]b[t])/(FAAH\[LetterSpace]D\[LetterSpace]b*Km\[LetterSpace]FAAH\[LetterSpace]P), vP\[LetterSpace]degr\[LetterSpace]m -> (a\[LetterSpace]FAAH\[LetterSpace]P*kcat\[LetterSpace]FAAH*MEC*FAAH\[LetterSpace]m[t]*P\[LetterSpace]m[t])/(FAAH\[LetterSpace]D\[LetterSpace]m*Km\[LetterSpace]FAAH\[LetterSpace]P), vP\[LetterSpace]degr\[LetterSpace]r -> (a\[LetterSpace]FAAH\[LetterSpace]P*kcat\[LetterSpace]FAAH*ROB*FAAH\[LetterSpace]r[t]*P\[LetterSpace]r[t])/(FAAH\[LetterSpace]D\[LetterSpace]r*Km\[LetterSpace]FAAH\[LetterSpace]P), vP\[LetterSpace]m\[LetterSpace]p -> ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]P*MEC*(P\[LetterSpace]m[t] - Ktr\[LetterSpace]p\[LetterSpace]m\[LetterSpace]P*P\[LetterSpace]p[t]), vP\[LetterSpace]r\[LetterSpace]p -> ktr\[LetterSpace]r\[LetterSpace]p*PLASMA*(-(Ktr\[LetterSpace]p\[LetterSpace]r\[LetterSpace]P*P\[LetterSpace]p[t]) + P\[LetterSpace]r[t]), vP\[LetterSpace]syn\[LetterSpace]b -> (BRAIN*k\[LetterSpace]NP\[LetterSpace]PE*PLD\[LetterSpace]b*NPPE\[LetterSpace]b[t])/(den\[LetterSpace]b*Km\[LetterSpace]NP\[LetterSpace]PE), vP\[LetterSpace]syn\[LetterSpace]r -> (k\[LetterSpace]NP\[LetterSpace]PE*PLD\[LetterSpace]r*ROB*NPPE\[LetterSpace]r[t])/(den\[LetterSpace]r*Km\[LetterSpace]NP\[LetterSpace]PE), vS\[LetterSpace]UE\[LetterSpace]b -> BRAIN*b\[LetterSpace]FAAH\[LetterSpace]Brain*kcl\[LetterSpace]S*S\[LetterSpace]b[t], vS\[LetterSpace]UE\[LetterSpace]r -> c\[LetterSpace]NAAA\[LetterSpace]ROB*kcl\[LetterSpace]S*S\[LetterSpace]r[t], vS\[LetterSpace]b\[LetterSpace]m -> ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]S*MEC*(S\[LetterSpace]b[t] - S\[LetterSpace]m[t]), vS\[LetterSpace]degr\[LetterSpace]b -> (a\[LetterSpace]FAAH\[LetterSpace]S*BRAIN*kcat\[LetterSpace]FAAH*FAAH\[LetterSpace]b[t]*S\[LetterSpace]b[t])/(FAAH\[LetterSpace]D\[LetterSpace]b*Km\[LetterSpace]FAAH\[LetterSpace]S), vS\[LetterSpace]degr\[LetterSpace]m -> (a\[LetterSpace]FAAH\[LetterSpace]S*kcat\[LetterSpace]FAAH*MEC*FAAH\[LetterSpace]m[t]*S\[LetterSpace]m[t])/(FAAH\[LetterSpace]D\[LetterSpace]m*Km\[LetterSpace]FAAH\[LetterSpace]S), vS\[LetterSpace]degr\[LetterSpace]r -> (a\[LetterSpace]FAAH\[LetterSpace]S*kcat\[LetterSpace]FAAH*ROB*FAAH\[LetterSpace]r[t]*S\[LetterSpace]r[t])/(FAAH\[LetterSpace]D\[LetterSpace]r*Km\[LetterSpace]FAAH\[LetterSpace]S), vS\[LetterSpace]m\[LetterSpace]p -> ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]S*MEC*(S\[LetterSpace]m[t] - Ktr\[LetterSpace]p\[LetterSpace]m\[LetterSpace]S*S\[LetterSpace]p[t]), vS\[LetterSpace]r\[LetterSpace]p -> ktr\[LetterSpace]r\[LetterSpace]p*PLASMA*(-(Ktr\[LetterSpace]p\[LetterSpace]r\[LetterSpace]S*S\[LetterSpace]p[t]) + S\[LetterSpace]r[t]), vS\[LetterSpace]syn\[LetterSpace]b -> (BRAIN*k\[LetterSpace]NS\[LetterSpace]PE*PLD\[LetterSpace]b*NSPE\[LetterSpace]b[t])/(den\[LetterSpace]b*Km\[LetterSpace]NS\[LetterSpace]PE), vS\[LetterSpace]syn\[LetterSpace]r -> (k\[LetterSpace]NS\[LetterSpace]PE*PLD\[LetterSpace]r*ROB*NSPE\[LetterSpace]r[t])/(den\[LetterSpace]r*Km\[LetterSpace]NS\[LetterSpace]PE), vo\[LetterSpace]m\[LetterSpace]p -> ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]O*MEC*(O\[LetterSpace]m[t] - Ktr\[LetterSpace]p\[LetterSpace]m\[LetterSpace]O*O\[LetterSpace]p[t])
};

parameters = {
AG2\[LetterSpace]b -> 0.0, Dose -> 10.0, ED50 -> 0.53, Emax\[LetterSpace]PFM -> 0.773, FAAH\[LetterSpace]t -> 78.0, Gut -> 1.65, H -> 0.48, Heart -> 0.31, Kd\[LetterSpace]AG2 -> 3424.0, Kd\[LetterSpace]CB1\[LetterSpace]A -> 239.2, Ki\[LetterSpace]A -> 230.0, Ki\[LetterSpace]L -> 1000.0, Ki\[LetterSpace]O -> 240.0, Ki\[LetterSpace]P -> 6700.0, Ki\[LetterSpace]S -> 840.0, Kidney -> 0.28, Km\[LetterSpace]FAAH\[LetterSpace]A -> 8200.0, Km\[LetterSpace]FAAH\[LetterSpace]L -> 10800.0, Km\[LetterSpace]FAAH\[LetterSpace]O -> 52200.0, Km\[LetterSpace]FAAH\[LetterSpace]P -> 543000.0, Km\[LetterSpace]FAAH\[LetterSpace]S -> 10000.0, Km\[LetterSpace]NA\[LetterSpace]PE -> 2800.0, Km\[LetterSpace]NL\[LetterSpace]PE -> 1000.0, Km\[LetterSpace]NO\[LetterSpace]PE -> 2900.0, Km\[LetterSpace]NP\[LetterSpace]PE -> 3300.0, Km\[LetterSpace]NS\[LetterSpace]PE -> 3400.0, Km\[LetterSpace]PFM -> 26.1, Km\[LetterSpace]p\[LetterSpace]m\[LetterSpace]A -> 1.0, Kp\[LetterSpace]b\[LetterSpace]PF -> 1.3, Kp\[LetterSpace]m\[LetterSpace]PF -> 1.3, Kp\[LetterSpace]r\[LetterSpace]PF -> 1.5, Ktr\[LetterSpace]p\[LetterSpace]m\[LetterSpace]A -> 1.89, Ktr\[LetterSpace]p\[LetterSpace]m\[LetterSpace]L -> 2.77, Ktr\[LetterSpace]p\[LetterSpace]m\[LetterSpace]O -> 9.07, Ktr\[LetterSpace]p\[LetterSpace]m\[LetterSpace]P -> 2.65, Ktr\[LetterSpace]p\[LetterSpace]m\[LetterSpace]S -> 30.01, Ktr\[LetterSpace]p\[LetterSpace]r\[LetterSpace]A -> 0.62, Ktr\[LetterSpace]p\[LetterSpace]r\[LetterSpace]L -> 0.89, Ktr\[LetterSpace]p\[LetterSpace]r\[LetterSpace]O -> 2.8, Ktr\[LetterSpace]p\[LetterSpace]r\[LetterSpace]P -> 0.85, Ktr\[LetterSpace]p\[LetterSpace]r\[LetterSpace]S -> 9.19, LEU -> 0.025, LIVER -> 1.69, Leucocytes -> 0.025, Lungs -> 1.172, M\[LetterSpace]A -> 347.5, M\[LetterSpace]L -> 323.5, M\[LetterSpace]O -> 325.5, M\[LetterSpace]P -> 299.5, M\[LetterSpace]PF -> 455.4, M\[LetterSpace]S -> 321.5, Muscles -> 35.0, PLD\[LetterSpace]b -> 10000000.0, PLD\[LetterSpace]r -> 10000000.0, Pancreas -> 0.077, R\[LetterSpace]PF -> 0.6, Spleen -> 0.192, Testis -> 0.036, Thymus -> 0.029, Vm\[LetterSpace]PFM -> 1511.0, Vmax\[LetterSpace]NAT -> 300.0, Vss\[LetterSpace]PFM -> 58.328, a\[LetterSpace]FAAH\[LetterSpace]A -> 1.0, a\[LetterSpace]FAAH\[LetterSpace]L -> 1.15, a\[LetterSpace]FAAH\[LetterSpace]O -> 5.7, a\[LetterSpace]FAAH\[LetterSpace]P -> 37.8, a\[LetterSpace]FAAH\[LetterSpace]S -> 1.0, a\[LetterSpace]NAT\[LetterSpace]A -> 1.0, a\[LetterSpace]NAT\[LetterSpace]L -> 8.6, a\[LetterSpace]NAT\[LetterSpace]O -> 13.0, a\[LetterSpace]NAT\[LetterSpace]P -> 0.42, a\[LetterSpace]NAT\[LetterSpace]S -> 1.0, b\[LetterSpace]FAAH\[LetterSpace]Brain -> 0.197, b\[LetterSpace]FAAH\[LetterSpace]Gut -> 0.034, b\[LetterSpace]FAAH\[LetterSpace]Kidney -> 0.069, b\[LetterSpace]FAAH\[LetterSpace]Leucocytes -> 0.0, b\[LetterSpace]FAAH\[LetterSpace]Liver -> 1.0, b\[LetterSpace]FAAH\[LetterSpace]Lungs -> 0.032, b\[LetterSpace]FAAH\[LetterSpace]MEC -> 0.137, b\[LetterSpace]FAAH\[LetterSpace]Spleen -> 0.03, b\[LetterSpace]FAAH\[LetterSpace]Testis -> 0.126, b\[LetterSpace]NAAA\[LetterSpace]Brain -> 0.6, b\[LetterSpace]NAAA\[LetterSpace]Gut -> 0.2, b\[LetterSpace]NAAA\[LetterSpace]Heart -> 0.2, b\[LetterSpace]NAAA\[LetterSpace]Kidney -> 0.6, b\[LetterSpace]NAAA\[LetterSpace]Liver -> 1.0, b\[LetterSpace]NAAA\[LetterSpace]Lungs -> 14.0, b\[LetterSpace]NAAA\[LetterSpace]Spleen -> 8.0, b\[LetterSpace]NAAA\[LetterSpace]Testis -> 0.6, b\[LetterSpace]NAAA\[LetterSpace]Thymus -> 4.0, b\[LetterSpace]NAT\[LetterSpace]Brain -> 1.667, b\[LetterSpace]NAT\[LetterSpace]Heart -> 1.0, b\[LetterSpace]NAT\[LetterSpace]Kidney -> 0.667, b\[LetterSpace]NAT\[LetterSpace]Leucocytes -> 0.0, b\[LetterSpace]NAT\[LetterSpace]Lungs -> 0.033, b\[LetterSpace]NAT\[LetterSpace]Muscles -> 0.333, b\[LetterSpace]NAT\[LetterSpace]Pancreas -> 0.333, b\[LetterSpace]NAT\[LetterSpace]Testis -> 0.667, e1 -> 0.0, k\[LetterSpace]NA\[LetterSpace]PE -> 202.0, k\[LetterSpace]NL\[LetterSpace]PE -> 100.0, k\[LetterSpace]NO\[LetterSpace]PE -> 230.0, k\[LetterSpace]NP\[LetterSpace]PE -> 270.0, k\[LetterSpace]NS\[LetterSpace]PE -> 280.0, k\[LetterSpace]deg\[LetterSpace]FAAH -> 0.0051, k\[LetterSpace]inh -> 1.1, kabs\[LetterSpace]PFM -> 2.2, kcat\[LetterSpace]FAAH -> 18000.0, kcl\[LetterSpace]A -> 1.74, kcl\[LetterSpace]L -> 1.25, kcl\[LetterSpace]O -> 2.5, kcl\[LetterSpace]P -> 2.61, kcl\[LetterSpace]S -> 1.2, kin\[LetterSpace]PFM -> 0.117, klinear\[LetterSpace]PFM -> 0.0803, kout\[LetterSpace]PFM -> 0.18, ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]A -> 150.0, ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]L -> 0.0, ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]O -> 10.0, ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]P -> 10.0, ktr\[LetterSpace]m\[LetterSpace]p\[LetterSpace]S -> 10.0, ktr\[LetterSpace]r\[LetterSpace]p -> 100.0, p\[LetterSpace]A -> 0.051, p\[LetterSpace]L -> 0.016, p\[LetterSpace]O -> 0.098, p\[LetterSpace]P -> 0.615, p\[LetterSpace]S -> 0.191, t -> 0.0, BRAIN -> 1.45000004768372, Default -> 1.0, MEC -> 1.49999996210681*^-05, PLASMA -> 2.6489999294281, ROB -> 65.3000030517578
};

assignments = {
tid -> 0.041666666666666664*t, CB1\[LetterSpace]b -> (AG2\[LetterSpace]b/Kd\[LetterSpace]AG2 + A\[LetterSpace]b[t]/Kd\[LetterSpace]CB1\[LetterSpace]A)/(1 + AG2\[LetterSpace]b/Kd\[LetterSpace]AG2 + A\[LetterSpace]b[t]/Kd\[LetterSpace]CB1\[LetterSpace]A), SG\[LetterSpace]p -> (M\[LetterSpace]S*S\[LetterSpace]p[t])/1000, LG\[LetterSpace]p -> (M\[LetterSpace]L*L\[LetterSpace]p[t])/1000, PG\[LetterSpace]p -> (M\[LetterSpace]P*P\[LetterSpace]p[t])/1000, PFG\[LetterSpace]p -> PFM\[LetterSpace]p[t]/(1000*Vss\[LetterSpace]PFM), F\[LetterSpace]r -> (1500.*FAAH\[LetterSpace]r[t])/(FAAHinh\[LetterSpace]r[t] + FAAH\[LetterSpace]r[t]), FAAH\[LetterSpace]D\[LetterSpace]b -> 1 + A\[LetterSpace]b[t]/Km\[LetterSpace]FAAH\[LetterSpace]A + L\[LetterSpace]b[t]/Km\[LetterSpace]FAAH\[LetterSpace]L + O\[LetterSpace]b[t]/Km\[LetterSpace]FAAH\[LetterSpace]O + P\[LetterSpace]b[t]/Km\[LetterSpace]FAAH\[LetterSpace]P + S\[LetterSpace]b[t]/Km\[LetterSpace]FAAH\[LetterSpace]S, PF\[LetterSpace]m -> Kp\[LetterSpace]m\[LetterSpace]PF*PF\[LetterSpace]p, PF\[LetterSpace]r -> Kp\[LetterSpace]r\[LetterSpace]PF*PF\[LetterSpace]p, den\[LetterSpace]b -> 1 + slag1\[LetterSpace]b + slag2\[LetterSpace]b, MD -> 1000000*Dose*F\[LetterSpace]PFM + PFM\[LetterSpace]gut[t], PF\[LetterSpace]b -> Kp\[LetterSpace]b\[LetterSpace]PF*PF\[LetterSpace]p, FAAH\[LetterSpace]D\[LetterSpace]r -> 1 + A\[LetterSpace]r[t]/Km\[LetterSpace]FAAH\[LetterSpace]A + L\[LetterSpace]r[t]/Km\[LetterSpace]FAAH\[LetterSpace]L + O\[LetterSpace]r[t]/Km\[LetterSpace]FAAH\[LetterSpace]O + P\[LetterSpace]r[t]/Km\[LetterSpace]FAAH\[LetterSpace]P + S\[LetterSpace]r[t]/Km\[LetterSpace]FAAH\[LetterSpace]S, den\[LetterSpace]r -> 1 + slag1\[LetterSpace]r + slag2\[LetterSpace]r, c\[LetterSpace]FAAH\[LetterSpace]ROB -> b\[LetterSpace]FAAH\[LetterSpace]Gut*Gut + b\[LetterSpace]FAAH\[LetterSpace]Kidney*Kidney + b\[LetterSpace]FAAH\[LetterSpace]Leucocytes*Leucocytes + b\[LetterSpace]FAAH\[LetterSpace]Liver*LIVER + b\[LetterSpace]FAAH\[LetterSpace]Lungs*Lungs + b\[LetterSpace]FAAH\[LetterSpace]Spleen*Spleen + b\[LetterSpace]FAAH\[LetterSpace]Testis*Testis, AG\[LetterSpace]p -> (M\[LetterSpace]A*A\[LetterSpace]p[t])/1000, OG\[LetterSpace]p -> (M\[LetterSpace]O*O\[LetterSpace]p[t])/1000, FAAH\[LetterSpace]D\[LetterSpace]m -> 1 + A\[LetterSpace]m[t]/Km\[LetterSpace]FAAH\[LetterSpace]A + L\[LetterSpace]m[t]/Km\[LetterSpace]FAAH\[LetterSpace]L + O\[LetterSpace]m[t]/Km\[LetterSpace]FAAH\[LetterSpace]O + P\[LetterSpace]m[t]/Km\[LetterSpace]FAAH\[LetterSpace]P + S\[LetterSpace]m[t]/Km\[LetterSpace]FAAH\[LetterSpace]S, c\[LetterSpace]NAAA\[LetterSpace]ROB -> b\[LetterSpace]NAAA\[LetterSpace]Gut*Gut + b\[LetterSpace]NAAA\[LetterSpace]Heart*Heart + b\[LetterSpace]NAAA\[LetterSpace]Kidney*Kidney + b\[LetterSpace]NAAA\[LetterSpace]Liver*LIVER + b\[LetterSpace]NAAA\[LetterSpace]Lungs*Lungs + b\[LetterSpace]NAAA\[LetterSpace]Spleen*Spleen + b\[LetterSpace]NAAA\[LetterSpace]Thymus*Testis^2 + b\[LetterSpace]NAAA\[LetterSpace]Thymus*Thymus, slag2\[LetterSpace]r -> A\[LetterSpace]r[t]/Ki\[LetterSpace]A + L\[LetterSpace]r[t]/Ki\[LetterSpace]L + O\[LetterSpace]r[t]/Ki\[LetterSpace]O + P\[LetterSpace]r[t]/Ki\[LetterSpace]P + S\[LetterSpace]r[t]/Ki\[LetterSpace]S, slag1\[LetterSpace]r -> NAPE\[LetterSpace]r[t]/Km\[LetterSpace]NA\[LetterSpace]PE + NLPE\[LetterSpace]r[t]/Km\[LetterSpace]NL\[LetterSpace]PE + NOPE\[LetterSpace]r[t]/Km\[LetterSpace]NO\[LetterSpace]PE + NPPE\[LetterSpace]r[t]/Km\[LetterSpace]NP\[LetterSpace]PE + NSPE\[LetterSpace]r[t]/Km\[LetterSpace]NS\[LetterSpace]PE, c\[LetterSpace]NAT\[LetterSpace]ROB -> b\[LetterSpace]NAT\[LetterSpace]Heart*Heart + b\[LetterSpace]NAT\[LetterSpace]Kidney*Kidney + b\[LetterSpace]NAT\[LetterSpace]Leucocytes*Leucocytes + b\[LetterSpace]NAT\[LetterSpace]Lungs*Lungs + b\[LetterSpace]NAT\[LetterSpace]Muscles*Muscles + b\[LetterSpace]NAT\[LetterSpace]Pancreas*Pancreas + b\[LetterSpace]NAT\[LetterSpace]Testis*Testis, slag2\[LetterSpace]b -> A\[LetterSpace]b[t]/Ki\[LetterSpace]A + L\[LetterSpace]b[t]/Ki\[LetterSpace]L + O\[LetterSpace]b[t]/Ki\[LetterSpace]O + P\[LetterSpace]b[t]/Ki\[LetterSpace]P + S\[LetterSpace]b[t]/Ki\[LetterSpace]S, slag1\[LetterSpace]b -> NAPE\[LetterSpace]b[t]/Km\[LetterSpace]NA\[LetterSpace]PE + NLPE\[LetterSpace]b[t]/Km\[LetterSpace]NL\[LetterSpace]PE + NOPE\[LetterSpace]b[t]/Km\[LetterSpace]NO\[LetterSpace]PE + NPPE\[LetterSpace]b[t]/Km\[LetterSpace]NP\[LetterSpace]PE + NSPE\[LetterSpace]b[t]/Km\[LetterSpace]NS\[LetterSpace]PE, PF\[LetterSpace]p -> PFM\[LetterSpace]p[t]/(M\[LetterSpace]PF*Vss\[LetterSpace]PFM), F\[LetterSpace]PFM -> (Dose*Emax\[LetterSpace]PFM)/(Dose + ED50)
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
A\[LetterSpace]b'[t] == 1.0*vA\[LetterSpace]syn\[LetterSpace]b -1.0*vA\[LetterSpace]degr\[LetterSpace]b -1.0*vA\[LetterSpace]UE\[LetterSpace]b -1.0*vA\[LetterSpace]b\[LetterSpace]m, A\[LetterSpace]m'[t] == 1.0*vA\[LetterSpace]b\[LetterSpace]m -1.0*vA\[LetterSpace]degr\[LetterSpace]m -1.0*vA\[LetterSpace]m\[LetterSpace]p, A\[LetterSpace]p'[t] == 1.0*vA\[LetterSpace]m\[LetterSpace]p +1.0*vA\[LetterSpace]r\[LetterSpace]p , A\[LetterSpace]r'[t] == 1.0*vA\[LetterSpace]syn\[LetterSpace]r -1.0*vA\[LetterSpace]degr\[LetterSpace]r -1.0*vA\[LetterSpace]UE\[LetterSpace]r -1.0*vA\[LetterSpace]r\[LetterSpace]p, FAAH\[LetterSpace]b'[t] == 1.0*vFAAH\[LetterSpace]syn\[LetterSpace]b -1.0*vFAAH\[LetterSpace]degr\[LetterSpace]b -1.0*vFAAH\[LetterSpace]inh\[LetterSpace]b, FAAH\[LetterSpace]m'[t] == 1.0*vFAAH\[LetterSpace]syn\[LetterSpace]m -1.0*vFAAH\[LetterSpace]degr\[LetterSpace]m -1.0*vFAAH\[LetterSpace]inh\[LetterSpace]m, FAAH\[LetterSpace]r'[t] == 1.0*vFAAH\[LetterSpace]syn\[LetterSpace]r -1.0*vFAAH\[LetterSpace]degr\[LetterSpace]r -1.0*vFAAH\[LetterSpace]inh\[LetterSpace]r, FAAHinh\[LetterSpace]b'[t] == 1.0*vFAAH\[LetterSpace]inh\[LetterSpace]b -1.0*vFAAH\[LetterSpace]inh\[LetterSpace]degr\[LetterSpace]b, FAAHinh\[LetterSpace]m'[t] == 1.0*vFAAH\[LetterSpace]inh\[LetterSpace]m -1.0*vFAAH\[LetterSpace]inh\[LetterSpace]degr\[LetterSpace]m, FAAHinh\[LetterSpace]r'[t] == 1.0*vFAAH\[LetterSpace]inh\[LetterSpace]r -1.0*vFAAH\[LetterSpace]inh\[LetterSpace]degr\[LetterSpace]r, L\[LetterSpace]b'[t] == 1.0*vL\[LetterSpace]syn\[LetterSpace]b -1.0*vL\[LetterSpace]degr\[LetterSpace]b -1.0*vL\[LetterSpace]UE\[LetterSpace]b -1.0*vL\[LetterSpace]b\[LetterSpace]m, L\[LetterSpace]m'[t] == 1.0*vL\[LetterSpace]b\[LetterSpace]m -1.0*vL\[LetterSpace]degr\[LetterSpace]m -1.0*vL\[LetterSpace]m\[LetterSpace]p, L\[LetterSpace]p'[t] == 1.0*vL\[LetterSpace]m\[LetterSpace]p +1.0*vL\[LetterSpace]r\[LetterSpace]p , L\[LetterSpace]r'[t] == 1.0*vL\[LetterSpace]syn\[LetterSpace]r -1.0*vL\[LetterSpace]degr\[LetterSpace]r -1.0*vL\[LetterSpace]UE\[LetterSpace]r -1.0*vL\[LetterSpace]r\[LetterSpace]p, NAPE\[LetterSpace]b'[t] == 1.0*vNAPE\[LetterSpace]syn\[LetterSpace]b -1.0*vA\[LetterSpace]syn\[LetterSpace]b, NAPE\[LetterSpace]r'[t] == 1.0*vNAPE\[LetterSpace]syn\[LetterSpace]r -1.0*vA\[LetterSpace]syn\[LetterSpace]r, NLPE\[LetterSpace]b'[t] == 1.0*vNLPE\[LetterSpace]syn\[LetterSpace]b -1.0*vL\[LetterSpace]syn\[LetterSpace]b, NLPE\[LetterSpace]r'[t] == 1.0*vNLPE\[LetterSpace]syn\[LetterSpace]r -1.0*vL\[LetterSpace]syn\[LetterSpace]r, NOPE\[LetterSpace]b'[t] == 1.0*vNOPE\[LetterSpace]syn\[LetterSpace]b -1.0*vO\[LetterSpace]syn\[LetterSpace]b, NOPE\[LetterSpace]r'[t] == 1.0*vNOPE\[LetterSpace]syn\[LetterSpace]r -1.0*vO\[LetterSpace]syn\[LetterSpace]r, NPPE\[LetterSpace]b'[t] == 1.0*vNPPE\[LetterSpace]syn\[LetterSpace]b -1.0*vP\[LetterSpace]syn\[LetterSpace]b, NPPE\[LetterSpace]r'[t] == 1.0*vNPPE\[LetterSpace]syn\[LetterSpace]r -1.0*vP\[LetterSpace]syn\[LetterSpace]r, NSPE\[LetterSpace]b'[t] == 1.0*vNSPE\[LetterSpace]syn\[LetterSpace]b -1.0*vS\[LetterSpace]syn\[LetterSpace]b, NSPE\[LetterSpace]r'[t] == 1.0*vNSPE\[LetterSpace]syn\[LetterSpace]r -1.0*vS\[LetterSpace]syn\[LetterSpace]r, O\[LetterSpace]b'[t] == 1.0*vO\[LetterSpace]syn\[LetterSpace]b -1.0*vO\[LetterSpace]degr\[LetterSpace]b -1.0*vO\[LetterSpace]UE\[LetterSpace]b -1.0*vO\[LetterSpace]b\[LetterSpace]m, O\[LetterSpace]m'[t] == 1.0*vO\[LetterSpace]b\[LetterSpace]m -1.0*vO\[LetterSpace]degr\[LetterSpace]m -1.0*vo\[LetterSpace]m\[LetterSpace]p, O\[LetterSpace]p'[t] == 1.0*vo\[LetterSpace]m\[LetterSpace]p +1.0*vO\[LetterSpace]r\[LetterSpace]p , O\[LetterSpace]r'[t] == 1.0*vO\[LetterSpace]syn\[LetterSpace]r -1.0*vO\[LetterSpace]degr\[LetterSpace]r -1.0*vO\[LetterSpace]UE\[LetterSpace]r -1.0*vO\[LetterSpace]r\[LetterSpace]p, PFM\[LetterSpace]gut'[t] ==  -1.0*absorp, PFM\[LetterSpace]p'[t] == 1.0*absorp -1.0*dist -1.0*elim, PFM\[LetterSpace]r'[t] == 1.0*dist , P\[LetterSpace]b'[t] == 1.0*vP\[LetterSpace]syn\[LetterSpace]b -1.0*vP\[LetterSpace]degr\[LetterSpace]b -1.0*vP\[LetterSpace]UE\[LetterSpace]b -1.0*vP\[LetterSpace]b\[LetterSpace]m, P\[LetterSpace]m'[t] == 1.0*vP\[LetterSpace]b\[LetterSpace]m -1.0*vP\[LetterSpace]degr\[LetterSpace]m -1.0*vP\[LetterSpace]m\[LetterSpace]p, P\[LetterSpace]p'[t] == 1.0*vP\[LetterSpace]m\[LetterSpace]p +1.0*vP\[LetterSpace]r\[LetterSpace]p , P\[LetterSpace]r'[t] == 1.0*vP\[LetterSpace]syn\[LetterSpace]r -1.0*vP\[LetterSpace]degr\[LetterSpace]r -1.0*vP\[LetterSpace]UE\[LetterSpace]r -1.0*vP\[LetterSpace]r\[LetterSpace]p, S\[LetterSpace]b'[t] == 1.0*vS\[LetterSpace]syn\[LetterSpace]b -1.0*vS\[LetterSpace]degr\[LetterSpace]b -1.0*vS\[LetterSpace]UE\[LetterSpace]b -1.0*vS\[LetterSpace]b\[LetterSpace]m, S\[LetterSpace]m'[t] == 1.0*vS\[LetterSpace]b\[LetterSpace]m -1.0*vS\[LetterSpace]degr\[LetterSpace]m -1.0*vS\[LetterSpace]m\[LetterSpace]p, S\[LetterSpace]p'[t] == 1.0*vS\[LetterSpace]m\[LetterSpace]p +1.0*vS\[LetterSpace]r\[LetterSpace]p , S\[LetterSpace]r'[t] == 1.0*vS\[LetterSpace]syn\[LetterSpace]r -1.0*vS\[LetterSpace]degr\[LetterSpace]r -1.0*vS\[LetterSpace]UE\[LetterSpace]r -1.0*vS\[LetterSpace]r\[LetterSpace]p
};
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]}]