Michael Somos on Mon, 15 Nov 1999 19:42:39 -0500 (EST) |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
more zeta() faults |
What is think is interesting is the variety of different faults encountered : Reading GPRC: /home1/faculty/somos/.gprc ...Done. GP/PARI CALCULATOR Version 2.0.18 (development) i686 running linux (ix86 kernel) 32-bit version (readline v2.2 enabled, extended help available) .. parisize = 4000000, primelimit = 500000 gp> for(n=1,30,print(n" "zeta(I+10^n))) 1 1.000759077829856768809901195 - 0.0006401280267855239745764942372*I 2 1.000000000000000000000000000 - 2.623314226356272336124975847 E-31*I 3 1.000000000000000000000000000 - 2.993549263496598102783842505 E-302*I 4 1.000000000000000000000000000 - 1.601996728198064104058581242 E-3011*I 5 1.000000000000000000000000000 - 3.198125608525491241840590787 E-30104*I 6 1.000000000000000000000000000 - 3.226876165408552144077256540 E-301031*I *** overflow in R*R gp> for(n=1,30,print(n" "zeta(0+10^n))) 1 1.000994575127818085337145958 2 1.000000000000000000000000000 3 1.000000000000000000000000000 4 1.000000000000000000000000000 5 1.000000000000000000000000000 6 1.000000000000000000000000000 7 1.000000000000000000000000000 8 1.000000000000000000000000000 9 1.000000000000000000000000000 *** impossible assignment I-->S gp> for(n=1,30,print(n" "zeta(1/3*10^n))) 1 1.147356236882742872645615288 2 1.000000000092399026901528289 3 1.000000000000000000000000000 4 1.000000000000000000000000000 5 1.000000000000000000000000000 6 1.000000000000000000000000000 7 1.000000000000000000000000000 *** exponent too large in exp. gp> for(n=1,10,print(n" "zeta(1/10000^n))) 1 -0.5000919038861036320095103010 2 -0.5000000091893854323645512140 3 -0.5000000000009189385332056772 4 -0.5000000000000000918938718545 5 -0.5000000000000000000748188871 6 -0.5000000000000005856439856969 7 -0.5000000000070312676220292861 8 -0.5000001009959118357014707048 9 -0.5013897738467078603232921279 10 -18.70857365099656195236164656 Sorry to bear bad news. Shalom, Michael