#
# This file is automatically generated with 
# the System Biology Format Converter (http://sbfc.sourceforge.net/)
# from an SBML file.
#

#
# Model name = Petri net - feinberg1
#
#

# some function definitions that are allowed in SBML but not valid in xpp
ceil(x)=flr(1+x)

@delay=50


# Compartment: id = reactionNetworkSystem_0c, name = ReactionNetworkSystem, constant
par reactionN=1.0

# Compartment: id = pot_petriNet_1c, name = pot, constant
par pot_petri=1.0

# Reaction: id = a_BB_13r, name = A_BB
# Local Parameter:   id =  K, name = K
par K=0.0039


a_BB_13r=pot_petri*K*a_7s

# Reaction: id = bB_A_18r, name = BB_A
# Local Parameter:   id =  K, name = K
par K_1=0.0038


bB_A_18r=pot_petri*K_1*b_8s^(2)

# Reaction: id = aC_D_23r, name = AC_D
# Local Parameter:   id =  K, name = K
par K_1_2=0.0037


aC_D_23r=pot_petri*K_1_2*a_7s*c_9s

# Reaction: id = d_AC_29r, name = D_AC
# Local Parameter:   id =  K, name = K
par K_1_2_3=0.0036


d_AC_29r=pot_petri*K_1_2_3*d_10s

# Reaction: id = d_BE_35r, name = D_BE
# Local Parameter:   id =  K, name = K
par K_1_2_3_4=0.0035


d_BE_35r=pot_petri*K_1_2_3_4*d_10s

# Reaction: id = bE_AC_41r, name = BE_AC
# Local Parameter:   id =  K, name = K
par K_1_2_3_5=0.0034


bE_AC_41r=pot_petri*K_1_2_3_5*b_8s*e_11s

# Species:   id = a_7s, name = A, affected by kineticLaw

init a_7s=140.0
da_7s/dt=(1/(pot_petri))*((-1.0 * a_BB_13r) + ( 1.0 * bB_A_18r) + (-1.0 * aC_D_23r) + ( 1.0 * d_AC_29r) + ( 1.0 * bE_AC_41r))

# Species:   id = b_8s, name = B, affected by kineticLaw

init b_8s=180.0
db_8s/dt=(1/(pot_petri))*(( 2.0 * a_BB_13r) + (-2.0 * bB_A_18r) + ( 1.0 * d_BE_35r) + (-1.0 * bE_AC_41r))

# Species:   id = c_9s, name = C, affected by kineticLaw

init c_9s=200.0
dc_9s/dt=(1/(pot_petri))*((-1.0 * aC_D_23r) + ( 1.0 * d_AC_29r) + ( 1.0 * bE_AC_41r))

# Species:   id = d_10s, name = D, affected by kineticLaw

init d_10s=25.0
dd_10s/dt=(1/(pot_petri))*(( 1.0 * aC_D_23r) + (-1.0 * d_AC_29r) + (-1.0 * d_BE_35r))

# Species:   id = e_11s, name = E, affected by kineticLaw

init e_11s=80.0
de_11s/dt=(1/(pot_petri))*(( 1.0 * d_BE_35r) + (-1.0 * bE_AC_41r))

@ meth=cvode, tol=1e-6, atol=1e-8
# @ maxstor=1e6
@ bound=40000, total=200
done