Gauss equations and partial derivatives

Greetings to the Orekit Community,

I am looking for a way to evaluate Gauss’ equations, a.k.a the partial derivatives of orbital elements (ideally SMA / INC / RAAN / True longitude / E_X (e cos (w)) / E_Y (e sin(w)) relative to velocity expressed in the local orbital frame, as well as the partial derivative of the previous relative to orbital elements.

Below is the symbolic depiction of what I am after.

Can these quantities be evaluated by Orekit provided osculating orbital elements or cartesian position & velocities of the state of interest and retrieved in Python ?

Thanks !

To my short experience in orekit, No, you can not directly find those partial derivatives. But with the help of hipparchus’s DerivativeStructure, which is the underlying math library of orekit, Yes, you can develop code with Python to evaluate the quantities.