A215510 -id:A215510 - cp's OEIS frontend

A108716 a(n) = tan(Pi/14)^(-2n) + tan(3Pi/14)^(-2n) + tan(5Pi/14)^(-2n).

3, 21, 371, 7077, 135779, 2606261, 50028755, 960335173, 18434276035, 353858266965, 6792546291251, 130387472704741, 2502874814474531, 48044357383337973, 922243598852422035, 17703083191185355397

Offset: 0

Philippe Deléham, Jun 20 2005

The Berndt-type sequence number 11 for the argument 2*Pi/7 defined by the relation a(n) = t(1)^(2*n) + t(2)^(2*n) + t(4)^(2*n) = (-sqrt(7) + 4*s(1))^(2*n) + (-sqrt(7) + 4*s(2))^(2*n) + (-sqrt(7) + 4*s(4))^(2*n), where t(j) = tan(2*Pi*j/7) and s(j) = sin(2*Pi*j/7) (the respective sum with odd powers are discussed in A215794). See also A215007, A215008, A215143, A215493, A215494, A215510, A215512, A215694, A215695, A215828 and especially A215575, where a(n) = B(2n) for the function B(n) defined in the comments. - Roman Witula, Aug 23 2012
The sequence a(n+1)/a(n) is increasing and convergent to (t(2))^2 = 19,195669... (we note that the sequence A215794(n+1)/A215794(n) is decreasing and converges to the same limit). - Roman Witula, Aug 24 2012
Let L(p) be the total length of all sides and diagonals of a regular p-sided polygon inscribed in a unit circle. Then (L(p)/p)^2 = cot(Pi/(2p))^2 is the largest root of the equation: C(p,k)-C(p,2+k)*x+C(p,4+k)*x^2-C(p,6+k)*x^3+ ... +(-1)^q*x^q = 0, where k=1 if p is odd, k=0 if p is even, q = floor(p/2), and where C denotes the binomial coefficient. The complete set of roots is: x(i) = cot((2*i-1)*Pi/(2p))^2, i=1,2,...,q. Then a(n) = x(1)^n+x(2)^n+...x(q)^n for p=7. - Seppo Mustonen, Mar 25 2014
Sum_{k=1..(m-1)/2} tan^(2n) (k*Pi/m) is an integer when m >= 3 is an odd integer (see AMM link and formula); this sequence is the particular case m = 7. All terms are odd. - Bernard Schott, Apr 22 2022

Robert Israel, Table of n, a(n) for n = 0..700
Michel Bataille and Li Zhou, A Combinatorial Sum Goes on Tangent, The American Mathematical Monthly, Vol. 112, No. 7 (Aug. - Sep., 2005), Problem 11044, pp. 657-659.
Seppo Mustonen, Lengths of edges and diagonals and sums of them in regular polygons as roots of algebraic equations (2013).
Seppo Mustonen, Lengths of edges and diagonals and sums of them in regular polygons as roots of algebraic equations [Local copy]
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2Pi/7, J. Integer Seq., 12 (2009), Article 09.8.5.
Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6
Index entries for linear recurrences with constant coefficients, signature (21,-35,7).

Similar to: A000244 (m=3), 2*A165225 (m=5), this sequence (m=7), A353410 (m=9), A275546 (m=11), A353411 (m=13).

A:= gfun:-rectoproc({-a(n+3)+21*a(n+2)-35*a(n+1)+7*a(n), a(0) = 3, a(1) = 21, a(2) = 371},a(n), remember):
seq(A(n),n=0..20); # Robert Israel, Aug 23 2015

Table[ Round[ Cot[Pi/14]^(2n) + Cot[3Pi/14]^(2n) + Cot[5Pi/14]^(2n)], {n, 0, 12}] (* Robert G. Wilson v, Jun 21 2005 *)
RecurrenceTable[{a[0]== 3, a[1]== 21, a[2]==371, a[n]== 21*a[n-1] - 35*a[n-2] + 7*a[n-3]}, a, {n,30}] (* G. C. Greubel, Aug 22 2015 *)

a(n)=round(tan(Pi/14)^(-2*n) + tan(3*Pi/14)^(-2*n) + tan(5*Pi/14)^(-2*n)); \\ Anders Hellström, Aug 22 2015

a(n) = 7^n*A(2n), where A(n) := A(n-1) + A(n-2) + A(n-3)/7, with A(0)=3, A(1)=1, and A(2)=3. - see Witula-Slota's (Section 6) and Witula's (Remark 11) papers for the proofs and details. In these papers A(n) denotes the value of the big omega function with index n for the argument 2*i/sqrt(7) (see also A215512). - Roman Witula, Aug 23 2012
Conjecture: a(n) = 21*a(n-1)-35*a(n-2)+7*a(n-3). G.f.: -(35*x^2-42*x+3) / (7*x^3-35*x^2+21*x-1). - Colin Barker, Jun 01 2013
To verify conjecture, note that the roots of 7*x^3-35*x^2+21*x-1 are tan(Pi/14)^2, tan(3*Pi/14)^2 and tan(5*Pi/14)^2. - Robert Israel, Aug 23 2015
E.g.f.: exp((tan(Pi/7))^2*x) + exp((cot(Pi/14))^2*x) + exp((cot(3*Pi/14))^2*x). - G. C. Greubel, Aug 22 2015
a(n) = A275195(2*n)/(7^n). - Kai Wang, Aug 02 2016
a(n) = (tan(1*Pi/7))^(2*n) + (tan(2*Pi/7))^(2*n) + (tan(3*Pi/7))^(2*n). - Bernard Schott, Apr 22 2022

More terms from Robert G. Wilson v, Jun 21 2005

A033304 Expansion of (2 + 2x - 3x^2) / (1 - 2*x - x^2 + x^3).

Original entry on oeis.org

2, 6, 11, 26, 57, 129, 289, 650, 1460, 3281, 7372, 16565, 37221, 83635, 187926, 422266, 948823, 2131986, 4790529, 10764221, 24186985, 54347662, 122118088, 274396853, 616564132, 1385407029, 3112981337, 6994805571, 15717185450

Offset: 0

N. J. A. Sloane

From L. Edson Jeffery, Mar 22 2011: (Start)
Let A be the unit-primitive matrix (see [Jeffery])
A=A_(7,2)=
(0 0 1)
(0 1 1)
(1 1 1).
Let B={b(n)} be this sequence shifted to the right one place and setting b(0)=3. Then B=(3,2,6,11,26,...) with generating function (3-4*x-x^2)/(1-2*x-x^2+x^3) and b(n)=Trace(A^n). (End)
The following identity hold true (a(n)^2 - a(2n+2))/2 = A094648(n+1) = (-1)^(n+1)*A096975(n+1) - for the proof see Witula et al.'s papers - Roman Witula, Jul 25 2012
We note that the joined sequences (-1)^(n+1)*a(n) and A094648(n) form a two-sided sequence defined either by the recurrence formula x(n+3) + x(n+2) - 2x(n+1) - x(n) = 0, n in Z, x(0)=3, x(-1)=-2, x(1)=-1, or by the following trigonometric identities: x(n) = (c(1))^n + (c(2))^n + (c(4))^n = (c(1)c(2))^(-n) + (c(1)c(4))^(-n) + (c(2)c(4))^(-n) = (s(2)/s(1))^n + (s(4)/s(2))^n + (s(1)/s(4))^n, for n in Z, where c(j) := 2*cos(2Pi*j/7) and s(j) := sin(2*Pi*j/7) - for the proof see Witula's and Witula et al.'s papers. - Roman Witula, Jul 25 2012
We have 4*a(n+2) - a(n) = 7*A077998(n+2). - Roman Witula, Aug 13 2012
Two very intriguing identities of trigonometric nature hold: (-1)^n*(a(n)-a(n-1)) = c(1)*c(2)^(-n) + c(2)*c(4)^(-n) + c(4)*c(1)^(-n), and (-1)^(n+1)*(a(n-1)-a(n+1)) = c(1)*c(4)^(-n-1) + c(2)*c(1)^(-n-1) + c(4)*c(2)^(-n-1), where a(-1):=3 and c(j) is defined as above. For the proof see Remark 6 in the first Witula's paper. - Roman Witula, Aug 14 2012
With respect to the form of the trigonometric formulas describing a(n), we call this sequence the Berndt-type sequence number 20 for the argument 2Pi/7. The A-numbers of other Berndt-type sequences numbers are given in below. - Roman Witula, Sep 30 2012

R. P. Stanley, Enumerative Combinatorics I, p. 244, Eq. (36).

G. C. Greubel, Table of n, a(n) for n = 0..1000
Mark W. Coffey, James L. Hindmarsh, Matthew C. Lettington, John Pryce, On Higher Dimensional Interlacing Fibonacci Sequences, Continued Fractions and Chebyshev Polynomials, arXiv:1502.03085 [math.NT], 2015 (see p. 40).
L. E. Jeffery, Unit-primitive matrices
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2Pi/7, J. Integer Seq., 12 (2009), Article 09.8.5.
Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6
Roman Witula, Damian Slota and Adam Warzynski, Quasi-Fibonacci Numbers of the Seventh Order, J. Integer Seq., 9 (2006), Article 06.4.3.
Index entries for linear recurrences with constant coefficients, signature (2,1,-1).

Cf. A066170, A006356, A096975.

Cf. A215007, A215008, A215143, A215493, A215494, A215510, A215512, A215575, A215694, A215695, A108716, A215794, A215828, A215817, A215877, A094429, A094430, A217274, A094648.

I:=[2,6,11]; [n le 3 select I[n] else 2*Self(n-1) +Self(n-2) - Self(n-3): n in [1..30]]; // G. C. Greubel, Apr 19 2018

CoefficientList[Series[(2+2x-3x^2)/(1-2x-x^2+x^3),{x,0,50}], x]  (* Harvey P. Dale, Mar 14 2011 *)
LinearRecurrence[{2, 1, -1}, {2, 6, 11}, 29] (* Jean-François Alcover, Sep 27 2017 *)

{a(n)=if(n<0, n=-n; polsym(x^3-x^2-2*x+1,n-1)[n], n+=2; polsym(1-x-2*x^2+x^3,n-1)[n])} /* Michael Somos, Aug 03 2006 */

x='x+O('x^99); Vec((2+2*x-3*x^2)/(1-2*x-x^2+x^3)) \\ Altug Alkan, Apr 19 2018

a(-1-n) = A096975(n).
a(n) = (1-2*cos(1/7*Pi))^(n+1)+(1+2*cos(2/7*Pi))^(n+1)+(1-2*cos(3/7*Pi))^(n+1). - Vladeta Jovovic, Jun 27 2001
a(n) = trace of (n+1)-th power of the 3 X 3 matrix (in the example of A066170): [1 1 1 / 1 1 0 / 1 0 0]. Alternatively, the sum of the (n+1)st powers of the roots of the corresponding characteristic polynomial: x^3 - 2*x^2 - x + 1 = 0. a(n) = A006356(n) + A006356(n-1) + 2*A006356(n-2). E.g., a(3) = 26 = the trace of M^4. The characteristic polynomial of this matrix (see A066170) is x^3 - 2*x^2 - x + 1 and the roots are 2.24697960372..., -0.8019377358... and 0.55495813208... = a, b, c. Then Sum(a^4 + b^4 + c^4) = 26. - Gary W. Adamson, Feb 01 2004
(-1)^(n+1)*a(n) = (c(1))^(-n-1) + (c(2))^(-n-1) + (c(3))^(-n-1) = (c(1)c(2))^(n+1) + (c(1)c(4))^(n+1) + (c(2)c(4))^(n+1) = (s(1)/s(2))^(n+1) + (s(2)/s(4))^(n+1) + (s(4)/s(1))^(n+1), where c(j) := 2*cos(2*Pi*j/7) and s(j) := sin(2*Pi*j/7) - for the proof see Witula's and Witula et al.'s papers. - Roman Witula, Jul 25 2012
a(n) = 3*A077998(n+1) - A006054(n+2) - A006054(n+1). - Roman Witula, Aug 13 2012
a(n)*(-1)^(n+1) = (A094648(n+1)^2 - A094648(2*(n+1)))/2. - Roman Witula, Sep 30 2012

A094648 An accelerator sequence for Catalan's constant.

Original entry on oeis.org

3, -1, 5, -4, 13, -16, 38, -57, 117, -193, 370, -639, 1186, -2094, 3827, -6829, 12389, -22220, 40169, -72220, 130338, -234609, 423065, -761945, 1373466, -2474291, 4459278, -8034394, 14478659, -26088169, 47011093, -84708772, 152642789, -275049240

Offset: 0

Paul Barry, May 18 2004

The pair A094648 and the alternating sequence A033304 when joined form a two-sided sequence defined by the recurrence formula x(n+3) + x(n+2) - 2x(n+1) - x(n) = 0, n in Z, x(-1)=-2, x(0)=3, x(1)=-1 - for details see Witula's comments to A033304. - Roman Witula, Jul 25 2012
From Roman Witula, Aug 09 2012: (Start)
There exist two interesting subsequences b(n) and c(n) of the given above sequence x(n) defined by the following relations: b(n)=a(2^n) and c(n)=x(-2^n). These subsequences satisfy the following system of recurrence equations:
b(n+1)=b(n)^2-2*c(n), and c(n+1)=c(n)^2-2*b(n),
which easily follow from the general identity: x(n)^2=x(2*n)-2*x(-n), n in Z. We note that b(0)=-1, b(1)=5, b(2)=13, b(3)=117, c(0)=-2, c(1)=6, c(2)=26, c(3)=650. From the above system we deduce that all b(n) are odd, whereas all c(n) are even. Moreover we obtain c(n+1)-b(n+1)=(c(n)-b(n))*(b(n)+c(n)+2), which yields b(n+1)-c(n+1)=product{k=1,..,n}(b(k)+c(k)+2)=13*product{k=2,..,n}(b(k)+c(k)+2)=13^2*41*product{k=3,..,n}(b(k)+c(k)+2). It follows that b(n)-c(n) is divisible by 13^2*41 for every n=3,4,..., and after using the above system again each b(n) and c(n), for n=2,3,..., is divisible by 13. (End)
If we set W(n):=3*A077998(n)-A006054(n+1)-A006054(n), n=0,1,..., then a(n)=(W(n)^2-W(2*n))/2 and W(n) = (-c(1))^(-n) + (-c(2))^(-n) + (-c(4))^(-n) = (-c(1)*c(2))^n + (-c(1)*c(4))^n + (-c(2)*c(4))^n = (-1-c(1))^n + (-1-c(2))^n + (-1-c(4))^n, where c(j):=2*cos(2*Pi*j/7) - for the proof see Witula-Slota-Warzynski's paper. Moreover it follows from the comment at the top and from comments to A033304 that W(n+1)=A033304(n)=(-1)^(n+1)*x(-n-1). - Roman Witula, Aug 11 2012
The following trigonometric type identitities hold true: (1) -a(n-1)-a(n) = c(1)*c(2)^n + c(2)*c(4)^n + c(4)*c(1)^n and (2) a(n)-a(n+2) = c(4)*c(2)^(n+1) + c(1)*c(4)^(n+1) + c(2)*c(1)^(n+1), where a(-1)=-2 and c(j) is defined as above (see also the respective comment to A033304). For the proof see Remark 6 in Witula's paper. - Roman Witula, Aug 14 2012
It can be proved that A033304(n-1)*(-1)^n = (a(n)^2 - a(2*n))/2, n=1,2,... - Roman Witula, Sep 30 2012
With respect to the form of the trigonometric formulas describing a(n), we call this sequence the Berndt-type sequence number 19 for the argument 2*Pi/7. The A-numbers of other Berndt-type sequences numbers are given in below. - Roman Witula, Sep 30 2012

			We have a(17) = a(19) + 50000, a(4) + a(5) = -3, 2*a(7) + a(8) = 3, and 2*a(9) + a(10) = a(5). - _Roman Witula_, Sep 14 2012

G. C. Greubel, Table of n, a(n) for n = 0..3900
A. Akbary and Q. Wang, On some permutation polynomials over finite fields, International Journal of Mathematics and Mathematical Sciences, 2005:16 (2005) 2631-2640.
A. Akbary and Q. Wang, A generalized Lucas sequence and permutation binomials, Proceeding of the American Mathematical Society, 134 (1) (2006), 15-22.
David M. Bradley, A Class of Series Acceleration Formulae for Catalan's Constant, The Ramanujan Journal, Vol. 3, Issue 2, 1999, pp. 159-173
David M. Bradley, A Class of Series Acceleration Formulae for Catalan's Constant, arXiv:0706.0356 [math.CA], 2007.
Kai Wang, Fibonacci Numbers And Trigonometric Functions Outline, (2019).
Q. Wang, On generalized Lucas sequences, Contemp. Math. 531 (2010) 127-141
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2Pi/7, J. Integer Seq., 12 (2009), Article 09.8.5.
Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6
Roman Witula, Damian Slota and Adam Warzynski, Quasi-Fibonacci Numbers of the Seventh Order, J. Integer Seq., 9 (2006), Article 06.4.3.
Index entries for linear recurrences with constant coefficients, signature (-1,2,1).

Cf. A000032, A094649, A094650.

Cf. A215007, A215008, A215143, A215493, A215494, A215510, A215512, A215575, A215694, A215695, A108716, A215794, A215828, A215817, A215877, A094429, A094430, A217274.

I:=[3,-1,5]; [n le 3 select I[n]  else -Self(n-1)+2*Self(n-2)+Self(n-3): n in [1..40]]; // Vincenzo Librandi, Jul 25 2015

CoefficientList[ Series[(3 + 2x - 2x^2)/(1 + x - 2x^2 - x^3), {x, 0, 33}], x] (* Robert G. Wilson v, May 24 2004 *)
a[n_] := Round[(2Sin[3Pi/14])^n + (-2Sin[Pi/14])^n + (-2Cos[Pi/7])^n]; Table[ a[n], {n, 0, 33}] (* Robert G. Wilson v, May 24 2004 *)
LinearRecurrence[{-1,2,1}, {3,-1,5}, 50] (* Roman Witula, Aug 09 2012 *)

x='x+O('x^30); Vec((3+2*x-2*x^2)/(1+x-2*x^2-x^3)) \\ G. C. Greubel, May 09 2018

G.f.: (3+2*x-2*x^2)/(1+x-2*x^2-x^3);
a(n) = (2*sin(3*Pi/14))^n+(-2*sin(Pi/14))^n+(-2*cos(Pi/7))^n.
a(p) == -1 mod(p), p prime. - Philippe Deléham, Oct 03 2009
a(n) = (2*cos(2*Pi/7))^n + (2*cos(4*Pi/7))^n + (2*cos(8*Pi/7))^n, which is equivalent to the formula given above (for analogous sums with sines see A215493 and A215494). Moreover we have a(n+3) + a(n+2) - 2a(n+1) - a(n) = 0 - for the proof see Witula-Slota's paper. - Roman Witula, Jul 24 2012
a(n) = 3*(-1)^n*A006053(n+2) +2*A078038(n-1). - R. J. Mathar, Nov 03 2020

A215575 a(n) = 7*(a(n-1) - a(n-2) - a(n-3)), with a(0)=3, a(1)=7, a(2)=35.

Original entry on oeis.org

3, 7, 35, 175, 931, 5047, 27587, 151263, 830403, 4560871, 25054435, 137642127, 756187747, 4154438295, 22824258947, 125395430335, 688917131651, 3784882096583, 20793986742179, 114241312597615, 627637106311971, 3448212648805239, 18944339609269571

Offset: 0

Roman Witula, Aug 16 2012

The Berndt-type sequence number 8 for the argument 2Pi/7 defined by trigonometric relations from "Formula" below.
We note that the following decompositions hold true: (X-cot(2*Pi/7)^n)*(X-cot(4*Pi/7)^n)*(X-cot(8*Pi/7)^n) = X^3 - sqrt(7)^(-n)*a(n)*X^2 + (-sqrt(7))^(-n)*B(n)*X
  - (-sqrt(7))^(-n), and (X-tan(2*Pi/7)^n)*(X-tan(4*Pi/7)^n)*(X-tan(8*Pi/7)^n) = X^3 - B(n)*X^2  + (-1)^n*a(n)*X - (-sqrt(7))^n, where B(n) := tan(2*Pi/7)^n + tan(4*Pi/7)^n + tan(8*Pi/7)^n = (-sqrt(7) + 4*sin(2*Pi/7))^n + (-sqrt(7) + 4*sin(4*Pi/7))^n + (-sqrt(7) + 4*sin(8*Pi/7))^n. Moreover we have 2*(-1)^n*B(n) = 7^(-n/2)*(a(n)^2 - a(2*n)). For the proof of these decompositions see Witula-Slota's (Section 6) and Witula's (Remark 11) reference.
We note that the numbers a(n)*7^(-ceiling(n/3)) are all integers. Moreover from the recurrence relation: a(n+3)+7*a(n+1)=7*(a(n+2)-a(n)) it can be easily obtained the following summations formulas: 8*sum{k=1,..,n} a(2*k) = 7*(a(2*n+1)-2)-a(2*n+2), which also means that the result is divisible by 8 for every n=1,2,..., and 8*sum{k=1,..,n} a(2*k-1) = 7*(a(2*n)-2)-a(2*n+1), which implies that 7*(a(n)-2)-a(n+1) is divisible by 8 for each n=0,1,...

			We have  cot(2*Pi/7)^2 + cot(4*Pi/7)^2 + cot(8*Pi/7)^2 = 5,  cot(2*Pi/7)^4 + cot(4*Pi/7)^4 + cot(8*Pi/7)^4 = 19, but cot(2*Pi/7)^6 + cot(4*Pi/7)^6 + cot(8*Pi/7)^6 = 563/7. Similarly the numbers sqrt(7)*(cot(2*Pi/7)^n + cot(4*Pi/7)^n + cot(8*Pi/7)^n) are integers for n=1,3,5,7 (equal to 7, 25, 103, 441, respectively), whereas for n=9 we obtain the rational value 13297/7.

E Hetmaniok, P Lorenc, S Damian, et al., Periodic orbits of boundary logistic map and new kind of modified Chebyshev polynomials in R. Witula, D. Slota, W. Holubowski (eds.), Monograph on the Occasion of 100th Birthday Anniversary of Zygmunt Zahorski. Wydawnictwo Politechniki Slaskiej, Gliwice 2015, pp. 325-343.

G. C. Greubel, Table of n, a(n) for n = 0..1350
Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2*Pi/7, Journal of Integer Sequences, Vol. 12 (2009), Article 09.8.5.
R. Witula, P. Lorenc, M. Rozanski, M. Szweda, Sums of the rational powers of roots of cubic polynomials, Zeszyty Naukowe Politechniki Slaskiej, Seria: Matematyka Stosowana z. 4, Nr. kol. 1920, 2014.
Index entries for linear recurrences with constant coefficients, signature (7, -7, -7).

Cf. A215007, A215008, A215143, A215493, A215494, A215510, A215512, A215455.

I:=[3,7,35]; [n le 3 select I[n] else 7*Self(n-1) - 7*Self(n-2) - 7*Self(n-3): n in [1..30]]; // G. C. Greubel, Apr 25 2017

LinearRecurrence[{7,-7,-7}, {3,7,35}, 50]

polsym(x^3 - 7*x^2 + 7*x + 7, 30) \\ Charles R Greathouse IV, Jul 20 2016

a(n) = (sqrt(7)^n)*(cot(2*Pi/7)^n + cot(4*Pi/7)^n + cot(8*Pi/7)^n) = (3 + 4*cos(2*Pi/7))^n + (3 + 4*cos(4*Pi/7))^n + (3 + 4*cos(8*Pi/7))^n = (-tan(2*Pi/7)*tan(4*Pi/7))^n + (-tan(2*Pi/7)*tan(8*Pi/7))^n + (-tan(4*Pi/7)*tan(8*Pi/7))^n.
G.f.: (3-14*x+7*x^2)/(1-7*x+7*x^2+7*x^3).
a(n) = x1^n + x2^n + x3^n, where x1, x2, x3 are the roots of the polynomial x^3 - 7*x^2 + 7*x + 7, that is, x1 = sqrt(7)/tan(Pi/7), x2 = sqrt(7)/tan(2*Pi/7), x3 = sqrt(7)/tan(4*Pi/7). - Kai Wang, Jul 19 2016

A215494 a(n) = 7a(n-1) - 14a(n-2) + 7*a(n-3) with a(1)=7, a(2)=21, a(3)=70.

Original entry on oeis.org

7, 21, 70, 245, 882, 3234, 12005, 44933, 169099, 638666, 2417807, 9167018, 34790490, 132119827, 501941055, 1907443237, 7249766678, 27557748813, 104759610858, 398257159370, 1514069805269, 5756205681709, 21884262613787, 83201447389466, 316323894905207

Offset: 1

Roman Witula, Aug 13 2012

The Berndt-type sequence number 5 for the argument 2*Pi/7; see also A215007, A215008, A215143, A215493 and A215510.
We note that if we set:
  x(n) := s(2)*s(1)^n + s(4)*s(2)^n + s(1)*s(4)^n,
  y(n) := s(4)*s(1)^n + s(1)*s(2)^n + s(2)*s(4)^n,
  z(n) := s(1)^(n+1) + s(2)^(n+1) + s(4)^(n+1),
  for every n=0,1,..., where s(j) := 2*sin(2*Pi*j/7), then the following system of recurrence equations holds true:
  x(n+2)=2*x(n)-y(n), y(n+2)=2*y(n)-x(n)+z(n), z(n+2)=y(n)+3*z(n).
Moreover we have a(n)=z(2*n-1), A215493(n)=z(2*n-2), A094429(n)=y(2n-1)-x(2n-1)=-x(2*n+2)/sqrt(7), A094430(n)=-x(2*n+3), y(2*n-2)=sqrt(7)*A215143(n), y(2*n-1)=A215510(n) and x(11)=-(y(10)+z(10))/sqrt(7)=-1078.
We can also deduce the following relations:
  x(n-1) = c(1)*s(1)^n + c(2)*s(2)^n + c(4)*s(4)^n,
  -y(n-1)-z(n-1) = c(2)*s(1)^n + c(4)*s(2)^n + c(1)*s(4)^n,
  y(n-1)-x(n-1) = c(4)*s(1)^n + c(1)*s(2)^n + c(2)*s(4)^n,
  for every n=1,2,..., where x(0)=y(0)=z(0)=sqrt(7), and c(j) := 2*cos(2*Pi*j/7).
All these sequences satisfy the following recurrence equation: Z(n+6)-7*Z(n+4)+14*Z(n+2)-7*Z(n)=0. The characteristic polynomial of this equation (after rescaling) has the form (X-s(1)^2)*(X-s(2)^2)*(X-s(3)^2)=X^3-7*X^2+14*X-7 and was recognized by Johannes Kepler (1571-1630); see the Savio-Suryanarayan paper.
We also have the following decomposition: (X-s(1)^(n+1))*(X-s(2)^(n+1))*(X-s(4)^(n+1)) = X^3 - z(n)*X^2 + (1/2)*(z(n)^2-z(2n+1))*X - (-sqrt(7))^(n+1).
Further we have a(n)=A146533(n) for n=1,...,6, and A146533(7)-a(7)=7. We note that all numbers 7^(-1-floor(n/3))*a(n) are integers.

			We have a(3)=5*7^2 and a(6)=5*7^4, which implies that s(1)^12 + s(2)^12 + s(4)^12 = 49*(s(1)^6 + s(2)^6 + s(4)^6). We also have a(9) = (a(1) + a(3))*7^49.

G. C. Greubel, Table of n, a(n) for n = 1..1000
B. C. Berndt, A. Zaharescu, Finite trigonometric sums and class numbers, Math. Ann. 330 (2004), 551-575.
B. C. Berndt, L.-C. Zhang, Ramanujan's identities for eta-functions, Math. Ann. 292 (1992), 561-573.
L. Bankoff and J. Garfunkel, The Heptagonal Triangle, Math. Magazine, 46 (1973), 7-19.
Z.-G. Liu, Some Eisenstein series identities related to modular equations of the seventh order, Pacific J. Math. 209 (2003), 103-130.
D. Y. Savio and E. R. Suryanarayan, Chebyshev polynomials and regular polygons, Amer. Math. Monthly, 100 (1993), 657-661.
Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6
Index entries for linear recurrences with constant coefficients, signature (7, -14, 7).

See A122068.

I:=[7,21,70]; [n le 3 select I[n] else 7*Self(n-1)-14*Self(n-2)+7*Self(n-3): n in [1..30]]; // Vincenzo Librandi, Dec 01 2016

LinearRecurrence[{7,-14,7}, {7,21,70}, 50]

polsym(x^3 - 7*x^2 + 14*x - 7, 30) \\ (includes a(0)=3) Joerg Arndt, May 31 2017

x='x+O('x^30); Vec((7-28*x+21*x^2)/(1-7*x+14*x^2-7*x^3)) \\ G. C. Greubel, Apr 23 2018

Equals 7*A122068. - M. F. Hasler, Aug 25 2012
a(n) = s(1)^(2n) + s(2)^(2n) + s(4)^(2n), where s(j) := 2*Sin(2*Pi*j/7) (for the sums of the respective odd powers see A215493, see also A094429, A115146). For the proof of these formula see Witula-Slota's paper.
G.f.: (7 - 28*x + 21*x^2)/(1 - 7*x + 14*x^2 - 7*x^3) = -d(log(1 - 7*x + 14*x^2 - 7*x^3))/dx.

A215512 a(n) = 5a(n-1) - 6a(n-2) + a(n-3), with a(0)=1, a(1)=3, a(2)=8.

Original entry on oeis.org

1, 3, 8, 23, 70, 220, 703, 2265, 7327, 23748, 77043, 250054, 811760, 2635519, 8557089, 27784091, 90213440, 292919743, 951102166, 3088205812, 10027335807, 32558546329, 105716922615, 343260670908, 1114560365179, 3618954723062, 11750672095144, 38154192502527

Offset: 0

Roman Witula, Aug 14 2012

The Berndt-type sequence number 7 for the argument 2Pi/7 defined by the relation: sqrt(7)*a(n) = s(1)*c(4)^(2*n) + s(2)*c(1)^(2*n) + s(4)*c(2)^(2*n), where c(j):=2*cos(2*Pi*j/7) and s(j):=2*sin(2*Pi*j/7). If we additionally defined the following sequences:
sqrt(7)*b(n) = s(2)*c(4)^(2*n) + s(4)*c(1)^(2*n) + s(1)*c(2)^(2*n),
sqrt(7)*c(n) = s(4)*c(4)^(2*n) + s(1)*c(1)^(2*n) + s(2)*c(2)^(2*n), and
sqrt(7)*a1(n) = s(1)*c(4)^(2*n+1) + s(2)*c(1)^(2*n+1) + s(4)*c(2)^(2*n+1),
sqrt(7)*b1(n) = s(2)*c(4)^(2*n+1) + s(4)*c(1)^(2*n+1) + s(1)*c(2)^(2*n+1),
sqrt(7)*c1(n) = s(4)*c(4)^(2*n+1) + s(1)*c(1)^(2*n+1) + s(2)*c(2)^(2*n+1), then the following simple relationships between elements of these sequences hold true: a(n)=c1(n), c(n+1)=a1(n), -a(n)-b(n)=b1(n), which means that the sequences a1(n), b1(n), and c1(n) are completely and in very simple way determined by the sequences a(n), b(n) and c(n). However the last one's satisfy the following system of recurrence equations: a(n+1) = 2*a(n) + b(n), b(n+1) = a(n) + 2*b(n) - c(n), c(n+1) = c(n) - b(n). We have b(n)=A215694(n) and c(n)=A215695(n).
We note that a(n)=A000782(n) for every n=0,1,...,4 and A000782(5)-a(5)=2.
From general recurrence relation: a(n) = 5*a(n-1) - 6*a(n-2) + a(n-3), i.e. a(n) = 5*(a(n-1)-a(n-2)) + (a(n-3)-a(n-2)) the following summation formula can be easily obtained: sum{k=3,..,n} a(k) = 5*a(n-1)-a(n-2)+a(0)-5*a(1). Hence in discussed sequence it follows that: sum{k=3,..,n} a(k) = 5*a(n-1) - a(n-2) - 14.

			We have a(6) = 10*a(4)+a(1), a(5) = 11*(a(3)-a(1)), a(10)-a(4)+a(3)+a(1)+a(0) = 77*10^3, and a(11)-a(4)+a(3)-a(2)+a(0) = 25*10^4 = (5^6)*(2^4).

G. C. Greubel, Table of n, a(n) for n = 0..1000
Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6
Index entries for linear recurrences with constant coefficients, signature (5, -6, 1).

Cf.A215694, A215695, A215007, A215008, A215143, A215493, A215494, A215510, A215575, A215455, A214683, A214699.

I:=[1,3,8]; [n le 3 select I[n] else 5*Self(n-1) - 6*Self(n-2) + Self(n-3): n in [1..30]]; // G. C. Greubel, Apr 23 2018

Mathematica
```
LinearRecurrence[{5,-6,1}, {1,3,8}, 50]
```

x='x+O('x^30); Vec((1-2*x-x^2)/(1-5*x+6*x^2-x^3)) \\ G. C. Greubel, Apr 23 2018

G.f.: (1-2*x-x^2)/(1-5*x+6*x^2-x^3).

A094429 Given the 3 X 3 matrix M = [0 1 0 / 0 0 1 / 7 -14 7], a(n) = (-) rightmost term of M^n * [1 1 1].

Original entry on oeis.org

0, 7, 42, 196, 833, 3381, 13377, 52136, 201341, 773122, 2958032, 11291903, 43042727, 163918671, 623875840, 2373568575, 9028148962, 34334213564, 130560389505, 496440779373, 1887579497489, 7176808297736, 27286630574917

Offset: 1

Gary W. Adamson, May 02 2004

M is derived from the Lucas polynomial: x^3 - 7*x^2 + 14*x - 7 with a root (and eigenvalue of the matrix): 3.801377358... = (2*sin(3*Pi/7))^2, the convergent of the sequence.
From Roman Witula, Sep 29 2012: (Start)
The Berndt-type sequence number 16 for the argument 2*Pi/7 (see Formula section and Crossrefs for other Berndt-type sequences for the argument 2*Pi/7 - for numbers from 1 to 18 without 16).
Note that all numbers of the form a(n)*7^(-1 - floor((n-1)/3)) are integers and even a(10) and a(11) are divisible by 7^5. (End)

			a(5) = 833. M^5 * [1 1 1] = [ -42 -196 -833].
We have 4*a(4) - a(5) = 4*a(5) - a(6) = 7*a(2) = 49, 88*a(10) = 23*a(11), and a(3) = 6*a(2), which implies the equalities c(4)*(s(1))^6 + c(2)*(s(4))^6 + c(1)*(s(2))^6 = 6*(c(4)*(s(1))^4 + c(2)*(s(4))^4 + c(1)*(s(2))^4) and
s(2)*(s(1))^8 + s(4)*(s(2))^8 + s(1)*(s(4))^8 = 6*( s(2)*(s(1))^6 + s(4)*(s(2))^6 + s(1)*(s(4))^6). - _Roman Witula_, Sep 29 2012

G. C. Greubel, Table of n, a(n) for n = 1..1700
Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6.
Index entries for linear recurrences with constant coefficients, signature (7,-14,7).

Cf. A215007, A215008, A215143, A215493, A215494, A215510, A215512, A215575, A215694, A215695, A108716, A215794, A215828, A215817, A215877, A094430, A217274.

I:=[0,7,42]; [n le 3 select I[n] else 7*Self(n-1) -14*Self(n-2) +7*Self(n-3): n in [1..30]]; // G. C. Greubel, May 09 2018

Table[(MatrixPower[{{0, 1, 0}, {0, 0, 1}, {7, -14, 7}}, n].{-1, -1, -1})[[3]], {n, 23}] (* Robert G. Wilson v, May 08 2004 *)
LinearRecurrence[{7,-14,7}, {0,7,42}, 50] (* Roman Witula, Aug 13 2012 *)

x='x+O('x^30); concat([0], Vec(7*x^2*(1-x)/(1-7*x+14*x^2-7*x^3))) \\ G. C. Greubel, May 09 2018

a(n) = -(([0, 1, 0; 0, 0, 1; 7, -14, 7]^n)*[1,1,1]~)[3]; \\ Michel Marcus, May 10 2018

From Colin Barker, Jun 19 2012: (Start)
a(n) = 7*a(n-1) - 14*a(n-2) + 7*a(n-3).
G.f.: 7*x^2*(1-x)/(1 - 7*x + 14*x^2 - 7*x^3). (End)
a(n) = c(4)*(s(1))^(2*n) + c(2)*(s(4))^(2*n) + c(1)*(s(2))^(2*n) = (-1/sqrt(7))*(c(1)*(s(1))^(2*n+3) + c(2)*(s(2))^(2*n+3) + c(3)*(s(3))^(2*n+3)) = (-1/sqrt(7))*(s(2)*(s(1))^(2*n+2) + s(4)*(s(2))^(2*n+2) + s(1)*(s(4))^(2*n+2)), where c(j) := 2*cos(2*Pi*j/7) and s(j) := 2*sin(2*Pi*j/7) (for the sums of the respective odd powers see A217274, see also A215493 and comments to A215494). For the proof of these formulas see Witula-Slota's paper. - Roman Witula, Jul 24 2012

More terms from Robert G. Wilson v, May 08 2004

A094430 a(n) is the rightmost term of M^n * [1 0 0], where M is the 3 X 3 matrix [0 1 0 / 0 0 1 / 7 -14 7].

Original entry on oeis.org

7, 49, 245, 1078, 4459, 17836, 69972, 271313, 1044435, 4002467, 15294370, 58337097, 222255768, 846131608, 3219700183, 12247849145, 46582062709, 177142452214, 673583231587, 2561162729076, 9737971026812, 37024601601729

Offset: 1

Gary W. Adamson, May 02 2004

In A094429 the multiplier is [1 1 1] instead of [1 0 0]. The matrix M is derived from the 3rd-order Lucas polynomial x^3 - 7x^2 + 14x - 7, with a convergent of the series = 3.801937735... = (2 sin 3*Pi/7)^2; (an eigenvalue of the matrix and a root of the polynomial).
From Roman Witula, Sep 29 2012: (Start)
This sequence is the Berndt-type sequence number 17 for the argument 2*Pi/7 (see Formula section and Crossrefs for other Berndt-type sequences for the argument 2*Pi/7 - for numbers from 1 to 18 without 17).
Note that all numbers of the form a(n)*7^(-floor((n+4)/3)) are integers. (End)

			a(4) = 1078 since M^4 * [1 0 0] = [49 245 1078] = [a(2), a(3), a(4)].
We have a(2)=7*a(1), a(3)=5*a(2), 22*a(3)=5*a(4), and a(6)=4*a(5), which implies s(2)*s(1)^15 + s(4)*s(2)^15 + s(1)*s(4)^15 = 4*(s(2)*s(1)^13 + s(4)*s(2)^13 + s(1)*s(4)^13). - _Roman Witula_, Sep 29 2012

G. C. Greubel, Table of n, a(n) for n = 1..1000
Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6.
Index entries for linear recurrences with constant coefficients, signature (7,-14,7).

Cf. A215007, A215008, A215143, A215493, A215494, A215510, A215512, A215575, A215694, A215695, A108716, A215794, A215828, A215817, A215877, A094429, A217274.

I:=[49,245,1078]; [7] cat [n le 3 select I[n] else 7*Self(n-1) -14*Self(n-2) + 7*Self(n-3): n in [1..30]]; // G. C. Greubel, May 09 2018

Table[(MatrixPower[{{0, 1, 0}, {0, 0, 1}, {7, -14, 7}}, n].{1, 0, 0})[[3]], {n, 22}] (* Robert G. Wilson v, May 08 2004 *)
Join[{7}, LinearRecurrence[{7,-14,7}, {49,245,1078}, 50]] (* Roman Witula, Aug 13 2012 *)(* corrected by G. C. Greubel, May 09 2018 *)

x='x+O('x^30); Vec(7*x/(1-7*x+14*x^2-7*x^3)) \\ G. C. Greubel, May 09 2018

a(n) = (([0, 1, 0; 0, 0, 1; 7, -14, 7]^n)*[1,0,0]~)[3]; \\ Michel Marcus, May 10 2018

From Colin Barker, Jun 19 2012: (Start)
a(n) = 7*a(n-1)-14*a(n-2)+7*a(n-3).
G.f.: 7*x/(1-7*x+14*x^2-7*x^3). (End)
-a(n) = s(2)*s(1)^(2*n+3) + s(4)*s(2)^(2*n+3) + s(1)*s(4)^(2*n+3), where s(j) := 2*sin(2*Pi*j/7); for the proof see A215494 and the Witula-Slota paper. This formula and the respective recurrence also give a(0)=a(-1)=0. - Roman Witula, Aug 13 2012

More terms from Robert G. Wilson v, May 08 2004

Name edited by Michel Marcus, May 10 2018

A217274 a(n) = 7a(n-1) - 14a(n-2) + 7*a(n-3) with a(0)=0, a(1)=1, a(2)=7.

Original entry on oeis.org

0, 1, 7, 35, 154, 637, 2548, 9996, 38759, 149205, 571781, 2184910, 8333871, 31750824, 120875944, 459957169, 1749692735, 6654580387, 25306064602, 96226175941, 365880389868, 1391138718116, 5289228800247, 20109822277181, 76457523763621, 290689756066542

Offset: 0

Roman Witula, Sep 29 2012

This is the Berndt-type sequence number 18 for the argument 2*Pi/7 defined by the relation
a(n)*sqrt(7) = c(4)*s(1)^(2n+1) + c(2)*s(4)^(2n+1) + c(1)*s(2)^(2n+1) = (1/s(4))*s(1)^(2n+2) + (1/s(2))*s(4)^(2n+2) + (1/s(1))*s(2)^(2n+2), where c(j) := 2*cos(2*Pi*j/7) and s(j) := 2*sin(2*Pi*j/7) (for the sums of the respective even powers see A094429). For the proof of this formula see the Witula/Slota and Witula references.
The definitions of the other Berndt-type sequences for the argument 2*Pi/7 (with numbers from 1 to 17) are in the cross references.
We note that all numbers of the form a(n)*7^(-floor((n+1)/3)) = A217444(n) are integers.
It can be proved that Sum_{k=2..n}a(k) = 7*(a(n-1) - a(n-2)).

			Writing c(j) as cj and s(k) as sk,
we have 7*sqrt(7) = c4*s1^5 + c2*s4^5 + c1*s2^5
and c4*s1^13 + c2*s4^13 + c1*s2^13 = 4(c4*s1^11 + c2*s4^11 + c1*s2^11).
We note that a(9) = 87*a(3)*a(2)^2 and a(11) = 2*a(3)*a(5)*a(2)^2.

G. C. Greubel, Table of n, a(n) for n = 0..1000
R. Witula, Ramanujan type trigonometric formulas, Demonstratio Math., Vol. XLV, No. 4, 2012, pp. 789-796.
Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6.
Index entries for linear recurrences with constant coefficients, signature (7,-14,7).

Cf. A033304, A094429, A094430, A094648, A108716, A215007, A215008, A215143, A215493, A215494, A215510, A215512, A215575, A215694, A215695, A215794, A215817, A215828, A215877, A217444.

I:=[0,1,7]; [n le 3 select I[n] else 7*Self(n-1)-14*Self(n-2)+7*Self(n-3): n in [1..30]]; // Vincenzo Librandi, Jul 26 2015

LinearRecurrence[{7,-14,7}, {0,1,7}, 30]
CoefficientList[Series[x/(1 - 7*x + 14*x^2 - 7*x^3), {x,0,50}], x] (* G. C. Greubel, Apr 16 2017 *)

a[0]:0$
a[1]:1$
a[2]:7$
a[n]:=7*a[n-1] - 14*a[n-2] + 7*a[n-3];
makelist(a[n], n, 0, 25); /* Martin Ettl, Oct 11 2012 */

concat(0, Vec(x/(1-7*x+14*x^2-7*x^3) + O(x^40))) \\ Michel Marcus, Jul 25 2015

G.f.: x/(1-7*x+14*x^2-7*x^3).

A215817 a(n) is the rational part of A(n) = (6-sqrt(7))A(n-1) - (12-4sqrt(7))A(n-2) + (8-3sqrt(7))A(n-3) with A(0)=3, A(1)=6-sqrt(7), A(2)=19-4sqrt(7).

Original entry on oeis.org

3, 6, 19, 66, 237, 866, 3202, 11948, 44917, 169914, 646134, 2467988, 9462498, 36398004, 140399901, 542894726, 2103745125, 8167514346, 31762430143, 123704647562, 482435457922, 1883712663668, 7363103647479, 28809291337986, 112820819490970, 442175629583316

Offset: 0

Roman Witula, Aug 25 2012

nonn

The Berndt-type sequence number 14 for the argument 2Pi/7 defined by requiring a(n) to be the rational part of the trigonometric sum A(n) := c(1)^(2*n) + c(2)^(2*n) + c(4)^(2*n), where c(j) := 2*cos(Pi/4 + 2*Pi*j/7) = 2*cos((7+8*j)*Pi/28). We note that (A(n)-a(n))/sqrt(7) = A215877(n) are all integers. We have A(n)=2^n*O(n;i/2), where O(n;d) denote the big omega function with index n for the argument d in C defined in comments to A215794 (see also Witula-Slota's paper - Section 6). From the respective recurrence relation for this function we generate the title recurrence for A(n).

Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2*Pi/7, Journal of Integer Sequences, Vol. 12 (2009), Article 09.8.5

Cf. A215493, A215494, A215143, A215510, A094429, A215794.

a(n) = rational part of c(1)^(2n) + c(2)^(2n) + c(4)^(2n) = (1-s(1))^n + (1-s(2))^n + (1-s(4))^n, where c(j) := 2*cos((7+8*j)/28) and s(j) := sin(2*Pi*j/7).

Empirical g.f.: -(2*x-1)*(6*x^4 -40*x^3 +58*x^2 -24*x +3) / (x^6 -24*x^5 +86*x^4 -104*x^3 +53*x^2 -12*x +1). - Colin Barker, Jun 01 2013

cp's OEIS Frontend

A108716 a(n) = tan(Pi/14)^(-2n) + tan(3*Pi/14)^(-2n) + tan(5*Pi/14)^(-2n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A033304 Expansion of (2 + 2*x - 3*x^2) / (1 - 2*x - x^2 + x^3).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Formula

A094648 An accelerator sequence for Catalan's constant.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A215575 a(n) = 7*(a(n-1) - a(n-2) - a(n-3)), with a(0)=3, a(1)=7, a(2)=35.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A215494 a(n) = 7*a(n-1) - 14*a(n-2) + 7*a(n-3) with a(1)=7, a(2)=21, a(3)=70.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Formula

A215512 a(n) = 5*a(n-1) - 6*a(n-2) + a(n-3), with a(0)=1, a(1)=3, a(2)=8.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

A108716 a(n) = tan(Pi/14)^(-2n) + tan(3Pi/14)^(-2n) + tan(5Pi/14)^(-2n).

A033304 Expansion of (2 + 2x - 3x^2) / (1 - 2*x - x^2 + x^3).

A215494 a(n) = 7a(n-1) - 14a(n-2) + 7*a(n-3) with a(1)=7, a(2)=21, a(3)=70.

A215512 a(n) = 5a(n-1) - 6a(n-2) + a(n-3), with a(0)=1, a(1)=3, a(2)=8.

A217274 a(n) = 7a(n-1) - 14a(n-2) + 7*a(n-3) with a(0)=0, a(1)=1, a(2)=7.

A215817 a(n) is the rational part of A(n) = (6-sqrt(7))A(n-1) - (12-4sqrt(7))A(n-2) + (8-3sqrt(7))A(n-3) with A(0)=3, A(1)=6-sqrt(7), A(2)=19-4sqrt(7).