09/04/05 20:31:00
>>823
(略解)
ベクトル (a_1, a_2, a_3, a_4, a_5)† に行列
[ 1/2, 1/3, 1/4, 1/5, 1/6 ]
[ 1/5, 1/6, 1/7. 1/8, 1/9 ]
[ 1/10, 1/11, 1/12, 1/13, 1/14 ]
[ 1/17, 1/18, 1/19, 1/20, 1/21 ]
[ 1/26, 1/27, 1/28, 1/29, 1/30 ]
を作用した結果が
(1, 1/4, 1/9, 1/16, 1/25)†
上記の逆行列は
[ (25*13*17)/(8*9), -(4*5*13*17)/3, (5*11*169*17)/16, -(2*25*13*17*19)/9, (125*13*17*29)/(16*3)]
[ -(81*11)/2, (4*243*11), -(2187*11*13)/4, (2*27*11*17*19), -(27*5*11*13*29)/4 ]
[ (49*19), -(64*3*7*19), (3*49*11*13*19)/2, -(32*49*17*19), (3*5*7*13*19*29)/2 ]
[ -(4*5*13*29)/9, (8*25*13*29)/3,-(2*25*11*13*29), (8*5*13*17*19*29)/9,-(2*125*169*29)/3 ]
[ (9*5*49)/8, -(4*9*5*49), (243*5*7*11*13)/16,-(2*9*5*7*17*19), (27*5*49*13*29)/16 ]
これを使ってa_kを求めると、
a_1 = 1105/72,
a_2 = -2673/40,
a_3 = 1862/15,
a_4 = -1885/18,
a_5 = 1323/40,
これより
b_6 = 187465/(3*37*38*39*41) ≒ 1.00061649483987･･･ / 36,
b_7 = 1197/(5*13*17*53) ≒ 1.00150260394436･･･ / 49,
b_8 = 85345/(16*13*17*23*67) ≒ 1.00240485551780･･･ / 64,
b_9 = 277289/(9*17*41*43*83) ≒ 1.00321917612728･･･ / 81,
b_10=12117378/(3*25*7*13*17*101*103) ≒ 1.00391855290609･･･ / 100,
b_0 = 13489 / 3600 ≒ 3.74694444444444･･･
ここに b_k = a_1/(1+k^2) + a_2/(2+k^2) + a_3/(3+k^2) + a_4/84+k^2) + a_5/(5+k^2),