The OEIS is supported by the many generous donors to the OEIS Foundation.
%I #8 Apr 07 2019 16:59:42
%S 1,1,7,65,859,14733,311111,7807241,227032211,7507677445,278238129887,
%T 11422778680257,514581516871339,25237339884056125,1338608176888717303,
%U 76351140762424998521,4660096802861702344387,303058994467535502382197,20920390811865673069060367,1527794647255129205046171185
%N a(0)=a(1)=1; thereafter a(n) = (4*n-3)*a(n-1) + 2*a(n-2).
%H D. H. Lehmer, <a href="http://www.jstor.org/stable/1968107">Arithmetical periodicities of Bessel functions</a>, Annals of Mathematics, 33 (1932): 143-150. The sequence is on page 148.
%p f:=proc(n) option remember; if n <= 1 then 1 else 2*f(n-2)+(4*n-3)*f(n-1); fi; end;
%p [seq(f(n),n=0..30)];
%t RecurrenceTable[{a[0]==a[1]==1,a[n]==(4n-3)a[n-1]+2a[n-2]},a,{n,20}] (* _Harvey P. Dale_, Apr 07 2019 *)
%K nonn
%O 0,3
%A _N. J. A. Sloane_, May 09 2016