OFFSET
0,3
COMMENTS
The triples (X(n) = A389120(n) * phi^(3/2) - a(n) * phi^(1/2), Y(n) = A138573(n) * phi - A138574(n), 1) where phi = (1+sqrt(5))/2 form right triangles with hypotenuse 1. X(n)^2 + Y(n)^2 = 1 for all n > 0.
Alternatively the sequences are related with this integer matrix
+- -+
| 1, -1, 1, -1 |
| -1, 0, -1, 0 |
M= | 1, 0, 1, -1 |
| 0, 1, -1, 0 |
+- -+
+- -+ +- -+
|- a(n+1) | |- a(n) |
+- -+ +- -+
LINKS
Matthew L. LaSelle, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (2,2,2,-1).
FORMULA
G.f.: x*(1 + 2*x + x^2) / (1 - 2*x - 2*x^2 - 2*x^3 + x^4).
EXAMPLE
When n = 1: (Kepler triangle)
X(1) = 1 * phi^(3/2) - 1 * phi^(1/2) = 0.78615137775742328606955858584295892953,
Y(1) = 1 * phi - 1 = 0.61803398874989484820458683436563811772.
When n = 2:
X(2) = 2 * phi^(3/2) - 4 * phi^(1/2) = -0.97173654351329135636572775178906512438,
Y(2) = 2 * phi - 3 = 0.23606797749978969640917366873127623544.
When n = 3:
X(3) = 7 * phi^(3/2) - 11 * phi^(1/2) = 0.41498104624568714547722025395074653983,
Y(3) = 5 * phi - 9 = -0.90983005625052575897706582817180941141.
MATHEMATICA
LinearRecurrence[{2, 2, 2, -1}, {0, 1, 4, 11}, 30] (* Paolo Xausa, Jan 03 2026 *)
PROG
(PARI) my(x='x+O('x^50)); concat([0], Vec( x*(1 + 2*x + x^2) / (1 - 2*x - 2*x^2 - 2*x^3 + x^4) ))
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Chittaranjan Pardeshi, Sep 23 2025
STATUS
approved
