OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..1306
FORMULA
G.f.: (1+3*x-sqrt(1-6*x+x^2))/2. - Michael Somos, Jun 08 2000
D-finite with recurrence n*a(n) +3*(-2*n+3)*a(n-1) +(n-3)*a(n-2)=0. - R. J. Mathar, Feb 25 2015
MAPLE
f:= gfun:-rectoproc({n*a(n) +3*(-2*n+3)*a(n-1) +(n-3)*a(n-2)=0, a(1)=3, a(2)=2}, a(n), remember):
map(f, [$1..30]); # Robert Israel, Sep 28 2025
MATHEMATICA
Rest[CoefficientList[Series[(1+3x-Sqrt[1-6x+x^2])/2, {x, 0, 40}], x]] (* Harvey P. Dale, Aug 26 2013 *)
PROG
(PARI) a(n)=polcoeff((1+3*x-sqrt(1-6*x+x^2+x*O(x^n)))/2, n)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved
