  Approximate analytical formula for density distribution in differentially
rotating stars is derived. Any barotropic EOS and conservative rotation law can
be handled with use of this method for wide range of differential rotation
strength. Results are in good qualitative agreement with comparison to the
other methods. Some applications are suggested and possible improvements of the
formula are discussed.
