cp's OEIS Frontend

A089077 A nonsense sequence.

%I A089077 #15 Dec 18 2022 08:38:49
%S A089077 0,0,1,0,0,-1,0,-1,-1,-2,-2,-1,-2,0,-2,2,0,1,2,1,4,-1,3,-4,-1,-5,-5,
%T A089077 -4,-8,-1,-8,4,-4,7,4,7,11,3,14,-5,10,-11,-1,-15,-16,-11,-26,0,-25,15,
%U A089077 -10,25,15,24,40,9,49,-16,33,-41,-7,-50,-57,-34,-90,6,-83,56,-27,89,62,82,145,26,171,-63,108,-146,-37,-172,-209,-109,-316
%N A089077 A nonsense sequence.
%t A089077 (* Editors' note: this is not a valid program. *)
%t A089077 (* Adamson's matrix functions alternating x^4-x^3-x^2-x-1 Pisot*)
%t A089077 (* and x^4-x^3-1 minimal Pisot theta1*)
%t A089077 digits=200
%t A089077 Solve[x^4-x^3-1==0, x]
%t A089077 k=theta1 real root
%t A089077 q=N[k-1/k^3, 20]
%t A089077 m0={{0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 1}, {1, 0, 0, -q}}
%t A089077 NSolve[x^4-x^3-x^2-x-1==0, x]
%t A089077 k1=1.9275619754829254
%t A089077 q1=k1^2-k1-1/k1-1/k1^2
%t A089077 m1={{0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 1}, {1, 1, 1, q}}
%t A089077 m[n_Integer?Positive] := If[Mod[n, 2]==0, m[n-1].m0, m[n-1].m1]
%t A089077 m[0] ={{1, 0, 0, 0}, {0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 1}}
%t A089077 a=Table[Floor[m[n][[4, 4]]], {n, 1, digits}]
%K A089077 sign,uned,less
%O A089077 1,10
%A A089077 _Roger L. Bagula_, Dec 04 2003