4.6. Custom equations

4.6.1. Complex flux rate equations

The simplest way to script custom rate equations is to use ExpressionFluxProcess. Here is an example taken from the Drosophila sample model. In this expression, Size * N_A of the supersystem of the Process is used to keep the unit of the expression [ num / second ].

Process ExpressionFluxProcess( R_toy1 )
{
    vs      0.76;
    KI      1;
    Expression "(vs*KI) / (KI + C0.MolarConc ^ 3) 
                                       * self.getSuperSystem().SizeN_A";

    VariableReferenceList   [ P0 :.:M 1 ] [ C0 :.:Pn 0 ];
}

FIXME: some more examples

4.6.2. Algebraic equations

Use of ExpressionAlgebraicProcess is the easiest method to describe algebraic equations.

Be careful about the coefficients of the VariableReferences. (Usually just set unities.)

FIXME: some more examples here