R equation is linearly dependent

[gasilva]

You do not need to solve for R.

See

Requirements.
All parameters must be positive, and S(0)+E(0)+I(0) ≤ 1.

theory

python code

R=1-S-I

The author uses this equation to guarantee consistency.