We consider the Lagrange interpolation employed in the multilevel fast multipole algorithm (MLFMA) as part of our efforts to obtain faster and more efficient solutions for large problems of computational electromagnetics. For the translation operator, we present the choice of the parameters for optimal interpolation. For the aggregation and dissaggregation processes, we discuss the interpolation matrices and introduce an efficient way of improving the accuracy by employing the poles.