Model Mathematics

Component: environment

Component: M

ddtimeM=- k1 M + k2 Mp + gx AM

Component: Mp

ddtimeMp= k1 M +- k2 - fp Mp + gp AMp

Component: AMp

ddtimeAMp= fp Mp +- k5 - gp AMp + k6 AM

Component: AM

ddtimeAM= k5 AMp +- k6 - gx AM

Component: model_parameters

k1 = 0.35 if time < 0.0 time < 5.0 0.060 otherwise k6 = k1 fp = 0.0 if x < 0.0 fp1 x h if x 0.0 x h 0.0 if h < x gp = gp2 if x < 0.0 gp1 x h if x 0.0 x h gp1 + gp3 x h if h < x gx = g2 if x < 0.0 g1 x h if x 0.0 x h g1 + g3 x h if h < x gp2 = 4.0 fp1 + gp1 g2 = 20.0 g1 g3 = 3.0 g1 gp3 = 3.0 gp1