simulates a CRNN ode model with extra work — gsl_odeiv2

This function calls a C function which solves an initial value problem, derived from a CRNN.

Usage

gsl_odeiv2_CRNN(
  name,
  experiments,
  l,
  nu,
  m,
  abs.tol = 1e-06,
  rel.tol = 1e-05,
  initial.step.size = 0.001,
  method = 0
)

name: the name of the ODE model to simulate (a shared library of the same name will be dynamically loaded and needs to be created first)
experiments: a list of N simulation experiments (time, parameters, initial value, events).
l: a matrix of parameters with M columns, in log-space.
nu: a stoichiometry matrix (N×R) where N is the number of state variables and R the number of reactions, all reactions are assumed to be reversible.
m: modifiers -- similar to stoichiometry, but indicates whether the species takes part in the reaction without being consumed.
abs.tol: absolute tolerance, real scalar.
rel.tol: relative tolerance, real scalar.
initial.step.size: initial value for the step size; the step size will adapt to a value that observes the tolerances, real scalar.

a list of the solution trajectories y(t;p) for all experiments (named like the experiments), as well as the output functions.