The post-Newtonian approximation of the general relativistic Liouville's eq
uation is presented. Two integrals of motion, generalizations of the classi
cal energy and angular momentum, are obtained. Polytropic models are constr
ucted as an application.