Yeah, I think there was work on deriving the sensitivity equations on a symbolic level, but as far as I know it is not done yet. See
You might already tried it, but you can get the adjoint equations with this package:
(The adjoint ODE computes the derivative of a loss with respect to parameters. It’s mathematically kind of the ‘symbolic derivative’, which you can get to compute this gradient.)