Return to program CVS log

Up to [local] / OpenXM_contrib / pari / src / test / 32

Annotation of OpenXM_contrib/pari/src/test/32/program, Revision 1.1

1.1     ! maekawa     1:    realprecision = 38 significant digits
        !             2:    echo = 1 (on)
        !             3: ? alias(ln,log)
        !             4: ? ln(2)
        !             5: 0.69314718055994530941723212145817656807
        !             6: ? for(x=1,5,print(x!))
        !             7: 1
        !             8: 2
        !             9: 6
        !            10: 24
        !            11: 120
        !            12: ? fordiv(10,x,print(x))
        !            13: 1
        !            14: 2
        !            15: 5
        !            16: 10
        !            17: ? forprime(p=1,30,print(p))
        !            18: 2
        !            19: 3
        !            20: 5
        !            21: 7
        !            22: 11
        !            23: 13
        !            24: 17
        !            25: 19
        !            26: 23
        !            27: 29
        !            28: ? forstep(x=0,Pi,Pi/12,print(sin(x)))
        !            29: 0.E-38
        !            30: 0.25881904510252076234889883762404832834
        !            31: 0.49999999999999999999999999999999999999
        !            32: 0.70710678118654752440084436210484903928
        !            33: 0.86602540378443864676372317075293618346
        !            34: 0.96592582628906828674974319972889736763
        !            35: 1.0000000000000000000000000000000000000
        !            36: 0.96592582628906828674974319972889736764
        !            37: 0.86602540378443864676372317075293618348
        !            38: 0.70710678118654752440084436210484903930
        !            39: 0.50000000000000000000000000000000000002
        !            40: 0.25881904510252076234889883762404832838
        !            41: 4.7019774032891500318749461488889827112 E-38
        !            42: ? forvec(x=[[1,3],[-2,2]],print1([x[1],x[2]]," "));print(" ");
        !            43: [1, -2] [1, -1] [1, 0] [1, 1] [1, 2] [2, -2] [2, -1] [2, 0] [2, 1] [2, 2] [3
        !            44: , -2] [3, -1] [3, 0] [3, 1] [3, 2]
        !            45: ? getheap
        !            46: [3, 29]
        !            47: ? getrand
        !            48: 1
        !            49: ? getstack
        !            50: 0
        !            51: ? if(3<2,print("bof"),print("ok"));
        !            52: ok
        !            53: ? kill(y);print(x+y);reorder([x,y]);print(x+y);
        !            54: x + y
        !            55: x + y
        !            56: ? f(u)=u+1;
        !            57: ? print(f(5));kill(f);
        !            58: 6
        !            59: ? f=12
        !            60: 12
        !            61: ? g(u)=if(u,,return(17));u+2
        !            62: ? g(2)
        !            63: 4
        !            64: ? g(0)
        !            65: 17
        !            66: ? setrand(10)
        !            67: 10
        !            68: ? n=33;until(n==1,print1(n," ");if(n%2,n=3*n+1,n=n/2));print(1)
        !            69: 33 100 50 25 76 38 19 58 29 88 44 22 11 34 17 52 26 13 40 20 10 5 16 8 4 2 1
        !            70: ? m=5;while(m<20,print1(m," ");m=m+1);print()
        !            71: 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
        !            72: ? default(realprecision,28)
        !            73:    realprecision = 28 significant digits
        !            74: ? default(seriesprecision,12)
        !            75:    seriesprecision = 12 significant terms
        !            76: ? print((x-12*y)/(y+13*x));
        !            77: (x - 12*y)/(13*x + y)
        !            78: ? print([1,2;3,4])
        !            79: [1, 2; 3, 4]
        !            80: ? print1(x+y);print(x+y);
        !            81: x + yx + y
        !            82: ? print((x-12*y)/(y+13*x));
        !            83: (x - 12*y)/(13*x + y)
        !            84: ? print([1,2;3,4])
        !            85: [1, 2; 3, 4]
        !            86: ? print1(x+y);print1(" equals ");print(x+y);
        !            87: x + y equals x + y
        !            88: ? print1("give a value for s? ");s=input();print(1/(s^2+1))
        !            89: give a value for s? printtex((x+y)^3/(x-y)^2)
        !            90: {{x^{3} + {{3}y}x^{2} + {{3}y^{2}}x + {y^{3}}}\over{x^{2} - {{2}y}x + {y^{2}
        !            91: }}}
        !            92: 1
        !            93: ? for(i=1,100,for(j=1,25,if(i+j==32,break(2)));print(i))
        !            94: 1
        !            95: 2
        !            96: 3
        !            97: 4
        !            98: 5
        !            99: 6
        !           100: ? u=v=p=q=1;for(k=1,400,w=u+v;u=v;v=w;p*=w;q=lcm(q,w);if(k%50==0,print(k" "l
        !           101: og(p)/log(q))));
        !           102: 50 1.561229126903099279206171725
        !           103: 100 1.601335375590875348711141031
        !           104: 150 1.606915548673659127523394774
        !           105: 200 1.618659998991528481508175175
        !           106: 250 1.626284706204746765086080988
        !           107: 300 1.627822776845103001192024532
        !           108: 350 1.632105905172986668189652273
        !           109: 400 1.632424285532931448171405619
        !           110: ? install(addii,GG)
        !           111: ? addii(1,2)
        !           112: 3
        !           113: ? kill(addii)
        !           114: ? getheap
        !           115: [21, 2992]
        !           116: ? print("Total time spent: ",gettime);
        !           117: Total time spent: 1151
        !           118: ? \q