Hi everyone,
I am using ode15s to solve a set of stiff odes and the final model is developed in the form of
where M(t,y) is the mass matrix, and y is the generalised coordinate vector. According to the syntax of ode15s, it is crucial to provide its Jacobian matrix to raise the accuracy.
I do not know too much numerical calculation of stiff odes, and here I want to know, is the Jacobian matrix JJ calculated in this manner:

If not, how to arrive at a correct one?
Many thanks!Hi everyone,
I am using ode15s to solve a set of stiff odes and the final model is developed in the form of
where M(t,y) is the mass matrix, and y is the generalised coordinate vector. According to the syntax of ode15s, it is crucial to provide its Jacobian matrix to raise the accuracy.
I do not know too much numerical calculation of stiff odes, and here I want to know, is the Jacobian matrix JJ calculated in this manner:

If not, how to arrive at a correct one?
Many thanks! Hi everyone,
I am using ode15s to solve a set of stiff odes and the final model is developed in the form of
where M(t,y) is the mass matrix, and y is the generalised coordinate vector. According to the syntax of ode15s, it is crucial to provide its Jacobian matrix to raise the accuracy.
I do not know too much numerical calculation of stiff odes, and here I want to know, is the Jacobian matrix JJ calculated in this manner:

If not, how to arrive at a correct one?
Many thanks! ode15s, jacobian matrix MATLAB Answers — New Questions

​