Lagrangian interpolation