Model Mathematics

GnRH = \{\begin{matrix} 0.5 if time \geq 0.0 \land time < 4.0 \\ 0.0 if time \geq 4.0 \land time < 24.0 \\ 0.5 if time \geq 24.0 \land time < 28.0 \\ 0.0 if time \geq 28.0 \land time < 148.0 \\ 0.5 if time \geq 148.0 \land time < 152.0 \\ 0.0 if time \geq 152.0 \land time < 157.0 \\ 0.5 if time \geq 157.0 \land time < 161.0 \\ 0.0 if time \geq 161.0 \land time < 281.0 \\ 0.5 if time \geq 281.0 \land time < 285.0 \\ 0.0 if time \geq 285.0 \land time < 295.0 \\ 0.5 if time \geq 295.0 \land time < 299.0 \\ 0.0 if time \geq 299.0 \land time < 419.0 \\ 0.5 if time \geq 419.0 \land time < 424.0 \\ 0.0 if time \geq 424.0 \land time < 464.0 \\ 0.5 if time \geq 464.0 \land time < 468.0 \\ 0.0 if time \geq 468.0 \land time < 588.0 \end{matrix}

\frac{d}{d time} (F) = kdf D - (kfb F GnRH)

\frac{d}{d time} (D) = kbd B - (kdf D)

\frac{d}{d time} (B) = kfb F GnRH - (kbd B)

\frac{d}{d time} (R) = s - ((a1 + a2 B) R)

\frac{d}{d time} (C) = (a1 + a2 B) R