Model Mathematics

GnRH = \{\begin{matrix} 0.5 if time \geq 0.0 \land time < 0.0666667 \\ 0.0 if time \geq 0.0666667 \land time < 0.4 \\ 0.5 if time \geq 0.4 \land time < 0.4666667 \\ 0.0 if time \geq 0.4666667 \land time < 2.4666667 \\ 0.5 if time \geq 2.4666667 \land time < 2.533333 \\ 0.0 if time \geq 2.533333 \land time < 2.6166667 \\ 0.5 if time \geq 2.6166667 \land time < 2.6833333 \\ 0.0 if time \geq 2.6833333 \land time < 4.6833333 \\ 0.5 if time \geq 4.6833333 \land time < 4.75 \\ 0.0 if time \geq 4.75 \land time < 4.916667 \\ 0.5 if time \geq 4.916667 \land time < 4.983333 \\ 0.0 if time \geq 4.983333 \land time < 6.983333 \\ 0.5 if time \geq 6.983333 \land time < 7.066667 \\ 0.0 if time \geq 7.066667 \land time < 7.733333 \\ 0.5 if time \geq 7.733333 \land time < 7.8 \\ 0.0 if time \geq 7.8 \land time < 9.8 \end{matrix}

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

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

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

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