%I #8 Jun 13 2015 00:52:07
%S 1,6,25,116,517,2338,10517,47400,213481,961726,4332145,19515036,
%T 87908397,395998298,1783838637,8035595600,36197658961,163058307446,
%U 734522939465,3308779311556,14904940203477,67141752851858
%N The (1,3)-entry in the matrix M^n, where M is the 3 X 3 matrix [0,2,1; 2,1,2; 1,2,2] (n>=1).
%C a(n)/a(n-1) tends to 4.50466435...an eigenvalue of M and a root to the characteristic polynomial x^3 - 3x^2 - 7x + 1.
%H Title?, <a href="http://www.maths.utas.edu.au/People/Jackson/tables/jezek.html">Title?</a>
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (3,7,-1).
%F a(n)=3a(n-1)+7a(n-2)-a(n-3) (follows from the minimal polynomial of the matrix M).
%F G.f. x*(1+3*x) / ( 1-3*x-7*x^2+x^3 ). - _R. J. Mathar_, Mar 03 2013
%e a(7)=10517 because M^7= [6682,9842,10517;9842,14401,15438;10517,15438,16524].
%p with(linalg): M[1]:=matrix(3,3,[0,2,1,2,1,2,1,2,2]): for n from 2 to 25 do M[n]:=multiply(M[1],M[n-1]) od: seq(M[n][3,1],n=1..25);
%Y Cf. A120757.
%K nonn
%O 1,2
%A _Gary W. Adamson_ & _Roger L. Bagula_, Jul 01 2006
%E Corrected by _T. D. Noe_, Nov 07 2006
%E Edited by _N. J. A. Sloane_, Dec 04 2006
|