Model Mathematics

\frac{d}{d time} (GH) = kr1 (\frac{{(\frac{GHRH + ghr_GHRH}{t1})}^{n1}}{{(\frac{GHRH + ghr_GHRH}{t1})}^{n1} + 1.0} + g0 \frac{{(\frac{GHS}{tg0})}^{ng0}}{1.0 + {(\frac{GHS}{tg0})}^{ng0}} + \frac{1.0 + F1_GHS}{1.0 + {(\frac{SRIF_PeV}{t2})}^{n2} + F1_GHS}) - (k1 GH)

\frac{d}{d time} (SRIF_PeV) = (- (k4 SRIF_PeV)) + kr4 \frac{{(\frac{GH}{t5})}^{n5}}{{(\frac{GH}{t5})}^{n5} + 1.0} + S_basal

\frac{d}{d time} (SRIF_ArC) = kr2 \frac{{(\frac{GHRH + ghr_GHRH}{t3})}^{n3}}{1.0 + {(\frac{GHRH + ghr_GHRH}{t3})}^{n3}} - (k2 SRIF_ArC)

\frac{d}{d time} (GHRH) = kr3 \frac{1.0 + F2_GHS}{1.0 + {(\frac{SRIF_PeV + SRIF_ArC}{t4})}^{n4} + F2_GHS} + ghr_GHRH 1.0 - (k3 GHRH)

F1_GHS = g1 \frac{{(\frac{GHS}{tg1})}^{ng1}}{1.0 + {(\frac{GHS}{tg1})}^{ng1}} F2_GHS = g2 \frac{{(\frac{GHS}{tg2})}^{ng2}}{1.0 + {(\frac{GHS}{tg2})}^{ng2}}

\frac{d}{d time} (ghr_GHRH) = inject - (kghr ghr_GHRH) dghr_GHRH_dt = inject - (kghr ghr_GHRH) inject = \{\begin{matrix} 0.0 if time < onset \\ C if time \geq onset \land time ≦ onset + duration \\ 0.0 if time > onset + duration \end{matrix}