@ t0=0,total=1,dt=0.01,METH=stiff

# XPPAUT content, generated by COPASI from a Xholon-generated SBML file

# Initial values:
#metabolite 'B': reactions
init B=180
#metabolite 'C': reactions
init C=200
#metabolite 'E': reactions
init E=80
 
# Fixed Model Entities:
#compartment 'ReactionNetworkSystem':fixed
param Reactio=1
#compartment 'pot':fixed
param pot=1
#reaction 'A_BB':  kinetic parameter 'k1'
param k1=0.0039
#reaction 'BB_A':  kinetic parameter 'k1'
param k11=0.0038
#reaction 'AC_D':  kinetic parameter 'k1'
param k12=0.0037
#reaction 'D_AC':  kinetic parameter 'k1'
param k13=0.0036
#reaction 'D_BE':  kinetic parameter 'k1'
param k14=0.0035
#reaction 'BE_AC':  kinetic parameter 'k1'
param k15=0.0034
 
# Assignment Model Entities:
#metabolite 'D': reactions
D=305-C-1*E
#metabolite 'A': reactions
A=-10.00000000000001-0.5*B+C+0.5*E
B_c=B/pot
C_c=C/pot
E_c=E/pot
D_c=D/pot
A_c=A/pot
 
 
 
#Kinetics:  
MassAct=k1*A_c
MassAc=k11*B_c*B_c
MassAc1=k12*A_c*C_c
MassAc2=k13*D_c
MassAc3=k14*D_c
MassAc4=k15*B_c*E_c
 
# Equations:
dB/dt=2*MassAct*pot-2*MassAc*pot+MassAc3*pot-MassAc4*pot
dC/dt=-MassAc1*pot+MassAc2*pot+MassAc4*pot
dE/dt=MassAc3*pot-MassAc4*pot
 done