%I #9 Feb 16 2025 08:33:46
%S 0,0,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,
%T 3,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,5,5,5,5,5,5,5,5,5,5,5,5,
%U 5,5,5,5,5,5,5,5,5,5,5,5,5,5
%N a(n) is the number of zeros of the Laguerre L(n, x) polynomial in the open interval (-1, +1).
%C The Laguerre polynomials are given by the sum: L(n,x) = Sum_{k=0..n} ((-1)^k)/k!*binomial(n,k)*x^k.
%C The first few Laguerre polynomials are:
%C L(0,x) = 1,
%C L(1,x) = -x + 1,
%C L(2,x) = ( x^2 - 4*x + 2)/2,
%C L(3,x) = (-x^3 + 9*x^2 - 18*x + 6)/6,
%C L(4,x) = ( x^4 - 16*x^3 + 72*x^2 - 96*x + 24)/24,
%C L(5,x) = (-x^5 + 25*x^4 - 200*x^3 + 600*x^2 - 600*x + 120)/120.
%C The number of occurrences a(n) = 0, 1, 2,.. is given by the sequence {2, 6, 11, 16, 21, ...}.
%H Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/LaguerrePolynomial.html">The World of Mathematics: Laguerre Polynomial</a>
%e a(3) = 1 because the zeros of L(3,x) = (-x^3 + 9*x^2 - 18*x + 6)/6 are x1=.4157745568..., x2=2.294280360... and x3=6.289945083... with the root x1 in the open interval (-1, +1). Hence, a(3) = 1.
%p for n from 0 to 90 do:it:=0:
%p y:=[fsolve(expand(LaguerreL(n,x)),x,real)]:for m from 1 to nops(y) do:if abs(y[m])<1 then it:=it+1:else fi:od: printf(`%d, `, it):od:
%Y Cf. A066667.
%K nonn,changed
%O 0,9
%A _Michel Lagneau_, May 22 2017