(* Generated by JWS Online *) (* This is an experimental feature of JWS Online. Please report any mistakes.*) (* Note that the following notable SBML entities or features are not supported in notebook outputyet: *) (* Events *) (* Constraints *) (* Units and UnitDefinitions *) (* AlgebraicRules *) (* conversionFactors *) variables = { AF1[t], AF2[t], AF3[t], AF4[t], Ago\[LetterSpace]R\[LetterSpace]a[t], Ago\[LetterSpace]R\[LetterSpace]i[t], Ago\[LetterSpace]c[t], Ago\[LetterSpace]d[t], Ant\[LetterSpace]R[t], Ant\[LetterSpace]c[t], Ant\[LetterSpace]d[t], Ant\[LetterSpace]p[t], E2[t], FSH\[LetterSpace]R[t], FSH\[LetterSpace]bld[t], FSH\[LetterSpace]pit[t], GnRH[t], GnRH\[LetterSpace]R\[LetterSpace]a[t], GnRH\[LetterSpace]R\[LetterSpace]i[t], InhA[t], InhA\[LetterSpace]delay[t], InhB[t], LH\[LetterSpace]Pit[t], LH\[LetterSpace]R[t], LH\[LetterSpace]bld[t], Lut1[t], Lut2[t], Lut3[t], Lut4[t], OvF[t], P4[t], PrF[t], R\[LetterSpace]FSH[t], R\[LetterSpace]FSH\[LetterSpace]des[t], R\[LetterSpace]Foll[t], R\[LetterSpace]GnRH\[LetterSpace]a[t], R\[LetterSpace]GnRH\[LetterSpace]i[t], R\[LetterSpace]LH[t], R\[LetterSpace]LH\[LetterSpace]des[t], Sc1[t], Sc2[t], s113[t], s114[t], s115[t], s116[t] }; initialValues = { AF1[0] == 2.26, AF2[0] == 19.92, AF3[0] == 0.504, AF4[0] == 1.604*^-05, Ago\[LetterSpace]R\[LetterSpace]a[0] == 0.0, Ago\[LetterSpace]R\[LetterSpace]i[0] == 0.0, Ago\[LetterSpace]c[0] == 0.0, Ago\[LetterSpace]d[0] == 0.0, Ant\[LetterSpace]R[0] == 0.0, Ant\[LetterSpace]c[0] == 0.0, Ant\[LetterSpace]d[0] == 0.0, Ant\[LetterSpace]p[0] == 0.0, E2[0] == 45.0, FSH\[LetterSpace]R[0] == 0.662, FSH\[LetterSpace]bld[0] == 4.1, FSH\[LetterSpace]pit[0] == 48627.2, GnRH[0] == 0.00326, GnRH\[LetterSpace]R\[LetterSpace]a[0] == 1.447*^-05, GnRH\[LetterSpace]R\[LetterSpace]i[0] == 1.294*^-05, InhA[0] == 1.07, InhA\[LetterSpace]delay[0] == 86.84, InhB[0] == 52.18, LH\[LetterSpace]Pit[0] == 261119.0, LH\[LetterSpace]R[0] == 0.263, LH\[LetterSpace]bld[0] == 2.667, Lut1[0] == 2.762*^-05, Lut2[0] == 0.0003772, Lut3[0] == 0.00334, Lut4[0] == 0.01385, OvF[0] == 7.652*^-19, P4[0] == 2.41, PrF[0] == 0.244, R\[LetterSpace]FSH[0] == 6.341, R\[LetterSpace]FSH\[LetterSpace]des[0] == 1.497, R\[LetterSpace]Foll[0] == 0.119, R\[LetterSpace]GnRH\[LetterSpace]a[0] == 0.00928, R\[LetterSpace]GnRH\[LetterSpace]i[0] == 0.0009409, R\[LetterSpace]LH[0] == 8.41, R\[LetterSpace]LH\[LetterSpace]des[0] == 0.699, Sc1[0] == 1.098*^-08, Sc2[0] == 2.171*^-06, s113[0] == 0.0, s114[0] == 0.0, s115[0] == 0.0, s116[0] == 0.0 }; rates = { re100, re101, re102, re103, re104, re105, re106, re107, re108, re109, re11, re15, re2, re24, re25, re26, re28, re29, re3, re35, re36, re37, re38, re39, re4, re40, re42, re43, re44, re45, re46, re49, re5, re50, re51, re52, re53, re54, re56, re57, re58, re59, re6, re60, re61, re62, re64, re65, re66, re67, re69, re71, re72, re73, re74, re75, re76, re78, re8, re82, re83, re84, re85, re87, re90, re93, re94, re95, re97, re98, re99 }; rateEquations = { re100 -> default*function\[LetterSpace]62[Ago\[LetterSpace]d[t], default, p274], re101 -> default*function\[LetterSpace]63[Ago\[LetterSpace]d[t], default, p273, p274], re102 -> default*function\[LetterSpace]64[Ant\[LetterSpace]d[t], default, p474], re103 -> default*function\[LetterSpace]65[Ant\[LetterSpace]d[t], default, p473, p474], re104 -> default*function\[LetterSpace]66[Ant\[LetterSpace]c[t], default, p475], re105 -> default*function\[LetterSpace]67[Ant\[LetterSpace]R[t], default, p513], re106 -> default*function\[LetterSpace]68[Ant\[LetterSpace]c[t], R\[LetterSpace]GnRH\[LetterSpace]a[t], default, p512], re107 -> default*function\[LetterSpace]69[Ant\[LetterSpace]R[t], default, p514], re108 -> default*function\[LetterSpace]70[Ant\[LetterSpace]c[t], default, p476], re109 -> default*function\[LetterSpace]71[Ant\[LetterSpace]p[t], default, p477], re11 -> default*function\[LetterSpace]7[R\[LetterSpace]GnRH\[LetterSpace]i[t], default, p308], re15 -> default*function\[LetterSpace]8[GnRH\[LetterSpace]R\[LetterSpace]a[t], default, p303], re2 -> default*function\[LetterSpace]1[GnRH\[LetterSpace]R\[LetterSpace]a[t], default, p309], re24 -> default*function\[LetterSpace]9[GnRH[t], default, p300], re25 -> default*function\[LetterSpace]10[E2[t], P4[t], default, facE2, facP4, p1, p2, p3, p4, p6, p7], re26 -> default*function\[LetterSpace]11[InhA\[LetterSpace]delay[t], InhB[t], default, freq, p11, p13, p21, p22, p23, p24, p25], re28 -> default*function\[LetterSpace]12[Ago\[LetterSpace]R\[LetterSpace]a[t], GnRH\[LetterSpace]R\[LetterSpace]a[t], LH\[LetterSpace]Pit[t], default, p16, p5, p8, p9], re29 -> default*function\[LetterSpace]13[Ago\[LetterSpace]R\[LetterSpace]a[t], FSH\[LetterSpace]pit[t], GnRH\[LetterSpace]R\[LetterSpace]a[t], default, p17, p18, p20, p28], re3 -> default*function\[LetterSpace]2[GnRH\[LetterSpace]R\[LetterSpace]i[t], default, p310], re35 -> default*function\[LetterSpace]14[LH\[LetterSpace]bld[t], R\[LetterSpace]LH[t], default, facLH, p230], re36 -> default*function\[LetterSpace]15[LH\[LetterSpace]R[t], default, p234], re37 -> default*function\[LetterSpace]16[R\[LetterSpace]LH\[LetterSpace]des[t], default, p232], re38 -> default*function\[LetterSpace]17[LH\[LetterSpace]bld[t], default, p231], re39 -> default*function\[LetterSpace]18[FSH\[LetterSpace]R[t], default, p244], re4 -> default*function\[LetterSpace]3[R\[LetterSpace]GnRH\[LetterSpace]i[t], default, p307], re40 -> default*function\[LetterSpace]19[R\[LetterSpace]FSH\[LetterSpace]des[t], default, p242], re42 -> default*function\[LetterSpace]20[AF1[t], FSH\[LetterSpace]R[t], default, p50], re43 -> default*function\[LetterSpace]21[FSH\[LetterSpace]R[t], default, p47, p48, p49], re44 -> default*function\[LetterSpace]22[AF2[t], LH\[LetterSpace]R[t], R\[LetterSpace]Foll[t], default, p46, p51, p52], re45 -> default*function\[LetterSpace]23[AF3[t], LH\[LetterSpace]R[t], R\[LetterSpace]Foll[t], default, p32, p43, p52], re46 -> default*function\[LetterSpace]24[AF4[t], LH\[LetterSpace]R[t], R\[LetterSpace]Foll[t], default, p34, p52], re49 -> default*function\[LetterSpace]25[Sc1[t], default, p37], re5 -> default*function\[LetterSpace]4[R\[LetterSpace]GnRH\[LetterSpace]a[t], default, p306], re50 -> default*function\[LetterSpace]26[Sc2[t], default, p38], re51 -> default*function\[LetterSpace]27[Ago\[LetterSpace]R\[LetterSpace]a[t], GnRH\[LetterSpace]R\[LetterSpace]a[t], Lut1[t], default, p39, p80, p83, p84], re52 -> default*function\[LetterSpace]28[Ago\[LetterSpace]R\[LetterSpace]a[t], GnRH\[LetterSpace]R\[LetterSpace]a[t], Lut2[t], default, p40, p80, p83, p84], re53 -> default*function\[LetterSpace]29[Ago\[LetterSpace]R\[LetterSpace]a[t], GnRH\[LetterSpace]R\[LetterSpace]a[t], Lut3[t], default, p41, p80, p83, p84], re54 -> default*function\[LetterSpace]30[FSH\[LetterSpace]bld[t], default, p90, p91, p94], re56 -> default*function\[LetterSpace]31[P4[t], R\[LetterSpace]Foll[t], default, facP4, p92, p93, p95], re57 -> default*function\[LetterSpace]32[AF3[t], FSH\[LetterSpace]R[t], default, p31, p55], re58 -> default*function\[LetterSpace]33[AF4[t], LH\[LetterSpace]R[t], default, p33, p44, p52, p55], re59 -> default*function\[LetterSpace]34[LH\[LetterSpace]R[t], PrF[t], R\[LetterSpace]Foll[t], default, p35, p45, p52], re6 -> default*function\[LetterSpace]5[GnRH\[LetterSpace]R\[LetterSpace]i[t], default, p305], re60 -> default*function\[LetterSpace]35[LH\[LetterSpace]R[t], PrF[t], R\[LetterSpace]Foll[t], default, p27, p45, p52, p53, p54], re61 -> default*function\[LetterSpace]36[OvF[t], default, p36], re62 -> default*function\[LetterSpace]37[OvF[t], default, p26, p56, p57], re64 -> default*function\[LetterSpace]38[Ago\[LetterSpace]R\[LetterSpace]a[t], GnRH\[LetterSpace]R\[LetterSpace]a[t], Lut4[t], default, p42, p80, p83, p84], re65 -> default*function\[LetterSpace]39[AF2[t], AF3[t], AF4[t], LH\[LetterSpace]bld[t], Lut1[t], Lut4[t], PrF[t], default, facE2, p152, p158, p159, p160, p161, p164, p165], re66 -> default*function\[LetterSpace]40[E2[t], default, p154], re67 -> default*function\[LetterSpace]41[Lut4[t], default, facP4, p166, p167], re69 -> default*function\[LetterSpace]42[Lut1[t], Lut2[t], Lut3[t], Lut4[t], PrF[t], Sc1[t], default, p168, p169, p170, p171, p172, p177, p178], re71 -> default*function\[LetterSpace]43[P4[t], default, p155], re72 -> default*function\[LetterSpace]44[AF2[t], Sc2[t], default, p173, p174, p175], re73 -> default*function\[LetterSpace]45[InhB[t], default, p157], re74 -> default*function\[LetterSpace]46[InhA\[LetterSpace]delay[t], default, p30], re75 -> default*function\[LetterSpace]47[default, p311], re76 -> default*function\[LetterSpace]48[GnRH\[LetterSpace]R\[LetterSpace]i[t], default, p304], re78 -> default*function\[LetterSpace]49[default, freq, mass], re8 -> default*function\[LetterSpace]6[GnRH[t], R\[LetterSpace]GnRH\[LetterSpace]a[t], default, p302], re82 -> default*function\[LetterSpace]50[FSH\[LetterSpace]bld[t], default, p241], re83 -> default*function\[LetterSpace]51[Ago\[LetterSpace]R\[LetterSpace]a[t], GnRH\[LetterSpace]R\[LetterSpace]a[t], LH\[LetterSpace]Pit[t], default, facLH, p12, p16, p5, p8, p9], re84 -> default*function\[LetterSpace]52[Ago\[LetterSpace]R\[LetterSpace]a[t], FSH\[LetterSpace]pit[t], GnRH\[LetterSpace]R\[LetterSpace]a[t], default, facFSH, p12, p17, p18, p20, p28], re85 -> default*function\[LetterSpace]53[FSH\[LetterSpace]bld[t], R\[LetterSpace]FSH[t], default, facFSH, p240], re87 -> default*function\[LetterSpace]54[InhA[t], default, p156], re90 -> default*function\[LetterSpace]55[Ago\[LetterSpace]c[t], default, p275], re93 -> default*function\[LetterSpace]56[Ago\[LetterSpace]R\[LetterSpace]a[t], default, p319], re94 -> default*function\[LetterSpace]57[Ago\[LetterSpace]R\[LetterSpace]i[t], default, p320], re95 -> default*function\[LetterSpace]58[Ago\[LetterSpace]R\[LetterSpace]i[t], default, p314], re97 -> default*function\[LetterSpace]59[Ago\[LetterSpace]c[t], R\[LetterSpace]GnRH\[LetterSpace]a[t], default, p312], re98 -> default*function\[LetterSpace]60[Ago\[LetterSpace]R\[LetterSpace]a[t], default, p313], re99 -> default*function\[LetterSpace]61[Ago\[LetterSpace]R\[LetterSpace]i[t], default, p315] }; 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.94934669769107*^-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, csa1\[LetterSpace]degraded -> 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, sa1\[LetterSpace]degraded -> 1.0, sa28\[LetterSpace]degraded -> 1.0, sa31\[LetterSpace]degraded -> 1.0, sa35\[LetterSpace]degraded -> 1.0, sa3\[LetterSpace]degraded -> 1.0, sa52\[LetterSpace]degraded -> 1.0, sa53\[LetterSpace]degraded -> 1.0, sa61\[LetterSpace]degraded -> 1.0, sa75\[LetterSpace]degraded -> 1.0, sa78\[LetterSpace]degraded -> 1.0, sa86\[LetterSpace]degraded -> 1.0, default -> 1.0 }; assignments = { function\[LetterSpace]44[AF2_,Sc2_,default_,p173_,p174_,p175_] -> (p173 + AF2*default*p174 + default*p175*Sc2)/default, function\[LetterSpace]45[InhB_,default_,p157_] -> InhB*p157, function\[LetterSpace]46[InhA\[LetterSpace]delay_,default_,p30_] -> InhA\[LetterSpace]delay*p30, function\[LetterSpace]47[default_,p311_] -> p311/default, function\[LetterSpace]49[default_,freq_,mass_] -> (freq*mass)/default, function\[LetterSpace]43[P4_,default_,p155_] -> p155*P4, function\[LetterSpace]53[FSH\[LetterSpace]bld_,R\[LetterSpace]FSH_,default_,facFSH_,p240_] -> (default*FSH\[LetterSpace]bld*p240*R\[LetterSpace]FSH)/facFSH, function\[LetterSpace]54[InhA_,default_,p156_] -> InhA*p156, function\[LetterSpace]55[Ago\[LetterSpace]c_,default_,p275_] -> Ago\[LetterSpace]c*p275, function\[LetterSpace]57[Ago\[LetterSpace]R\[LetterSpace]i_,default_,p320_] -> Ago\[LetterSpace]R\[LetterSpace]i*p320, function\[LetterSpace]58[Ago\[LetterSpace]R\[LetterSpace]i_,default_,p314_] -> Ago\[LetterSpace]R\[LetterSpace]i*p314, function\[LetterSpace]59[Ago\[LetterSpace]c_,R\[LetterSpace]GnRH\[LetterSpace]a_,default_,p312_] -> Ago\[LetterSpace]c*default*p312*R\[LetterSpace]GnRH\[LetterSpace]a, function\[LetterSpace]62[Ago\[LetterSpace]d_,default_,p274_] -> Ago\[LetterSpace]d*p274, function\[LetterSpace]63[Ago\[LetterSpace]d_,default_,p273_,p274_] -> (Ago\[LetterSpace]d*p274)/p273, function\[LetterSpace]64[Ant\[LetterSpace]d_,default_,p474_] -> (Ant\[LetterSpace]d*p474)/default, function\[LetterSpace]65[Ant\[LetterSpace]d_,default_,p473_,p474_] -> (Ant\[LetterSpace]d*p474)/(default*p473), function\[LetterSpace]66[Ant\[LetterSpace]c_,default_,p475_] -> (Ant\[LetterSpace]c*p475)/default, function\[LetterSpace]67[Ant\[LetterSpace]R_,default_,p513_] -> (Ant\[LetterSpace]R*p513)/default, function\[LetterSpace]69[Ant\[LetterSpace]R_,default_,p514_] -> (Ant\[LetterSpace]R*p514)/default, function\[LetterSpace]71[Ant\[LetterSpace]p_,default_,p477_] -> (Ant\[LetterSpace]p*p477)/default, function\[LetterSpace]70[Ant\[LetterSpace]c_,default_,p476_] -> (Ant\[LetterSpace]c*p476)/default, function\[LetterSpace]68[Ant\[LetterSpace]c_,R\[LetterSpace]GnRH\[LetterSpace]a_,default_,p512_] -> Ant\[LetterSpace]c*p512*R\[LetterSpace]GnRH\[LetterSpace]a, function\[LetterSpace]61[Ago\[LetterSpace]R\[LetterSpace]i_,default_,p315_] -> Ago\[LetterSpace]R\[LetterSpace]i*p315, function\[LetterSpace]60[Ago\[LetterSpace]R\[LetterSpace]a_,default_,p313_] -> Ago\[LetterSpace]R\[LetterSpace]a*p313, function\[LetterSpace]56[Ago\[LetterSpace]R\[LetterSpace]a_,default_,p319_] -> Ago\[LetterSpace]R\[LetterSpace]a*p319, function\[LetterSpace]52[Ago\[LetterSpace]R\[LetterSpace]a_,FSH\[LetterSpace]pit_,GnRH\[LetterSpace]R\[LetterSpace]a_,default_,facFSH_,p12_,p17_,p18_,p20_,p28_] -> (facFSH*FSH\[LetterSpace]pit*(p17 + (((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p18)^p20*p28)/(1 + ((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p18)^p20)))/p12, function\[LetterSpace]51[Ago\[LetterSpace]R\[LetterSpace]a_,GnRH\[LetterSpace]R\[LetterSpace]a_,LH\[LetterSpace]Pit_,default_,facLH_,p12_,p16_,p5_,p8_,p9_] -> (facLH*LH\[LetterSpace]Pit*(p16 + (p5*((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p8)^p9)/(1 + ((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p8)^p9)))/p12, function\[LetterSpace]50[FSH\[LetterSpace]bld_,default_,p241_] -> FSH\[LetterSpace]bld*p241, function\[LetterSpace]48[GnRH\[LetterSpace]R\[LetterSpace]i_,default_,p304_] -> GnRH\[LetterSpace]R\[LetterSpace]i*p304, function\[LetterSpace]42[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, function\[LetterSpace]38[Ago\[LetterSpace]R\[LetterSpace]a_,GnRH\[LetterSpace]R\[LetterSpace]a_,Lut4_,default_,p42_,p80_,p83_,p84_] -> Lut4*p42*(1 + (p80*((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p83)^p84)/(1 + ((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p83)^p84)), function\[LetterSpace]31[P4_,R\[LetterSpace]Foll_,default_,facP4_,p92_,p93_,p95_] -> (((default*P4)/(facP4*p92))^p93*p95*R\[LetterSpace]Foll)/(1 + ((default*P4)/(facP4*p92))^p93), function\[LetterSpace]29[Ago\[LetterSpace]R\[LetterSpace]a_,GnRH\[LetterSpace]R\[LetterSpace]a_,Lut3_,default_,p41_,p80_,p83_,p84_] -> Lut3*p41*(1 + (p80*((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p83)^p84)/(1 + ((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p83)^p84)), function\[LetterSpace]27[Ago\[LetterSpace]R\[LetterSpace]a_,GnRH\[LetterSpace]R\[LetterSpace]a_,Lut1_,default_,p39_,p80_,p83_,p84_] -> Lut1*p39*(1 + (p80*((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p83)^p84)/(1 + ((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p83)^p84)), function\[LetterSpace]21[FSH\[LetterSpace]R_,default_,p47_,p48_,p49_] -> (((default*FSH\[LetterSpace]R)/p48)^p47*p49)/(default*(1 + ((default*FSH\[LetterSpace]R)/p48)^p47)), function\[LetterSpace]13[Ago\[LetterSpace]R\[LetterSpace]a_,FSH\[LetterSpace]pit_,GnRH\[LetterSpace]R\[LetterSpace]a_,default_,p17_,p18_,p20_,p28_] -> FSH\[LetterSpace]pit*(p17 + (((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p18)^p20*p28)/(1 + ((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p18)^p20)), function\[LetterSpace]12[Ago\[LetterSpace]R\[LetterSpace]a_,GnRH\[LetterSpace]R\[LetterSpace]a_,LH\[LetterSpace]Pit_,default_,p16_,p5_,p8_,p9_] -> LH\[LetterSpace]Pit*(p16 + (p5*((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p8)^p9)/(1 + ((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p8)^p9)), function\[LetterSpace]1[GnRH\[LetterSpace]R\[LetterSpace]a_,default_,p309_] -> GnRH\[LetterSpace]R\[LetterSpace]a*p309, function\[LetterSpace]2[GnRH\[LetterSpace]R\[LetterSpace]i_,default_,p310_] -> GnRH\[LetterSpace]R\[LetterSpace]i*p310, function\[LetterSpace]3[R\[LetterSpace]GnRH\[LetterSpace]i_,default_,p307_] -> p307*R\[LetterSpace]GnRH\[LetterSpace]i, function\[LetterSpace]4[R\[LetterSpace]GnRH\[LetterSpace]a_,default_,p306_] -> p306*R\[LetterSpace]GnRH\[LetterSpace]a, function\[LetterSpace]5[GnRH\[LetterSpace]R\[LetterSpace]i_,default_,p305_] -> GnRH\[LetterSpace]R\[LetterSpace]i*p305, function\[LetterSpace]6[GnRH_,R\[LetterSpace]GnRH\[LetterSpace]a_,default_,p302_] -> default*GnRH*p302*R\[LetterSpace]GnRH\[LetterSpace]a, function\[LetterSpace]7[R\[LetterSpace]GnRH\[LetterSpace]i_,default_,p308_] -> p308*R\[LetterSpace]GnRH\[LetterSpace]i, function\[LetterSpace]8[GnRH\[LetterSpace]R\[LetterSpace]a_,default_,p303_] -> GnRH\[LetterSpace]R\[LetterSpace]a*p303, function\[LetterSpace]9[GnRH_,default_,p300_] -> GnRH*p300, function\[LetterSpace]10[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)), function\[LetterSpace]11[InhA\[LetterSpace]delay_,InhB_,default_,freq_,p11_,p13_,p21_,p22_,p23_,p24_,p25_] -> p21/(default*(1 + (freq/p11)^p13)*(1 + ((default*InhA\[LetterSpace]delay)/p22)^p24 + ((default*InhB)/p23)^p25)), function\[LetterSpace]14[LH\[LetterSpace]bld_,R\[LetterSpace]LH_,default_,facLH_,p230_] -> (default*LH\[LetterSpace]bld*p230*R\[LetterSpace]LH)/facLH, function\[LetterSpace]15[LH\[LetterSpace]R_,default_,p234_] -> LH\[LetterSpace]R*p234, function\[LetterSpace]16[R\[LetterSpace]LH\[LetterSpace]des_,default_,p232_] -> p232*R\[LetterSpace]LH\[LetterSpace]des, function\[LetterSpace]17[LH\[LetterSpace]bld_,default_,p231_] -> LH\[LetterSpace]bld*p231, function\[LetterSpace]18[FSH\[LetterSpace]R_,default_,p244_] -> FSH\[LetterSpace]R*p244, function\[LetterSpace]19[R\[LetterSpace]FSH\[LetterSpace]des_,default_,p242_] -> p242*R\[LetterSpace]FSH\[LetterSpace]des, function\[LetterSpace]20[AF1_,FSH\[LetterSpace]R_,default_,p50_] -> AF1*default*FSH\[LetterSpace]R*p50, function\[LetterSpace]25[Sc1_,default_,p37_] -> p37*Sc1, function\[LetterSpace]26[Sc2_,default_,p38_] -> p38*Sc2, function\[LetterSpace]22[AF2_,LH\[LetterSpace]R_,R\[LetterSpace]Foll_,default_,p46_,p51_,p52_] -> AF2*default*p51*((default*LH\[LetterSpace]R)/p52)^p46*R\[LetterSpace]Foll, function\[LetterSpace]23[AF3_,LH\[LetterSpace]R_,R\[LetterSpace]Foll_,default_,p32_,p43_,p52_] -> AF3*default*p32*((default*LH\[LetterSpace]R)/p52)^p43*R\[LetterSpace]Foll, function\[LetterSpace]24[AF4_,LH\[LetterSpace]R_,R\[LetterSpace]Foll_,default_,p34_,p52_] -> (AF4*default^2*LH\[LetterSpace]R*p34*R\[LetterSpace]Foll)/p52, function\[LetterSpace]28[Ago\[LetterSpace]R\[LetterSpace]a_,GnRH\[LetterSpace]R\[LetterSpace]a_,Lut2_,default_,p40_,p80_,p83_,p84_] -> Lut2*p40*(1 + (p80*((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p83)^p84)/(1 + ((Ago\[LetterSpace]R\[LetterSpace]a*default + default*GnRH\[LetterSpace]R\[LetterSpace]a)/p83)^p84)), function\[LetterSpace]30[FSH\[LetterSpace]bld_,default_,p90_,p91_,p94_] -> (((default*FSH\[LetterSpace]bld)/p90)^p91*p94)/(default*(1 + ((default*FSH\[LetterSpace]bld)/p90)^p91)), function\[LetterSpace]32[AF3_,FSH\[LetterSpace]R_,default_,p31_,p55_] -> AF3*default*FSH\[LetterSpace]R*p31*(1 - (AF3*default)/p55), function\[LetterSpace]33[AF4_,LH\[LetterSpace]R_,default_,p33_,p44_,p52_,p55_] -> AF4*p33*((default*LH\[LetterSpace]R)/p52)^p44*(1 - (AF4*default)/p55), function\[LetterSpace]34[LH\[LetterSpace]R_,PrF_,R\[LetterSpace]Foll_,default_,p35_,p45_,p52_] -> default*p35*((default*LH\[LetterSpace]R)/p52)^p45*PrF*R\[LetterSpace]Foll, function\[LetterSpace]35[LH\[LetterSpace]R_,PrF_,R\[LetterSpace]Foll_,default_,p27_,p45_,p52_,p53_,p54_] -> (p27*((default*LH\[LetterSpace]R)/p52)^p45*((default*PrF)/p53)^p54*R\[LetterSpace]Foll)/(1 + ((default*PrF)/p53)^p54), function\[LetterSpace]36[OvF_,default_,p36_] -> OvF*p36, function\[LetterSpace]37[OvF_,default_,p26_,p56_,p57_] -> (p26*((default*OvF)/p56)^p57)/(default*(1 + ((default*OvF)/p56)^p57)), function\[LetterSpace]39[AF2_,AF3_,AF4_,LH\[LetterSpace]bld_,Lut1_,Lut4_,PrF_,default_,facE2_,p152_,p158_,p159_,p160_,p161_,p164_,p165_] -> (facE2*(AF2*default*p152 + p158 + AF3*default^2*LH\[LetterSpace]bld*p159 + AF4*default*p160 + default*Lut1*p164 + default*Lut4*p165 + default^2*LH\[LetterSpace]bld*p161*PrF))/default, function\[LetterSpace]40[E2_,default_,p154_] -> E2*p154, function\[LetterSpace]41[Lut4_,default_,facP4_,p166_,p167_] -> (facP4*(p166 + default*Lut4*p167))/default, 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) }; events = { }; speciesAnnotations = { }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { AF1'[t] == 1.0*re43 -1.0*re42, AF2'[t] == 1.0*re42 -1.0*re44, AF3'[t] == 1.0*re44 +1.0*re57 -1.0*re45, AF4'[t] == 1.0*re45 +1.0*re58 -1.0*re46, Ago\[LetterSpace]R\[LetterSpace]a'[t] == 1.0*re94 +1.0*re97 -1.0*re93 -1.0*re98, Ago\[LetterSpace]R\[LetterSpace]i'[t] == 1.0*re93 -1.0*re94 -1.0*re95 -1.0*re99, Ago\[LetterSpace]c'[t] == 1.0*re98 +1.0*re101 -1.0*re90 -1.0*re97, Ago\[LetterSpace]d'[t] == -1.0*re100, Ant\[LetterSpace]R'[t] == 1.0*re106 -1.0*re105 -1.0*re107, Ant\[LetterSpace]c'[t] == 1.0*re103 +1.0*re105 +1.0*re109 -1.0*re104 -1.0*re106 -1.0*re108, Ant\[LetterSpace]d'[t] == -1.0*re102, Ant\[LetterSpace]p'[t] == 1.0*re108 -1.0*re109, E2'[t] == 1.0*re65 -1.0*re66, FSH\[LetterSpace]R'[t] == 1.0*re85 -1.0*re39, FSH\[LetterSpace]bld'[t] == 1.0*re84 -1.0*re82 -1.0*re85, FSH\[LetterSpace]pit'[t] == 1.0*re26 -1.0*re29, GnRH'[t] == 1.0*re15 +1.0*re78 -1.0*re8 -1.0*re24, GnRH\[LetterSpace]R\[LetterSpace]a'[t] == 1.0*re3 +1.0*re8 -1.0*re2 -1.0*re15, GnRH\[LetterSpace]R\[LetterSpace]i'[t] == 1.0*re2 -1.0*re3 -1.0*re6 -1.0*re76, InhA'[t] == 1.0*re69 -1.0*re87, InhA\[LetterSpace]delay'[t] == 1.0*re87 -1.0*re74, InhB'[t] == 1.0*re72 -1.0*re73, LH\[LetterSpace]Pit'[t] == 1.0*re25 -1.0*re28, LH\[LetterSpace]R'[t] == 1.0*re35 -1.0*re36, LH\[LetterSpace]bld'[t] == 1.0*re83 -1.0*re35 -1.0*re38, Lut1'[t] == 1.0*re50 -1.0*re51, Lut2'[t] == 1.0*re51 -1.0*re52, Lut3'[t] == 1.0*re52 -1.0*re53, Lut4'[t] == 1.0*re53 -1.0*re64, OvF'[t] == 1.0*re60 -1.0*re61, P4'[t] == 1.0*re67 -1.0*re71, PrF'[t] == 1.0*re46 -1.0*re59, R\[LetterSpace]FSH'[t] == 1.0*re40 -1.0*re85, R\[LetterSpace]FSH\[LetterSpace]des'[t] == 1.0*re39 -1.0*re40, R\[LetterSpace]Foll'[t] == 1.0*re54 -1.0*re56, R\[LetterSpace]GnRH\[LetterSpace]a'[t] == 1.0*re4 +1.0*re15 +1.0*re98 +1.0*re105 -1.0*re5 -1.0*re8 -1.0*re97 -1.0*re106, R\[LetterSpace]GnRH\[LetterSpace]i'[t] == 1.0*re5 +1.0*re6 +1.0*re75 +1.0*re99 -1.0*re4 -1.0*re11, R\[LetterSpace]LH'[t] == 1.0*re37 -1.0*re35, R\[LetterSpace]LH\[LetterSpace]des'[t] == 1.0*re36 -1.0*re37, Sc1'[t] == 1.0*re62 -1.0*re49, Sc2'[t] == 1.0*re49 -1.0*re50, s113'[t] == 1.0*re102 , s114'[t] == -1.0*re103, s115'[t] == 1.0*re104 , s116'[t] == 1.0*re107 }; timeCourse = NDSolve[Join[odes, initialValues]//.rateEquations//.assignments//.parameters, variables, {t, 0, 100}]; (* Steady-state solution initialized with result of time evolution *) findRootEquations = odes /.D[_[t],t]->0; findRootVariables = Partition[Flatten[{#, #/.timeCourse/.t->100} &/@variables],2]; steadyStateVariables = FindRoot[findRootEquations//.rateEquations//.assignments//.parameters, findRootVariables, MaxIterations->100] fluxes = #//.assignments//.parameters/.steadyStateVariables&/@rateEquations (* Plot the time evolution of the variables *) plotTable=Table[Plot[variables[[i]]/.parameters/.timeCourse,{t,0,100},PlotLegends->variables[[i]],PlotRange->Full],{i,Length[variables]}]