Lagrange interpolation formula