J1
w > IFNb_mRNA
J10
w > TNFenv
J11
w > STATm
J12
w > STAT
J2
w > IFNb_env
J3
w > STATP2n
J4
w > SOCS1m
J5
w > IRF7m
J6
w > IRF7Pn
J7
w > IFNa_mRNA
J8
w > IFNa_env
J9
w > TNFam
NS = r5 * pow(time, n3) / (pow(bm, n3) + pow(time, n3))
IC1 = (1.0 + sp * pow(NS / delta1, n1)) / (1.0 + pow(NS / delta1, n1))
IC2 = (1.0 + sv * pow(NS / delta2, n2)) / (1.0 + pow(NS / delta2, n2))
TJ = TJtot * (IFNb_env + IFNa_env) / (K3 + IFNb_env + IFNa_env) * 1.0 / (1.0 + K9 * SOCS1m / delta)
IC2ifa = (1.0 + 3.0 * sv * pow(NS / delta2, n2)) / (1.0 + pow(NS / delta2, n2))
Note that constraints are not enforced in simulations. It remains the responsibility of the user to verify that simulation results satisfy these constraints.