@@ -16,11 +16,11 @@ j₅ = ConstantRateJump(k1 * I, [R ~ R + 1])
1616j₆ = VariableRateJump (k1 * k2 / (1 + t) * S, [S ~ S - 1 , R ~ R + 1 ])
1717eqs = [j₁, j₂, j₃, j₄, j₅, j₆]
1818@named js = JumpSystem (eqs, t, [S, I, R], [k1, k2])
19- S = value (S);
20- I = value (I);
21- R = value (R);
22- k1 = value (k1);
23- k2 = value (k2);
19+ S = value (S)
20+ I = value (I)
21+ R = value (R)
22+ k1 = value (k1)
23+ k2 = value (k2)
2424# eq to vars they depend on
2525eq_sdeps = [Variable[], [S], [S, I], [S, R], [I], [S]]
2626eq_sidepsf = [Int[], [1 ], [1 , 2 ], [1 , 3 ], [2 ], [1 ]]
8484 j₆ = VariableRateJump (k1 * k2 / (1 + t) * S, [S ~ S - 1 , R ~ R + 1 ])
8585 eqs = [j₁, j₂, j₃, j₄, j₅, j₆]
8686 @named js = JumpSystem (eqs, t, [S, I, R], [k1, k2])
87- S = value (S);
88- I = value (I);
89- R = value (R);
90- k1 = value (k1);
91- k2 = value (k2);
87+ S = value (S)
88+ I = value (I)
89+ R = value (R)
90+ k1 = value (k1)
91+ k2 = value (k2)
9292 # eq to vars they depend on
9393 eq_sdeps = [Variable[], [S], [S, I], [S, R], [I]]
9494 eq_sidepsf = [Int[], [1 ], [1 , 2 ], [1 , 3 ], [2 ]]
0 commit comments