We apply the Hartle formalism to study equilibrium configurations in the framework of Newtonian gravity. This approach allows one to study in a simple manner the properties of the interior gravitational field in the case of static as well as stationary rotating stars in hydrostatic equilibrium. It is shown that the gravitational equilibrium conditions reduce to a system of ordinary differential equations which can be integrated numerically. We derive all the relevant equations up to the second order in the angular velocity. Moreover, we find explicitly the total mass, the moment of inertia, the quadrupole moment, the polar and equatorial radii, and the eccentricity of the rotating body. We also present the procedure to calculate the gravitational Love number. We test the formalism in the case of white dwarfs and show its compatibility with the known results in the literature.

We construct mass-radius relations of white dwarfs taking into account the effects of rotation and finite temperatures. We compare and contrast the theoretical mass-radius relations with observational data.

Recently, it has been shown by S.~M. de Carvalho et al. (2014) that the deviations between the degenerate case and observations were already evident for 0.7-0.8 M$_{\odot}$ white dwarfs. Such deviations were related to the neglected effects of finite temperatures on the structure of a white dwarf. Therefore, in this work by employing the Chandrasekhar equation of state taking into account the effects of temperature we show how the total pressure of the white dwarf matter depends on the mass density at different temperatures. Afterwards we construct equilibrium configurations of white dwarfs at finite temperatures. We obtain the mass-radius relations of white dwarfs for different temperatures by solving the Tolman-Oppenheimer-Volkoff equation, and compare them with the estimated masses and radii inferred from the Sloan Digital Sky Survey Data Release 4.