The theorem states that the map sending a function to its interpolation polynomial is linear. It follows because interpolants are determined by their node values, so linear combinations of functions give the corresponding linear combinations of interpolating polynomials.