We obtain an approximate global stationary and axisymmetric solution of Einstein's equations which can be thought of as a simple star model: a self-gravitating perfect fluid ball with a differential rotation motion pattern. Using the post-Minkowskian formalism (weak-field approximation) and considering rotation as a perturbation (slow-rotation approximation), we find approximate interior and exterior (asymptotically flat) solutions to this problem in harmonic coordinates. Interior and exterior solutions are matched, in the sense described by Lichnerowicz, on the surface of zero pressure, to obtain a global solution. The resulting metric depends on four arbitrary constants: mass density; rotational velocity at r = 0; a parameter that accounts for the change in rotational velocity through the star; and the star radius in the non-rotation limit. The mass, angular momentum, quadrupole moment and other constants of the exterior metric are determined in terms of these four parameters. PACS number(s) 04.40.Nr, 04.20.Jb