OFFSET
0,3
COMMENTS
Also a(n) = -sin(2n*arccos(sqrt(n)))^2 = -sin(2n*arcsin(sqrt(n)))^2.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..170
MAPLE
A173150 := proc(n) sinh(2*n*arccosh(sqrt(n))) ; %^2 ; expand(%) ; simplify(%) ; end proc: # R. J. Mathar, Feb 26 2011
MATHEMATICA
Table[Round[-Sin[2 n ArcCos[Sqrt[n]]]^2], {n, 0, 20}] (* Artur Jasinski, Feb 11 2010 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Feb 11 2010
EXTENSIONS
Minor edits by Vaclav Kotesovec, Apr 05 2016
STATUS
approved