|
|
A073889
|
|
a(1) = 1, a(2) = 3; for n>2, a(n) = {a(n-1)*a(n+1)}^(1/n) or a(n+1) = a(n)^n/a(n-1).
|
|
3
|
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
MAPLE
|
if n<3 then 2*n-1;
else a(n-1)^(n-1)/a(n-2);
fi; end;
|
|
MATHEMATICA
|
a[n_]:= a[n]= If[n<3, 2*n-1, a[n-1]^(n-1)/a[n-2]]; Table[a[n], {n, 7}] (* G. C. Greubel, May 18 2020 *)
|
|
PROG
|
(Magma) I:=[1, 3]; [n le 2 select I[n] else Self(n-1)^(n-1)/Self(n-2): n in [1..7]]; # G. C. Greubel, May 18 2020
(Sage)
def a(n):
if (n<3): return 2*n-1
else: return a(n-1)^(n-1)/a(n-2)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
The next term is too large to include.
|
|
STATUS
|
approved
|
|
|
|