T :=
proc(h,m)
option remember;
local n;

    if m=1 then if h=0 then return 0 else return 1 fi fi;

    n := m-1;

    (n+1)/n*
    add(binomial(n,q)*T(h,q)*T(h-1,n-q+1), q=1..n);
end;

TT := (h,m)-> T(h,m)-T(h-1,m);