Вложенные суммы в SymPy

Question

Очевидно, что SymPy не может вычислить это, давая в результате NaN

t=Symbol('t')
s=Symbol('s')
a=Symbol('a')
k=Symbol('k')
k=2
Sum(binomial(k-t+s,a), (a,s+1,k-1), (s,0,k-2), (t,0,k-1) ).doit()

Что я делаю не так?