Expand[InterpolatingPolynomial[{{1, 3}, {2, -1}, {3, 1}, {4, 2}}, x]]