{D[Power[x, n], x], DifferenceDelta[FactorialPower[x, n], x]} // TraditionalForm
{Integrate[Power[x, n], x], Sum[FactorialPower[x, n], x]} // TraditionalForm
{Sum[(D[1 + 2 x + 3 x^2, {x, k}] /. x -> x0) Power[x - x0, k]/k!, {k, 0, 2}], Sum[(DifferenceDelta[1 + 2 x + 3 x^2, {x, k}] /. x -> x0) FactorialPower[x - x0, k]/k!, {k, 0, 2}]}
{Sum[(D[1 + 2 x + 3 x^2, {x, k}] /. x -> x0) Power[x - x0, k]/k!, {k, 0, 2}], Sum[(DifferenceDelta[1 + 2 x + 3 x^2, {x, k}] /. x -> x0) FactorialPower[x - x0, k]/k!, {k, 0, 2}]}; First[%] - Last[%] // FunctionExpand // Expand
Plot[Evaluate@Table[FactorialPower[x, k], {k, 1, 4}], {x, 0, 5}, PlotRange -> {-1.1, 2}]