OFFSET
0,2
MATHEMATICA
(* Solve[x^3-x-1==0, x] *)
Phi=1.3247179572447460259609088544780973407;
Digits=18;
a=Table[If[m==0, Phi, Phi*Product[Prime[n], {n, 1, m}]], {m, 0, Digits}];
b=Table[Exp[Sum[n, {n, 0, m}]], {m, 0, Digits}];
c=Table[Floor[a[[Floor[(n+1)/2]]]*b[[Floor[n/2]]]], {n, 2, Digits+1}]
CROSSREFS
KEYWORD
nonn,uned,less
AUTHOR
Roger L. Bagula, Nov 20 2003
STATUS
approved