\frac{dx_{1}}{dt} = \left(1 \cdot k_{21} \cdot 2060 / 365 + -1 \cdot k_{21} \cdot 1 / \left(2 \cdot 365\right) \cdot x_{1} + -1 \cdot k_{21} \cdot 45685 / \left(439197 \cdot 30\right) \cdot x_{1} \cdot x_{2}\right) / k_{21}\\
\frac{dx_{2}}{dt} = \left(1 \cdot k_{21} \cdot 45685 / \left(439197 \cdot 30\right) \cdot x_{1} \cdot x_{2} + -1 \cdot k_{21} \cdot 45682 / \left(439197 \cdot 30\right) \cdot x_{2} + -1 \cdot k_{21} \cdot 1 / \left(2 \cdot 365\right) \cdot x_{2}\right) / k_{21}\\
\frac{dx_{3}}{dt} = \left(1 \cdot k_{22} \cdot \frac{237}{10} \cdot k_{20} / 365000 + -1 \cdot k_{22} \cdot 1 / \left(\frac{346}{5} \cdot 365\right) \cdot x_{3} + -1 \cdot 16 / \left(108 \cdot 1000000 \cdot \frac{1825}{2}\right) \cdot x_{3} \cdot x_{2}\right) / k_{22}\\
\frac{dx_{4}}{dt} = \left(1 \cdot 16 / \left(108 \cdot 1000000 \cdot \frac{1825}{2}\right) \cdot x_{3} \cdot x_{2} + -1 \cdot k_{22} \cdot 6 / \left(108 \cdot 1000000 \cdot \frac{1825}{2}\right) \cdot x_{4} + -1 \cdot k_{22} \cdot 1 / \left(\frac{346}{5} \cdot 365\right) \cdot x_{4}\right) / k_{22}