Expand[InterpolatingPolynomial[{0, 20, 26, 31, 40}, 6]]