The system of coupled nonlinear differential equations for a heavy symmetric top:
The solutions can be expressed through Weierstrass sigma and zeta functions:
Numerically check the correctness of the solutions:
Form any elliptic function with given periods, poles, and zeros as a rational function of
WeierstrassSigma:
Form an elliptic function with a single and a double zero and a triple pole:
Plot the resulting elliptic function: