By means of a piecewise cubic Hermitian polynomial approximation and a theorem for the asymptotics used in our previous papers, an accurate expression for the generalized Fermi-Dirac integral F ( z ) , derived from a Kane nonparabolic dispersion relation, for any values of reduced Fermi energy z and for the dimensionless conduction-band nonparabolicity parameter B in the interval [0,0.1493], is investigated with a precision of the order of 2 x lo-'. Then, two short series representations for F ( z ) are also presented and discussed. These accurate values of F ( z ) (or reduced carrier density u(z)) are used to evaluate the error of our two approximate expressions for z(u), appropriate to the lightly degenerate and degenerate cases, z 5 4 and z 2 4, respectively, and accurate within 0.5%.Finally, taking account of both the electron-electron interaction effect and the spin-orbit splitting interaction effect, the reduced Fermi energy z(u) and the electron effective mass in n-type heavily doped GaAs crystals at 4.2 and 300K are determined. Their numerical results are also compared with existing experimental results.