%I #25 Dec 31 2023 10:22:58
%S 1,3,7,14,25,40,56,63,37,-71,-350,-945,-2064,-3952,-6783,-10381,
%T -13625,-13330,-2359,33208,117672,288959,598325,1099385,1812546,
%U 2640543,3197152,2497824,-1541375,-12816925,-37865849,-86422322,-170718343,-301444536,-476474600,-655816385,-713055419,-351058887,1028750562,4501424879,11797832400,25361896880,47988600961
%N a(n) = 3*a(n-1) - 2*a(n-2) - a(n-3), with a(0)=1, a(1)=3 and a(2)=7.
%C Whittaker's Root Series Formula is applied to the polynomial equation -1+2x+3x^2+x^3. The following infinite series involving the Plastic Ratio (rho) is obtained: rho - 1 = 1/2 - 3/(2*7) + 7/(7*21) - 14/(21*65) + 25/(65*200) - 40/(200*616) + 56/(616*1897) - ...
%C The terms of the sequence appear in the numerators of the infinite sequence (with alternating signs). The denominators of the sequence are formed by multiplying consecutive terms from the sequence A218836.
%H E. T. Whittaker and G. Robinson, <a href="https://archive.org/details/calculusofobserv031400mbp/page/n139/mode/2up">The Calculus of Observations, London: Blackie & Son, Ltd. 1924, pp. 120-123.
%F a(n) = 3*a(n-1) - 2*a(n-2) - a(n-3).
%F a(n) = determinant of the n X n Toeplitz Matrix((3,2,-1,0,0,...,0),(3,1,0,0,0,...,0)).
%e a(0) = 1,
%e a(1) = 3*a(0) = 3*1 = 3,
%e a(2) = 3*a(1) - 2*a(0) = 3*3 - 2*1 = 7,
%e a(3) = 3*a(2) - 2*a(1) - a(0) = 3*7 - 2*3 - 1 = 14.
%p a:=proc(n) local c1,c2,c3;
%p option remember;
%p c1:=3; c2:=2; c3:=1;
%p if n=0 then 1
%p elif n=1 then 3
%p elif n=2 then 7
%p else c1*a(n-1)-c2*a(n-2)-c3*a(n-3); fi;
%p end; # _N. J. A. Sloane_, Dec 31 2023
%p [seq(a(n),n=0..30)];
%Y Cf. A218836 (denominator), A060006.
%K sign
%O 0,2
%A _Raul Prisacariu_, Dec 18 2023
|