nmax := 14; mmax := nmax: for n from 1 to nmax do BETA(n, n) := 0 end do: m := 1: for n from m+1 to nmax do BETA(n, m) := (2*n-3)^2*BETA(n-1, m)-(2*n-4)! od: for m from 2 to mmax do for n from m+1 to nmax do BETA(n, m) := (2*n-3)^2*BETA(n-1, m) - BETA(n-1, m-1) od: od: for n from 2 to nmax do s1(n) := 0: for m from 1 to n-1 do s1(n) := s1(n) + BETA(n, m) od: od: seq(s1(n), n=2..nmax);
# Alternative:
nmax := nmax;
A120778 := proc(n): numer(sum(binomial(2*k1, k1)/(k1+1) / 4^k1, k1=0..n)) end proc:
A000165 := proc(n): 2^n*n! end proc:
A049606 := proc(n): denom(2^n/n!) end proc: for n from 2 to nmax do s2(n) := (-1)*
A120778(n-2)*
A000165(n-2)*
A049606(n-1) end do: seq(s2(n), n=2..nmax);