Expand third order sensitivity

Builds the symengine code for the third order forward sensitivity d/dt(d^3 state / d s1 d s2 d s3) by total-differentiating the stored second order sensitivity RHS (for the pair s2,s3) with respect to the newly added variable s1. Mirrors rxExpandSens2_; assumes the first and second order sensitivities (and the Jacobian) are already loaded in the symengine model environment.

Usage

rxExpandSens3_(state, s1, s2, s3)

Arguments

state

is the state to expand

s1

is the newly added (outer) sensitivity variable

s2

is the second sensitivity variable

s3

is the third sensitivity variable

Value

data frame of symengine third order sensitivity code