InterpolatingPolynomial[{3,7,13,21,31, pi!},x]