begin name
goldbeter6
end name

begin reactions
reaction0 {1.0}$EXTERNAL = {1.0}P
reaction1 {1.0}P = {1.0}$EXTERNAL
reaction2 {1.0}$EXTERNAL = {1.0}Q
reaction3 {1.0}Q = {1.0}$EXTERNAL
reaction4 {1.0}$EXTERNAL = {1.0}R
reaction5 {1.0}R = {1.0}$EXTERNAL
end reactions

begin rate equations
reaction0 = body*function0[reaction0a, Q[t]]
reaction1 = function1[reaction1V, P[t], reaction1Km]
reaction2 = function2[reaction2V1, Q[t], reaction2K1]
reaction3 = function3[reaction3V2, R[t], Q[t], reaction3K2]
reaction4 = function4[P[t], reaction4V3, R[t], reaction4k3]
reaction5 = function1[reaction5V, R[t], reaction5Km]
end rate equations

begin parameters
reaction5Km = 0.01
reaction4k3 = 0.01
reaction0a = 0.1
reaction1V = 0.1
reaction1Km = 0.2
reaction2V1 = 1.0
reaction2K1 = 0.01
reaction3V2 = 1.5
reaction3K2 = 0.01
reaction4V3 = 6.0
reaction5V = 2.5
body = 1.0
end parameters

begin initial conditions
P[0] = Pi
Q[0] = Qi
R[0] = Ri
end initial conditions

begin initial values
Pi = 0.0
Qi = 0.0
Ri = 0.0
end initial values

begin assignment rules
end assignment rules

begin function definitions
function4[P_,V3_,R_,k3_] := (P*(1 - R)*V3)/(1 + k3 - R)
function0[a_,Q_] := a*Q
function1[V_,substrate_,Km_] := (substrate*V)/(Km + substrate)
function2[V1_,Q_,K1_] := ((1 - Q)*V1)/(1 + K1 - Q)
function3[V2_,R_,Q_,K2_] := (Q*R*V2)/(K2 + Q)
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units