This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A070190 #32 Feb 06 2024 08:07:15 %S A070190 1,0,10,0,270,240,10900,25200,551950,2116800,32458860,169092000, %T A070190 2120787900,13427013600,149506414200,1075081207200,11143223412750, %U A070190 87198375264000,865743970019500,7171730187336000,69416724049550020 %N A070190 Expansion of e.g.f. I_0(2*x)^5 + 2*Sum_{k>=1} I_k(2*x)^5, where I_n(z) are modified Bessel functions of order n. %C A070190 A modification of e.g.f. of A002898, where the exponent of I, which is 3, is here replaced by 5. %C A070190 U_10(n), in Labelle-Lacasse paper, number of closed paths of length n whose steps are 10th roots of unity. %H A070190 Andrew Howroyd, <a href="/A070190/b070190.txt">Table of n, a(n) for n = 0..200</a> %H A070190 Gilbert Labelle and Annie Lacasse, <a href="https://doi.org/10.46298/dmtcs.2937">Closed paths whose steps are roots of unity</a>, in FPSAC 2011, Reykjavik, Iceland DMTCS proc. AO, 2011, 599-610. %F A070190 a(n) ~ 5^(3/2) * 10^n / (4*Pi^2*n^2). - _Vaclav Kotesovec_, Jun 08 2021 %t A070190 With[{nmax = 25}, CoefficientList[Series[BesselI[0, 2*x]^5 + 2*Sum[BesselI[k, 2*x]^5, {k, 1, 2*nmax}], {x, 0, nmax}], x]*Range[0, nmax]!] (* _G. C. Greubel_, Nov 05 2018 *) %o A070190 (PARI) seq(n)={Vec(serlaplace(sum(k=0, n, if(k,2,1)*(x^k*besseli(k, 2*x + O(x^(n-k+1)))/k!)^5)))} \\ _Andrew Howroyd_, Nov 01 2018 %Y A070190 Cf. A002898. %K A070190 nonn %O A070190 0,3 %A A070190 _Karol A. Penson_, Apr 26 2002