########## NAME
M1
########## METADATA
time_unit = m
########## MACROS
########## STATES
d/dt(R) = r3+r2-r1
d/dt(Rp) = r1-r2-r3

R(0) = 1.0
Rp(0) = 0.0

########## PARAMETERS
k1 = 1.0
k2 = 0.0001
kfeed = 1000000.0

########## VARIABLES
r1 = R*A*k1
r2 = Rp*k2
r3 = Rp*kfeed

########## FUNCTIONS
########## EVENTS
########## OUTPUTS
########## INPUTS
A = A_in @ 1
########## FEATURES
y_sim = Rp