The second order field equations for the interior of slowly and rigidly rotating stars are solved numerically for different classes of possible equations of state and these solutions are then matched to the general asymptotically flat axisymmetric vacuum metric to second order in the rotational parameter. For these solutions we find that the quadrupole moment differs from that of the Kerr metric. Further we consider the post-Minkowskian limit analytically.