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.