Given the Jacobian matrix for the RHS of the system of ordinary differential equations is already computed symbolically using SymPy, I should be able to work out a method for computing the Jacobian matrix for the system of residual equations whose root is the set of optimal coefficients.
Given the Jacobian matrix for the RHS of the system of ordinary differential equations is already computed symbolically using SymPy, I should be able to work out a method for computing the Jacobian matrix for the system of residual equations whose root is the set of optimal coefficients.