Hi!<br><br><br>I am trying to implement a DAE index reduction algorithm using CppAD + CppADCodeGen which requires time differentiation of equations.<br><br>Consider the following equation:<br><br>dx(t)/dt = x(t) * y(t)<br>
<br>Such equation might have to be differentiated (an unknown number of times) relative to time (t). <br>For instance, the first time derivative  of the equation would be:<br><br>dx(t)/dt^2 = dx(t)/dt * y(t) + x(t) * dy(t)/dt<br>
<br>Any ideas regarding what would be the best strategy to implement this using CppAD?<br>I am considering using CPPAD_USER_ATOMIC to make the relation between x - dxdt - dxdt^2 - ...<br><br>Is there a better alternative?<br>
<br><br>Best regards,<br><br>Joćo Leal<br>