read `SternCF.txt`:


LIMIT1:=1000:

gu:=[]:

R:=20:
print(`This is a table of  J[r]^(k)(1,x) for SYMBOLIC k, defined in Richard Stanley's paper arXiv:2101.02131v2`):

print(`Let F^(k)_i be the k-bonacci numbers and let `):
print(`Prod( 1+t*x^F^(k)_{i+k-1), i=1..n)= Sum(a(n,k)*x^k,k=0..infinity)`):
print(`and let v_r(n,t) be the sum of the r-th power of the coefficients `):
print(``):
print(`Let J[r,k] be the (ordinary) generating function with respect to x (no relation to the previous x)`):
print(``):

for r from 2 to R do
 lu:=JrkS(r,k,x,LIMIT1,3):
 gu:=[op(gu),lu]:
 print(J[r,k](1,x), `equals `):
 print(``):
 print(lu):
 print(``):
 print(`and in Maple format `):
 print(``):
 lprint(lu):
 print(``):
od:
print(`----------------------------`):
print(`To sum-up, the first`, R-1, `terms starting at r=2 are`):

lprint(gu):
print(``):

print(`----------------------------`):
print(``):
print(`This took`, time(), `seconds. `):




quit: