I would like to report what is not (to me ) a major problem, but a lack of reality . As reported in many reports, and in the great magazine "Nature" (here), the 4 minor satellites of Pluto have chaotic rotation, what is the result of the gravity variation generated by the duet "Pluto-Charon".
I suggest to change the rotation axis (I assume this is the "z" rotating axis) from still (immobile) compared to the "x" and "y" axis, and make them rotating with a Pi-factor on each of them the create a illusion of chaotic rotation. My english is poor so i made a drawing to explain:
- I assume the current "X" direction is a vectorial speed, calculated with some tangential functions to take into considerations parameters like elliptic orbit and variations of the acceleration.
- n1 and n2 should be absolutely different, and values aroud 1e-4°/s to 1e-3°/s
If i had any programing capabilities, i could help, but....