## Not run:
#
# ## Generate a compiled ODE model from an equation vector
# ## The model will not return sensitivities for "switch"
# ## Files will be generated in your working directory!
#
# f <- eqnvec(A = "-k*A + switch*F")
# model <- odemodel(f, forcings = "F", fixed = "switch")
# print(model)
#
# ## Generate the same model from an equation list
# f <- addReaction(NULL, from = "", to = "A", rate = "switch*F", description = "production")
# f <- addReaction(f , from = "A", to = "", rate = "k*A", description = "degradation")
# print(f)
#
# model <- odemodel(f, forcings = "F", fixed = "switch")
# print(model)
#
# ## End(Not run)
Run the code above in your browser using DataLab