Join[{1, 4, 48}, LinearRecurrence[{7, -21, 35, -35, 21, -7, 1}, {240, 1000, 3300, 9114, 21952, 47520, 94500}, 20]]
Table[Piecewise[{{1, n == 1}, {48, n == 3}}, 6 Binomial[n, 3]^2 + n^2 (n - 1)^2], {n, 20}]
CoefficientList[Series[(-1 + 3 x - 41 x^2 + 47 x^3 - 223 x^4 + 221 x^5 - 217 x^6 + 127 x^7 - 42 x^8 + 6 x^9)/(-1 + x)^7, {x, 0, 20}], x]