@@ -16,11 +16,11 @@ j₅ = ConstantRateJump(k1 * I, [R ~ R + 1])
16
16
j₆ = VariableRateJump (k1 * k2 / (1 + t) * S, [S ~ S - 1 , R ~ R + 1 ])
17
17
eqs = [j₁, j₂, j₃, j₄, j₅, j₆]
18
18
@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)
24
24
# eq to vars they depend on
25
25
eq_sdeps = [Variable[], [S], [S, I], [S, R], [I], [S]]
26
26
eq_sidepsf = [Int[], [1 ], [1 , 2 ], [1 , 3 ], [2 ], [1 ]]
84
84
j₆ = VariableRateJump (k1 * k2 / (1 + t) * S, [S ~ S - 1 , R ~ R + 1 ])
85
85
eqs = [j₁, j₂, j₃, j₄, j₅, j₆]
86
86
@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)
92
92
# eq to vars they depend on
93
93
eq_sdeps = [Variable[], [S], [S, I], [S, R], [I]]
94
94
eq_sidepsf = [Int[], [1 ], [1 , 2 ], [1 , 3 ], [2 ]]
0 commit comments