A systematric method of exploring the 'geometrical' and 'non-geometrical' constants of the motion for an arbitrary spacetime is presented. This is done by introducing a series of coupled differential equation for the generators of the symmetry group of Vlasov's equation. The method is applied to the case of the maximaly symmetric spectime, and the geometrical and non-geometrical constants of motion are obtained.