2024年6月22日土曜日

240622

PARI


Continued fraction for Pi^n(2)

Pi^3の連分数展開も興味深い。

for(n=1, 10, print([n, Pi^n, contfrac(Pi^n)]))

出力結果
[1, 3.1415926535897932384626433832795028842, [3, 7, 15, 1, 292, 1, 1, 1, 2, 1, 3, 1, 14, 2, 1, 1, 2, 2, 2, 2, 1, 84, 2, 1, 1, 15, 3, 13, 1, 4, 2, 6, 6]]
[2, 9.8696044010893586188344909998761511353, [9, 1, 6, 1, 2, 47, 1, 8, 1, 1, 2, 2, 1, 1, 8, 3, 1, 10, 5, 1, 3, 1, 2, 1, 1, 3, 15, 1, 1, 2, 2, 1, 3, 2, 7, 1, 9, 18]]
[3, 31.006276680299820175476315067101395202, [31, 159, 3, 7, 1, 13, 2, 1, 3, 1, 12, 2, 2, 4, 34, 2, 43, 3, 1, 3, 2, 1, 1, 5, 1, 1, 4, 1, 5, 4, 2, 4]]
[4, 97.409091034002437236440332688705111250, [97, 2, 2, 3, 1, 16539, 1, 6, 7, 6, 8, 6, 3, 9, 1, 1, 1, 18, 1, 4, 1, 13, 1, 2, 1, 127, 2]]
[5, 306.01968478528145326274131004343560648, [306, 50, 1, 4, 60, 1, 4, 23, 1, 2, 1, 1, 1, 17, 1, 5, 1, 1, 1, 2, 4, 10, 4, 1, 3, 2, 1, 36, 7, 1, 27]]
[6, 961.38919357530443703021944365241989886, [961, 2, 1, 1, 3, 9, 1, 5, 1, 3, 6, 1, 3, 1, 77, 11, 3, 1, 11, 14, 4, 4, 27, 5, 178, 6]]
[7, 3020.2932277767920675142064930720418319, [3020, 3, 2, 2, 3, 2, 10, 4, 1, 1, 2, 1, 3, 13, 58, 2, 50, 1, 2, 21, 1, 1, 2, 1, 14, 1, 182, 2, 1, 1, 2]]
[8, 9488.5310160705740071285755039067657967, [9488, 1, 1, 7, 1, 1, 3, 1, 1, 1, 4, 10, 13, 1, 1, 2, 3, 6, 1, 1, 1, 135, 1, 2, 2, 3, 1, 2, 1, 7, 3, 1, 1, 4, 2, 7, 1, 2]]
[9, 29809.099333446211666509402401239655364, [29809, 10, 14, 1, 9, 3, 1, 5, 8, 16, 3, 15, 7, 39, 3, 3, 1, 1, 2, 1, 3, 4, 1, 1, 7, 2, 1, 2, 1, 3]]
[10, 93648.047476083020973716690184919345635, [93648, 21, 15, 1, 4, 2, 1, 3, 2, 1, 4, 2, 74, 26, 5, 3, 1, 5, 17, 1, 2, 1, 1, 2, 1, 14, 2, 2, 8]]