begin name
tyson13
end name

begin reactions
r1 {1.0}$EXTERNAL = {1.0}R
r2 {1.0}$EXTERNAL = {1.0}R
r3 {1.0}R = {1.0}$EXTERNAL
r4 {1.0}$EXTERNAL = {1.0}$EXTERNAL
r5 {1.0}$EXTERNAL = {1.0}$EXTERNAL
end reactions

begin rate equations
r1 = env*Ep*k0
r2 = env*k1*S
r3 = env*k2*R[t]
r4 = (E*env*k3*R[t])/(E + J3)
r5 = (env*Ep*k4)/(Ep + J4)
end rate equations

begin parameters
Et = 1.0
J3 = 0.05
J4 = 0.05
k0 = 0.4
k1 = 0.01
k2 = 1.0
k3 = 1.0
k4 = 0.2
S = 0.0
env = 1.0
end parameters

begin initial conditions
R[0] = Ri
end initial conditions

begin initial values
Ri = 0.0
end initial values

begin assignment rules
E := -Ep + Et
Ep := goldbeterkoshland[k3*R[t], k4, J3, J4]
end assignment rules

begin function definitions
goldbeterkoshland[v1_,v2_,J1_,J2_] := (2*J2*v1)/(-v1 + J2*v1 + v2 + J1*v2 + Sqrt[-4*J2*v1*(-v1 + v2) + (-v1 + J2*v1 + v2 + J1*v2)^2])
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units