結果(2<N<51、プログラムは後半にあります)
************* test  Max X & 2*3^N **************************************  
   N     4    2                      Max X                       fx/2*3^N 
    2    1    0                                            1    0.00000
    3    1    1                                            5    0.60818
    4    2    1                                            3    0.27894
    5    2    2                                           37    0.71459
    6    3    2                                            7    0.32589
    7    4    2                                            5    0.22156
    8    4    3                                           18    0.35055
    9    5    3                                            9    0.24200
   10    5    4                                           62    0.39092
   11    6    4                                           19    0.25561
   12    7    4                                           13    0.20285
   13    7    5                                           42    0.26978
   14    8    5                                           24    0.21225
   15    8    6                                          108    0.29168
   16    9    6                                           44    0.22098
   17    9    7                                          638    0.35355
   18   10    7                                           87    0.23096
   19   11    7                                           53    0.19452
   20   11    8                                          197    0.24505
   21   12    8                                           95    0.20113
   22   12    9                                          648    0.27248
   23   13    9                                          179    0.20871
   24   14    9                                          116    0.18327
   25   14   10                                          371    0.21869
   26   15   10                                          202    0.18858
   27   15   11                                          964    0.23488
   28   16   11                                          367    0.19460
   29   16   12                                         7818    0.28499
   30   17   12                                          720    0.20212
   31   18   12                                          428    0.18011
   32   18   13                                         1640    0.21304
   33   19   13                                          759    0.18502
   34   19   14                                         5879    0.23491
   35   20   14                                         1428    0.19094
   36   21   14                                          911    0.17410
   37   21   15                                         2995    0.19891
   38   22   15                                         1583    0.17820
   39   22   16                                         8115    0.21210
   40   23   16                                         2886    0.18300
   41   23   17                                       129538    0.26372
   42   24   17                                         5711    0.18913
   43   25   17                                         3326    0.17316
   44   25   18                                        13357    0.19819
   45   26   18                                         5918    0.17714
   46   26   19                                        53958    0.21735
   47   27   19                                        11222    0.18204
   48   28   19                                         7040    0.16931
   49   28   20                                        23948    0.18874
   50   29   20                                        12277    0.17268

結果(499<N<551、プログラムは後半にあります)
************* test  Max X & 2*3^N *************************************  
  N     4    2             Max X                                fx/2*3^N 
  500  292  207       14026236000412862538919629802084044914    0.15583
  501  293  207        8890023569007029884400803102605923727    0.15469
  502  293  208       29586016787206543851874158988788414567    0.15656
  503  294  208       15457680686733768470502839737732255518    0.15507
  504  294  209       81729405809263481892220397576757744425    0.15778
  505  295  209       28244395839386400355286669877687781298    0.15555
  506  295  210     2305491456877524756060777854109321875445    0.16316
  507  296  210       56144878417751026267355496252976037050    0.15617
  508  297  210       32421162792915826576670875359394456435    0.15487
  509  297  211      132807182312122619116814549762499010084    0.15709
  510  298  211       57788590710972403342110364868481234320    0.15530
  511  298  212      567678013660465357129157285953123087020    0.15907
  512  299  212      109979333238429819663817767177803686976    0.15584
  513  300  212       68535826029484848544523944893730635777    0.15469
  514  300  213      236463578630192919214216901958610834151    0.15659
  515  301  213      119721770000534682027712096968342658667    0.15508
  516  301  214      692162315670388774331374170641657802881    0.15788
  517  302  214      220383885725405458384330545935785604933    0.15555
  518  303  214      145833636587771060357802977334945486835    0.15453
  519  303  215      444083788128952671262482466521751395037    0.15618
  520  304  215      250531055352578652618201670691846545407    0.15488
  521  304  216     1085398899255849360389855962667346584220    0.15715
  522  305  216      449218008987287108319984737939956110978    0.15530
  523  305  217     5577554232159194337707048590909470019398    0.15939
  524  306  217      863656434676446385209900365825660127180    0.15585
  525  307  217      528596779694672487224398927547688629953    0.15470
  526  307  218     1896570735757433864590400055496564427501    0.15662
  527  308  218      927886998342813580752927967540443049819    0.15508
  528  308  219     5951378595655953177961827984620572498106    0.15800
  529  309  219     1721568064548901632064021003592400385825    0.15556
  530  310  219     1122972297268171145288004329358255629580    0.15453
  531  310  220     3520699221131060036402445344143112964350    0.15620
  532  311  220     1937025065851562145357961800231751609550    0.15489
  533  311  221     8932615488707306292573094088402219263157    0.15721
  534  312  221     3495096733914577702264850769762745521078    0.15531
  535  312  222    60322400354239345283521723847614415484789    0.15987
  536  313  222     6793388066324953157993138553127867766743    0.15586
  537  314  222     4078798672664674077440400556399973609320    0.15471
  538  314  223    15272178071681201072619610566941232585980    0.15665
  539  315  223     7196587368016130181602668156124964669430    0.15509
  540  315  224    52191263256488099399560121808050607709110    0.15815
  541  316  224    13464849210867608826988430748294464552457    0.15557
  542  317  224     8650682598027389731040900352241308307075    0.15454
  543  317  225    27983891497510713550418778409205817349096    0.15623
  544  318  225    14985175582915341546679612981818915080733    0.15489
  545  318  226    74125478430770618003436639315223777634948    0.15728
  546  319  226    27219054821227439643130379297509766678917    0.15532
  547  319  227   803654899675972022635316839289241513202233    0.16067
  548  320  227    53531576920845925767848209763165161625659    0.15588
  549  321  227    31488377829549806890847920902414515295274    0.15471
  550  321  228   123538436033753313911981570577554487952080    0.15670

プログラム(UBASIC)
   10   'asave "REVERS2"
   20   print "************* test  Max X & 2*3^N ***********  "
   30   print "        N     4    2      Max X                      fx/2*3^N "
   40   for N=2 to 50
   45   '4の指数の最大値L4とそのときの2の指数L2
   50   L4=int(N*(log(3)/log(2)-1)):L2=N-1-L4
   55   print using(5,0),N; print using(5,0),L4;L2;
   57   R=0:L3=0 'R:定数項 L3:3の指数
   58   'Rの最大値のため2を3に変更してゆく
   60   for L=L2 to 0 step -1
   80   R=R+(4^L4*2^L*3^L3)/3^N
   90   inc L3
  100   next L:L=0
  105   '4を3に変更する
  110   for M=L4-1 to 0 step -1
  130   R=R+(2^L*4^M*3^L3)/3^N
  140   inc L3
  150   next M
  160   Q=(4*4^L4*2^L2)/3^N
  170   FX=R/(Q-1)
  210   print using(40,0),FX;:print using(5,5),log(FX)/((N-1)*log(3)+log(2))
  220   next


戻る