bash-2.05$ gp
Reading GPRC: /home/his/.gprc ...Done.

                  GP/PARI CALCULATOR Version 2.1.1 (released)
                i386 running netbsd (ix86 kernel) 32-bit version
                (readline v1.0 enabled, extended help available)

                       Copyright (C) 2000 The PARI Group

PARI/GP is free software, covered by the GNU General Public License, and 
comes WITHOUT ANY WARRANTY WHATSOEVER.

Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.

   realprecision = 28 significant digits
   seriesprecision = 16 significant terms
   format = g0.28

parisize = 4000000, primelimit = 500000
(21:47) gp > read("bernoulli.gp")
(21:47) gp > b[1]
%1 = 1/6
(21:47) gp > b[2]
%2 = -1/30
(21:47) gp > b[3]
%3 = 1/42
(21:48) gp > b[4]
%4 = -1/30
(21:48) gp > b[5]
%5 = 5/66
(21:48) gp > b[6]
%6 = -691/2730
(21:48) gp > b[7]
%7 = 7/6
(21:48) gp > b[8]
%8 = -3617/510
(21:48) gp > b[9]
%9 = 43867/798
(21:48) gp > b[10]
%10 = -174611/330
(21:48) gp > b[11]
%11 = 854513/138
(21:48) gp > b[12]
%12 = -236364091/2730
(21:48) gp > b[13]
%13 = 8553103/6
(21:48) gp > b[14]
%14 = -23749461029/870
(21:48) gp > b[15]
%15 = 8615841276005/14322
(21:48) gp > b[16]
%16 = -7709321041217/510
(21:48) gp > b[17]
%17 = 2577687858367/6
(21:48) gp > b[18]
%18 = -26315271553053477373/1919190
(21:48) gp > b[19]
%19 = 2929993913841559/6
(21:49) gp > b[20]
%20 = -261082718496449122051/13530
(21:49) gp > b[100]
%21 = -498384049428333414764928632140399662108495887457206674968055822617263669621523687568865802302210999132601412697613279391058654527145340515840099290478026350382802884371712359337984274122861159800280019110197888555893671151/1366530
(21:49) gp > z[1]
%22 = 1/6
(21:49) gp > factor(%22)
%23 = 
[2 -1]

[3 -1]

(21:49) gp > z[2]
%24 = 1/90
(21:49) gp > factor(%24)
%25 = 
[2 -1]

[3 -2]

[5 -1]

(21:49) gp > z[3]
%26 = 1/945
(21:49) gp > factor(%26)
%27 = 
[3 -3]

[5 -1]

[7 -1]

(21:49) gp > z[4]
%28 = 1/9450
(21:49) gp > factor(%28)
%29 = 
[2 -1]

[3 -3]

[5 -2]

[7 -1]

(21:50) gp > z[5]
%30 = 1/93555
(21:50) gp > factor(%30)
%31 = 
[3 -5]

[5 -1]

[7 -1]

[11 -1]

(21:50) gp > z[6]
%32 = 691/638512875
(21:50) gp > factor(%32)
%33 = 
[3 -6]

[5 -3]

[7 -2]

[11 -1]

[13 -1]

[691 1]

(21:50) gp > z[7]
%34 = 2/18243225
(21:50) gp > factor(%34)
%35 = 
[2 1]

[3 -6]

[5 -2]

[7 -1]

[11 -1]

[13 -1]

(21:50) gp > z[8]
%36 = 3617/325641566250
(21:50) gp > factor(%36)
%37 = 
[2 -1]

[3 -7]

[5 -4]

[7 -2]

[11 -1]

[13 -1]

[17 -1]

[3617 1]

(21:51) gp > z[9]
%38 = 43867/38979295480125
(21:51) gp > factor(%38)
%39 = 
[3 -9]

[5 -3]

[7 -3]

[11 -1]

[13 -1]

[17 -1]

[19 -1]

[43867 1]

(21:51) gp > z[10]
%40 = 174611/1531329465290625
(21:51) gp > factor(%40)
%41 = 
[3 -9]

[5 -5]

[7 -2]

[11 -2]

[13 -1]

[17 -1]

[19 -1]

[283 1]

[617 1]

(21:51) gp > z[11]
%42 = 155366/13447856940643125
(21:51) gp > factor(%42)
%43 = 
[2 1]

[3 -10]

[5 -4]

[7 -3]

[11 -1]

[13 -1]

[17 -1]

[19 -1]

[23 -1]

[131 1]

[593 1]

(21:51) gp > z[12]
%44 = 236364091/201919571963756521875
(21:51) gp > factor(%44)
%45 = 
[3 -11]

[5 -5]

[7 -4]

[11 -2]

[13 -2]

[17 -1]

[19 -1]

[23 -1]

[103 1]

[2294797 1]

(21:52) gp > z[15]
%46 = 6892673020804/5660878804669082674070015625
(21:52) gp > factor(%46)
%47 = 
[2 2]

[3 -15]

[5 -6]

[7 -5]

[11 -3]

[13 -2]

[17 -1]

[19 -1]

[23 -1]

[29 -1]

[31 -1]

[1721 1]

[1001259881 1]

(21:52) gp > z[100]
%48 = 996768098856666829529857264280799324216991774914413349936111645234527339243047375137731604604421998265202825395226558782117309054290681031680198580956052700765605768743424718675968548245722319600560038220395777111787342302/2682678748792657844957504192313280657551803049278355275671666881580642758576467817615493645217977237214155689404787155170845497733836863647685885197919191727452679238952541411298115541287013688972773507748859386210346035176197388875022427877722880764252312145081723341902733317236524547144682628641021437942981719970703125
(21:52) gp > quit
Good bye!
bash-2.05$