  A self-consistent field method is developed, which can be used to construct
models of differentially rotating stars to first post-Newtonian order. The
rotation law is specified by the specific angular momentum distribution j(m),
where m is the baryonic mass fraction inside the surface of constant specific
angular momentum. The method is then used to compute models of the nascent
neutron stars resulting from the accretion induced collapse of white dwarfs.
The result shows that the ratios of kinetic energy to gravitational binding
energy of the relativistic models are slightly smaller than the corresponding
values of the Newtonian models.
