>>824
正解です。素晴らしい。
ちなみに用意の解答
―-
f(n,x) = (-1/x d/dx)^n (exp x/x)
とおけば
x^2 f(n,x) = (2n-1)f(n-1,x) + f(n-2,x)。
とくに p[n] = f(n,1)、q[n] = f(n,-1)とおけば
p[n] = (2n-1)p[n-1] + p[n-2]、q[n] = (2n-1)q[n-1] + q[n-2]。
これとp[1] = 0、p[2] = e、q[1] = -2/e、q[2] = -7/eにより
c[n] = p[n]/e、d[n] = (-7p[n]/e + 2e q[n])2。
一方で (-1/x d/dx)^n (exp x/x)をマクローリン展開して lim[n→∞] f(n,±1)/(2n)!! = ±1。
以上により
lim[n→∞] c[n]/(2n)!! = 1/e、lim[n→∞] d[n]/(2n)!! = (-7/e+e)/2。

前わかスレに出てた変形ベッセル関数による表示を利用しています。
(本来のベッセル関数だとx=-1を代入できないのでちょっと一工夫してますが。)