Abstract
The dynamical simulation method of the Lagrange's top proposed in [1] can be applied to nonsymmetric Euler's top in free motion. In application of the method, the invariant angular momentum is used to derive non-dimensional equations. To cancel integral error which central difference scheme accumulates, total derivative of the first integral is applied. For the evaluation of numerical precision, analytical solution of a symmetric Euler's top is utilized. The study case of the simulation is free rotation around near z-axis which is stable.
