Lehmer's Second Interesting Series(1)
S0(2) = 1 * π/2 + 1,
S1(2) = 2 * π/2 + 3,
S2(2) = 7 * π/2 + 11,
S3(2) = 35 * π/2 + 55,
...
となるが、
数列 1, 2, 7, 35, ...
と
数列 1, 3, 11, 55, ...
を出力してみた。
(01:26) gp > N=20; x='x+O('x^N); f=sqrt(exp(x)/(2-exp(x))); Vec(serlaplace(f/(2-exp(x))))
%1 = [1, 2, 7, 35, 226, 1787, 16717, 180560, 2211181, 30273047, 458186752, 7596317885, 136907048461, 2665084902482, 55726440112987, 1245661569161135, 29642264728189066, 748158516941653967, 19962900431638852297, 561472467839585937560]
(01:26) gp > Vec(serlaplace(deriv(f*intformal(f))))
%2 = [1, 3, 11, 55, 355, 2807, 26259, 283623, 3473315, 47552791, 719718067, 11932268231, 215053088835, 4186305575415, 87534887434835, 1956680617267879, 46561960552921315, 1175204650272267479, 31357650670190565363, 881958890078887314567]
(01:26) gp >
0 件のコメント:
コメントを投稿
注: コメントを投稿できるのは、このブログのメンバーだけです。