(PARI) c(n, q) = prod(j=1, n-1, (q^(j+1) -1)/(q-1));
T(n, k, q) = c(k, q) + c(n-k, q) - c(n, q);
for(n=0, 10, for(k=0, n, print1(T(n, k, 3), ", "))) \\
G. C. Greubel, Dec 05 2019
(Magma) c:= func< n, q | n lt 2 select 1 else &*[(q^(j+1) -1)/(q-1): j in [1..n-1]] >;
T:= func< n, k, q | c(k, q) + c(n-k, q) - c(n, q) >;
[T(n, k, 3): k in [0..n], n in [0..10]]; //
G. C. Greubel, Dec 05 2019
(SageMath)
def c(n, q): return product( (q^(j+1) -1)/(q-1) for j in (1..n-1))
def T(n, k, q): return c(k, q) + c(n-k, q) - c(n, q)
[[T(n, k, 3) for k in (0..n)] for n in (0..10)] #
G. C. Greubel, Dec 05 2019