begin name
roblitz1
end name

begin reactions
re100 {1.0}Agod = {1.0}$s107
re101 {1.0}$s108 = {1.0}Agoc
re102 {1.0}Antd = {1.0}s113
re103 {1.0}s114 = {1.0}Antc
re104 {1.0}Antc = {1.0}s115
re105 {1.0}AntR = {1.0}RGnRHa + {1.0}Antc
re106 {1.0}RGnRHa + {1.0}Antc = {1.0}AntR
re107 {1.0}AntR = {1.0}s116
re108 {1.0}Antc = {1.0}Antp
re109 {1.0}Antp = {1.0}Antc
re11 {1.0}RGnRHi = {1.0}$sa1degraded
re15 {1.0}GnRHRa = {1.0}GnRH + {1.0}RGnRHa
re2 {1.0}GnRHRa = {1.0}GnRHRi
re24 {1.0}GnRH = {1.0}$sa3degraded
re25 {1.0}$s33 = {1.0}LHPit
re26 {1.0}$s38 = {1.0}FSHpit
re28 {1.0}LHPit = {1.0}$s92
re29 {1.0}FSHpit = {1.0}$s94
re3 {1.0}GnRHRi = {1.0}GnRHRa
re35 {1.0}LHbld + {1.0}RLH = {1.0}LHR
re36 {1.0}LHR = {1.0}RLHdes
re37 {1.0}RLHdes = {1.0}RLH
re38 {1.0}LHbld = {1.0}$sa28degraded
re39 {1.0}FSHR = {1.0}RFSHdes
re4 {1.0}RGnRHi = {1.0}RGnRHa
re40 {1.0}RFSHdes = {1.0}RFSH
re42 {1.0}AF1 = {1.0}AF2
re43 {1.0}$s62 = {1.0}AF1
re44 {1.0}AF2 = {1.0}AF3
re45 {1.0}AF3 = {1.0}AF4
re46 {1.0}AF4 = {1.0}PrF
re49 {1.0}Sc1 = {1.0}Sc2
re5 {1.0}RGnRHa = {1.0}RGnRHi
re50 {1.0}Sc2 = {1.0}Lut1
re51 {1.0}Lut1 = {1.0}Lut2
re52 {1.0}Lut2 = {1.0}Lut3
re53 {1.0}Lut3 = {1.0}Lut4
re54 {1.0}$s64 = {1.0}RFoll
re56 {1.0}RFoll = {1.0}$sa61degraded
re57 {1.0}$s66 = {1.0}AF3
re58 {1.0}$s67 = {1.0}AF4
re59 {1.0}PrF = {1.0}$sa52degraded
re6 {1.0}GnRHRi = {1.0}RGnRHi
re60 {1.0}$s69 = {1.0}OvF
re61 {1.0}OvF = {1.0}$sa53degraded
re62 {1.0}$s71 = {1.0}Sc1
re64 {1.0}Lut4 = {1.0}$s72
re65 {1.0}$s74 = {1.0}E2
re66 {1.0}E2 = {1.0}$sa75degraded
re67 {1.0}$s76 = {1.0}P4
re69 {1.0}$s78 = {1.0}InhA
re71 {1.0}P4 = {1.0}$sa78degraded
re72 {1.0}$s82 = {1.0}InhB
re73 {1.0}InhB = {1.0}$sa86degraded
re74 {1.0}InhAdelay = {1.0}$sa35degraded
re75 {1.0}$s85 = {1.0}RGnRHi
re76 {1.0}GnRHRi = {1.0}$csa1degraded
re78 {1.0}$s87 = {1.0}GnRH
re8 {1.0}GnRH + {1.0}RGnRHa = {1.0}GnRHRa
re82 {1.0}FSHbld = {1.0}$sa31degraded
re83 {1.0}$s93 = {1.0}LHbld
re84 {1.0}$s95 = {1.0}FSHbld
re85 {1.0}FSHbld + {1.0}RFSH = {1.0}FSHR
re87 {1.0}InhA = {1.0}InhAdelay
re90 {1.0}Agoc = {1.0}$s102
re93 {1.0}AgoRa = {1.0}AgoRi
re94 {1.0}AgoRi = {1.0}AgoRa
re95 {1.0}AgoRi = {1.0}$s106
re97 {1.0}RGnRHa + {1.0}Agoc = {1.0}AgoRa
re98 {1.0}AgoRa = {1.0}Agoc + {1.0}RGnRHa
re99 {1.0}AgoRi = {1.0}RGnRHi
end reactions

begin rate equations
re100 = default*function62[Agod[t], default, p274]
re101 = default*function63[Agod[t], default, p273, p274]
re102 = default*function64[Antd[t], default, p474]
re103 = default*function65[Antd[t], default, p473, p474]
re104 = default*function66[Antc[t], default, p475]
re105 = default*function67[AntR[t], default, p513]
re106 = default*function68[Antc[t], RGnRHa[t], default, p512]
re107 = default*function69[AntR[t], default, p514]
re108 = default*function70[Antc[t], default, p476]
re109 = default*function71[Antp[t], default, p477]
re11 = default*function7[RGnRHi[t], default, p308]
re15 = default*function8[GnRHRa[t], default, p303]
re2 = default*function1[GnRHRa[t], default, p309]
re24 = default*function9[GnRH[t], default, p300]
re25 = default*function10[E2[t], P4[t], default, facE2, facP4, p1, p2, p3, p4, p6, p7]
re26 = default*function11[InhAdelay[t], InhB[t], default, freq, p11, p13, p21, p22, p23, p24, p25]
re28 = default*function12[AgoRa[t], GnRHRa[t], LHPit[t], default, p16, p5, p8, p9]
re29 = default*function13[AgoRa[t], FSHpit[t], GnRHRa[t], default, p17, p18, p20, p28]
re3 = default*function2[GnRHRi[t], default, p310]
re35 = default*function14[LHbld[t], RLH[t], default, facLH, p230]
re36 = default*function15[LHR[t], default, p234]
re37 = default*function16[RLHdes[t], default, p232]
re38 = default*function17[LHbld[t], default, p231]
re39 = default*function18[FSHR[t], default, p244]
re4 = default*function3[RGnRHi[t], default, p307]
re40 = default*function19[RFSHdes[t], default, p242]
re42 = default*function20[AF1[t], FSHR[t], default, p50]
re43 = default*function21[FSHR[t], default, p47, p48, p49]
re44 = default*function22[AF2[t], LHR[t], RFoll[t], default, p46, p51, p52]
re45 = default*function23[AF3[t], LHR[t], RFoll[t], default, p32, p43, p52]
re46 = default*function24[AF4[t], LHR[t], RFoll[t], default, p34, p52]
re49 = default*function25[Sc1[t], default, p37]
re5 = default*function4[RGnRHa[t], default, p306]
re50 = default*function26[Sc2[t], default, p38]
re51 = default*function27[AgoRa[t], GnRHRa[t], Lut1[t], default, p39, p80, p83, p84]
re52 = default*function28[AgoRa[t], GnRHRa[t], Lut2[t], default, p40, p80, p83, p84]
re53 = default*function29[AgoRa[t], GnRHRa[t], Lut3[t], default, p41, p80, p83, p84]
re54 = default*function30[FSHbld[t], default, p90, p91, p94]
re56 = default*function31[P4[t], RFoll[t], default, facP4, p92, p93, p95]
re57 = default*function32[AF3[t], FSHR[t], default, p31, p55]
re58 = default*function33[AF4[t], LHR[t], default, p33, p44, p52, p55]
re59 = default*function34[LHR[t], PrF[t], RFoll[t], default, p35, p45, p52]
re6 = default*function5[GnRHRi[t], default, p305]
re60 = default*function35[LHR[t], PrF[t], RFoll[t], default, p27, p45, p52, p53, p54]
re61 = default*function36[OvF[t], default, p36]
re62 = default*function37[OvF[t], default, p26, p56, p57]
re64 = default*function38[AgoRa[t], GnRHRa[t], Lut4[t], default, p42, p80, p83, p84]
re65 = default*function39[AF2[t], AF3[t], AF4[t], LHbld[t], Lut1[t], Lut4[t], PrF[t], default, facE2, p152, p158, p159, p160, p161, p164, p165]
re66 = default*function40[E2[t], default, p154]
re67 = default*function41[Lut4[t], default, facP4, p166, p167]
re69 = default*function42[Lut1[t], Lut2[t], Lut3[t], Lut4[t], PrF[t], Sc1[t], default, p168, p169, p170, p171, p172, p177, p178]
re71 = default*function43[P4[t], default, p155]
re72 = default*function44[AF2[t], Sc2[t], default, p173, p174, p175]
re73 = default*function45[InhB[t], default, p157]
re74 = default*function46[InhAdelay[t], default, p30]
re75 = default*function47[default, p311]
re76 = default*function48[GnRHRi[t], default, p304]
re78 = default*function49[default, freq, mass]
re8 = default*function6[GnRH[t], RGnRHa[t], default, p302]
re82 = default*function50[FSHbld[t], default, p241]
re83 = default*function51[AgoRa[t], GnRHRa[t], LHPit[t], default, facLH, p12, p16, p5, p8, p9]
re84 = default*function52[AgoRa[t], FSHpit[t], GnRHRa[t], default, facFSH, p12, p17, p18, p20, p28]
re85 = default*function53[FSHbld[t], RFSH[t], default, facFSH, p240]
re87 = default*function54[InhA[t], default, p156]
re90 = default*function55[Agoc[t], default, p275]
re93 = default*function56[AgoRa[t], default, p319]
re94 = default*function57[AgoRi[t], default, p320]
re95 = default*function58[AgoRi[t], default, p314]
re97 = default*function59[Agoc[t], RGnRHa[t], default, p312]
re98 = default*function60[AgoRa[t], default, p313]
re99 = default*function61[AgoRi[t], default, p315]
end rate equations

begin parameters
facE2 = 1.0
facFSH = 1.0
facLH = 1.0
facP4 = 1.0
p1 = 7309.91587614104
p11 = 10.0
p12 = 5.0
p13 = 3.0
p152 = 2.09450510112762
p154 = 5.23500984428137
p155 = 5.12958654018257
p156 = 4.28718453552891
p157 = 172.453910864507
p158 = 51.558081260068
p159 = 9.28
p16 = 0.00476024700196886
p160 = 3480.27
p161 = 0.972
p164 = 1713.71039914086
p165 = 8675.13871487382
p166 = 0.9426346876678
p167 = 761.643100053696
p168 = 1.44522999821013
p169 = 2.28494719885448
p17 = 0.0569894397708967
p170 = 28.2110255951316
p171 = 216.85
p172 = 114.247359942724
p173 = 89.4934669769107
p174 = 447.467334884553
p175 = 134240.200465366
p177 = 60.0
p178 = 180.0
p18 = 0.0003
p2 = 7309.91587614104
p20 = 2.0
p201 = 1.0
p202 = 16.0
p203 = 1.2
p204 = 2.0
p205 = 220.0
p206 = 10.0
p208 = 220.0
p209 = 2.0
p21 = 22129.0495793807
p210 = 9.6
p211 = 1.0
p22 = 95.81
p23 = 70.0
p230 = 2.1430105602291
p231 = 74.8505459101486
p232 = 68.9493466976911
p234 = 183.363164488992
p24 = 5.0
p240 = 3.52890638983354
p241 = 114.247359942724
p242 = 61.0291748702345
p244 = 138.303203866118
p25 = 2.0
p26 = 1.20834079112225
p269 = 91.0
p27 = 7.98433864327904
p272 = 100.0
p273 = 38.12
p274 = 54.2
p275 = 2.65
p28 = 0.27201539287632
p3 = 192.2041
p30 = 0.1989
p300 = 0.447467334884553
p301 = 0.005593
p302 = 322.176481116879
p303 = 644.352962233757
p304 = 0.00894934669769107
p305 = 32.2176481116879
p306 = 3.22176481116878
p307 = 32.2176481116879
p308 = 0.0894934669769107
p309 = 32.2176481116879
p31 = 0.122060139609808
p310 = 3.222
p311 = "8.94934669769107e-05"
p312 = 322.18
p313 = 644.35
p314 = 0.009
p315 = 32.22
p319 = 32.22
p32 = 122.060139609808
p320 = 3.22
p33 = 12.2060139609808
p34 = 332.754608913549
p35 = 122.060139609808
p36 = 12.2060139609808
p37 = 1.22060139609808
p38 = 0.958117057454806
p39 = 0.924735994272418
p4 = 2.3708
p40 = 0.756765706103455
p41 = 0.610291748702345
p42 = 0.542840522641847
p43 = 5.0
p44 = 2.0
p45 = 6.0
p46 = 3.68896
p469 = 34.0
p47 = 3.0
p472 = 500.0
p473 = 34.9
p474 = 45.56
p475 = 5.0
p476 = 3.216
p477 = 4.76
p48 = 0.608121
p49 = 3.66180418829425
p5 = 0.190415249686773
p50 = 1.22060139609808
p51 = 4.88231609092536
p512 = 322.18
p513 = 644.35
p514 = 0.009
p52 = 2.7262
p53 = 3.0
p54 = 10.0
p55 = 10.0
p56 = 0.02
p57 = 10.0
p6 = 10.0
p7 = 1.0
p8 = 0.0003
p80 = 20.0
p83 = 0.0008
p84 = 5.0
p9 = 5.0
p90 = 3.0
p91 = 5.0
p92 = 1.2348
p93 = 5.0
p94 = 0.2186056917845
p95 = 1.34267048505459
csa1degraded = 1.0
s102 = 1.0
s106 = 1.0
s107 = 1.0
s108 = 1.0
s33 = 1.0
s38 = 1.0
s62 = 1.0
s64 = 1.0
s66 = 1.0
s67 = 1.0
s69 = 1.0
s71 = 1.0
s72 = 1.0
s74 = 1.0
s76 = 1.0
s78 = 1.0
s82 = 1.0
s85 = 1.0
s87 = 1.0
s92 = 1.0
s93 = 1.0
s94 = 1.0
s95 = 1.0
sa1degraded = 1.0
sa28degraded = 1.0
sa31degraded = 1.0
sa35degraded = 1.0
sa3degraded = 1.0
sa52degraded = 1.0
sa53degraded = 1.0
sa61degraded = 1.0
sa75degraded = 1.0
sa78degraded = 1.0
sa86degraded = 1.0
default = 1.0
end parameters

begin initial conditions
AF1[0] = AF1i
AF2[0] = AF2i
AF3[0] = AF3i
AF4[0] = AF4i
AgoRa[0] = AgoRai
AgoRi[0] = AgoRii
Agoc[0] = Agoci
Agod[0] = Agodi
AntR[0] = AntRi
Antc[0] = Antci
Antd[0] = Antdi
Antp[0] = Antpi
E2[0] = E2i
FSHR[0] = FSHRi
FSHbld[0] = FSHbldi
FSHpit[0] = FSHpiti
GnRH[0] = GnRHi
GnRHRa[0] = GnRHRai
GnRHRi[0] = GnRHRii
InhA[0] = InhAi
InhAdelay[0] = InhAdelayi
InhB[0] = InhBi
LHPit[0] = LHPiti
LHR[0] = LHRi
LHbld[0] = LHbldi
Lut1[0] = Lut1i
Lut2[0] = Lut2i
Lut3[0] = Lut3i
Lut4[0] = Lut4i
OvF[0] = OvFi
P4[0] = P4i
PrF[0] = PrFi
RFSH[0] = RFSHi
RFSHdes[0] = RFSHdesi
RFoll[0] = RFolli
RGnRHa[0] = RGnRHai
RGnRHi[0] = RGnRHii
RLH[0] = RLHi
RLHdes[0] = RLHdesi
Sc1[0] = Sc1i
Sc2[0] = Sc2i
s113[0] = s113i
s114[0] = s114i
s115[0] = s115i
s116[0] = s116i
end initial conditions

begin initial values
AF1i = 2.26
AF2i = 19.92
AF3i = 0.504
AF4i = 1.604e-05
AgoRai = 0.0
AgoRii = 0.0
Agoci = 0.0
Agodi = 0.0
AntRi = 0.0
Antci = 0.0
Antdi = 0.0
Antpi = 0.0
E2i = 45.0
FSHRi = 0.662
FSHbldi = 4.1
FSHpiti = 48627.2
GnRHi = 0.00326
GnRHRai = 1.447e-05
GnRHRii = 1.294e-05
InhAi = 1.07
InhAdelayi = 86.84
InhBi = 52.18
LHPiti = 261119.0
LHRi = 0.263
LHbldi = 2.667
Lut1i = 2.762e-05
Lut2i = 0.0003772
Lut3i = 0.00334
Lut4i = 0.01385
OvFi = 7.652e-19
P4i = 2.41
PrFi = 0.244
RFSHi = 6.341
RFSHdesi = 1.497
RFolli = 0.119
RGnRHai = 0.00928
RGnRHii = 0.0009409
RLHi = 8.41
RLHdesi = 0.699
Sc1i = 1.098e-08
Sc2i = 2.171e-06
s113i = 0.0
s114i = 0.0
s115i = 0.0
s116i = 0.0
end initial values

begin assignment rules
mass := p301*(((default*p201*E2[t])/(facE2*p208*p7))^p209/(1 + ((default*p201*E2[t])/(facE2*p208*p7))^p209) + (1 + ((default*p201*E2[t])/(facE2*p210*p7))^p211)^(-1))
freq := (p202*(1 + (p201*((default*p201*E2[t])/(facE2*p205*p7))^p206)/(1 + ((default*p201*E2[t])/(facE2*p205*p7))^p206)))/(1 + ((default*p201*P4[t])/(facP4*p203*p7))^p204)
end assignment rules

begin function definitions
function44[AF2_,Sc2_,default_,p173_,p174_,p175_] := (p173 + AF2*default*p174 + default*p175*Sc2)/default
function45[InhB_,default_,p157_] := InhB*p157
function46[InhAdelay_,default_,p30_] := InhAdelay*p30
function47[default_,p311_] := p311/default
function49[default_,freq_,mass_] := (freq*mass)/default
function43[P4_,default_,p155_] := p155*P4
function53[FSHbld_,RFSH_,default_,facFSH_,p240_] := (default*FSHbld*p240*RFSH)/facFSH
function54[InhA_,default_,p156_] := InhA*p156
function55[Agoc_,default_,p275_] := Agoc*p275
function57[AgoRi_,default_,p320_] := AgoRi*p320
function58[AgoRi_,default_,p314_] := AgoRi*p314
function59[Agoc_,RGnRHa_,default_,p312_] := Agoc*default*p312*RGnRHa
function62[Agod_,default_,p274_] := Agod*p274
function63[Agod_,default_,p273_,p274_] := (Agod*p274)/p273
function64[Antd_,default_,p474_] := (Antd*p474)/default
function65[Antd_,default_,p473_,p474_] := (Antd*p474)/(default*p473)
function66[Antc_,default_,p475_] := (Antc*p475)/default
function67[AntR_,default_,p513_] := (AntR*p513)/default
function69[AntR_,default_,p514_] := (AntR*p514)/default
function71[Antp_,default_,p477_] := (Antp*p477)/default
function70[Antc_,default_,p476_] := (Antc*p476)/default
function68[Antc_,RGnRHa_,default_,p512_] := Antc*p512*RGnRHa
function61[AgoRi_,default_,p315_] := AgoRi*p315
function60[AgoRa_,default_,p313_] := AgoRa*p313
function56[AgoRa_,default_,p319_] := AgoRa*p319
function52[AgoRa_,FSHpit_,GnRHRa_,default_,facFSH_,p12_,p17_,p18_,p20_,p28_] := (facFSH*FSHpit*(p17 + (((AgoRa*default + default*GnRHRa)/p18)^p20*p28)/(1 + ((AgoRa*default + default*GnRHRa)/p18)^p20)))/p12
function51[AgoRa_,GnRHRa_,LHPit_,default_,facLH_,p12_,p16_,p5_,p8_,p9_] := (facLH*LHPit*(p16 + (p5*((AgoRa*default + default*GnRHRa)/p8)^p9)/(1 + ((AgoRa*default + default*GnRHRa)/p8)^p9)))/p12
function50[FSHbld_,default_,p241_] := FSHbld*p241
function48[GnRHRi_,default_,p304_] := GnRHRi*p304
function42[Lut1_,Lut2_,Lut3_,Lut4_,PrF_,Sc1_,default_,p168_,p169_,p170_,p171_,p172_,p177_,p178_] := (p168 + default*Lut2*p170 + default*Lut3*p171 + default*Lut4*p172 + default*Lut1*p178 + default*p169*PrF + default*p177*Sc1)/default
function38[AgoRa_,GnRHRa_,Lut4_,default_,p42_,p80_,p83_,p84_] := Lut4*p42*(1 + (p80*((AgoRa*default + default*GnRHRa)/p83)^p84)/(1 + ((AgoRa*default + default*GnRHRa)/p83)^p84))
function31[P4_,RFoll_,default_,facP4_,p92_,p93_,p95_] := (((default*P4)/(facP4*p92))^p93*p95*RFoll)/(1 + ((default*P4)/(facP4*p92))^p93)
function29[AgoRa_,GnRHRa_,Lut3_,default_,p41_,p80_,p83_,p84_] := Lut3*p41*(1 + (p80*((AgoRa*default + default*GnRHRa)/p83)^p84)/(1 + ((AgoRa*default + default*GnRHRa)/p83)^p84))
function27[AgoRa_,GnRHRa_,Lut1_,default_,p39_,p80_,p83_,p84_] := Lut1*p39*(1 + (p80*((AgoRa*default + default*GnRHRa)/p83)^p84)/(1 + ((AgoRa*default + default*GnRHRa)/p83)^p84))
function21[FSHR_,default_,p47_,p48_,p49_] := (((default*FSHR)/p48)^p47*p49)/(default*(1 + ((default*FSHR)/p48)^p47))
function13[AgoRa_,FSHpit_,GnRHRa_,default_,p17_,p18_,p20_,p28_] := FSHpit*(p17 + (((AgoRa*default + default*GnRHRa)/p18)^p20*p28)/(1 + ((AgoRa*default + default*GnRHRa)/p18)^p20))
function12[AgoRa_,GnRHRa_,LHPit_,default_,p16_,p5_,p8_,p9_] := LHPit*(p16 + (p5*((AgoRa*default + default*GnRHRa)/p8)^p9)/(1 + ((AgoRa*default + default*GnRHRa)/p8)^p9))
function1[GnRHRa_,default_,p309_] := GnRHRa*p309
function2[GnRHRi_,default_,p310_] := GnRHRi*p310
function3[RGnRHi_,default_,p307_] := p307*RGnRHi
function4[RGnRHa_,default_,p306_] := p306*RGnRHa
function5[GnRHRi_,default_,p305_] := GnRHRi*p305
function6[GnRH_,RGnRHa_,default_,p302_] := default*GnRH*p302*RGnRHa
function7[RGnRHi_,default_,p308_] := p308*RGnRHi
function8[GnRHRa_,default_,p303_] := GnRHRa*p303
function9[GnRH_,default_,p300_] := GnRH*p300
function10[E2_,P4_,default_,facE2_,facP4_,p1_,p2_,p3_,p4_,p6_,p7_] := (p1 + (p2*((default*E2)/(facE2*p3))^p6)/(1 + ((default*E2)/(facE2*p3))^p6))/(default*(1 + ((default*P4)/(facP4*p4))^p7))
function11[InhAdelay_,InhB_,default_,freq_,p11_,p13_,p21_,p22_,p23_,p24_,p25_] := p21/(default*(1 + (freq/p11)^p13)*(1 + ((default*InhAdelay)/p22)^p24 + ((default*InhB)/p23)^p25))
function14[LHbld_,RLH_,default_,facLH_,p230_] := (default*LHbld*p230*RLH)/facLH
function15[LHR_,default_,p234_] := LHR*p234
function16[RLHdes_,default_,p232_] := p232*RLHdes
function17[LHbld_,default_,p231_] := LHbld*p231
function18[FSHR_,default_,p244_] := FSHR*p244
function19[RFSHdes_,default_,p242_] := p242*RFSHdes
function20[AF1_,FSHR_,default_,p50_] := AF1*default*FSHR*p50
function25[Sc1_,default_,p37_] := p37*Sc1
function26[Sc2_,default_,p38_] := p38*Sc2
function22[AF2_,LHR_,RFoll_,default_,p46_,p51_,p52_] := AF2*default*p51*((default*LHR)/p52)^p46*RFoll
function23[AF3_,LHR_,RFoll_,default_,p32_,p43_,p52_] := AF3*default*p32*((default*LHR)/p52)^p43*RFoll
function24[AF4_,LHR_,RFoll_,default_,p34_,p52_] := (AF4*default^2*LHR*p34*RFoll)/p52
function28[AgoRa_,GnRHRa_,Lut2_,default_,p40_,p80_,p83_,p84_] := Lut2*p40*(1 + (p80*((AgoRa*default + default*GnRHRa)/p83)^p84)/(1 + ((AgoRa*default + default*GnRHRa)/p83)^p84))
function30[FSHbld_,default_,p90_,p91_,p94_] := (((default*FSHbld)/p90)^p91*p94)/(default*(1 + ((default*FSHbld)/p90)^p91))
function32[AF3_,FSHR_,default_,p31_,p55_] := AF3*default*FSHR*p31*(1 - (AF3*default)/p55)
function33[AF4_,LHR_,default_,p33_,p44_,p52_,p55_] := AF4*p33*((default*LHR)/p52)^p44*(1 - (AF4*default)/p55)
function34[LHR_,PrF_,RFoll_,default_,p35_,p45_,p52_] := default*p35*((default*LHR)/p52)^p45*PrF*RFoll
function35[LHR_,PrF_,RFoll_,default_,p27_,p45_,p52_,p53_,p54_] := (p27*((default*LHR)/p52)^p45*((default*PrF)/p53)^p54*RFoll)/(1 + ((default*PrF)/p53)^p54)
function36[OvF_,default_,p36_] := OvF*p36
function37[OvF_,default_,p26_,p56_,p57_] := (p26*((default*OvF)/p56)^p57)/(default*(1 + ((default*OvF)/p56)^p57))
function39[AF2_,AF3_,AF4_,LHbld_,Lut1_,Lut4_,PrF_,default_,facE2_,p152_,p158_,p159_,p160_,p161_,p164_,p165_] := (facE2*(AF2*default*p152 + p158 + AF3*default^2*LHbld*p159 + AF4*default*p160 + default*Lut1*p164 + default*Lut4*p165 + default^2*LHbld*p161*PrF))/default
function40[E2_,default_,p154_] := E2*p154
function41[Lut4_,default_,facP4_,p166_,p167_] := (facP4*(p166 + default*Lut4*p167))/default
end function definitions

begin events
event_1 = (t >= p269, {Agod[t] -> p272/default}, 0)
event_2 = (t >= p469, {Antd[t] -> p472 + Antd[t]}, 0)
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units