|
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}]
|