∅ > P1
P1 > ∅
∅ > P2
P2 > ∅
∅ > Target
Target > ∅
dsp1p2kd = Kd / 2.0 * (1.0 + (dsp1ksp + P2 * univ) / Kd - pow(pow(1.0 + (dsp1ksp + P2 * univ) / Kd, 2.0) - 4.0 * dsp1ksp * P2 * univ / pow(Kd, 2.0), 0.5))
dsp1ksp = Ksp / 2.0 * (1.0 + (s + P1 * univ) / Ksp - pow(pow(1.0 + (s + P1 * univ) / Ksp, 2.0) - 4.0 * s * P1 * univ / pow(Ksp, 2.0), 0.5))
Note that constraints are not enforced in simulations. It remains the responsibility of the user to verify that simulation results satisfy these constraints.