An efficient method for the analytic evaluation of the plasma dispersion function for the Fermi-Dirac distribution is proposed. The new method has been developed using the binomial expansion theorem and the Gamma functions. The general formulas obtained for the plasma dispersion function are utilized for the evaluation of the response function. The resulting series present better convergence rates. Several acceleration techniques are combined to further improve the efficiency. The obtained results for the plasma dispersion function are in good agreement with the known numerical data.