Return to program CVS log

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

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

1.1     ! maekawa     1: \p 38
        !             2: \e
        !             3: alias(ln,log)
        !             4: ln(2)
        !             5: for(x=1,5,print(x!))
        !             6: fordiv(10,x,print(x))
        !             7: forprime(p=1,30,print(p))
        !             8: forstep(x=0,Pi,Pi/12,print(sin(x)))
        !             9: forvec(x=[[1,3],[-2,2]],print1([x[1],x[2]]," "));print(" ");
        !            10: getheap
        !            11: getrand
        !            12: getstack
        !            13: if(3<2,print("bof"),print("ok"));
        !            14: kill(y);print(x+y);reorder([x,y]);print(x+y);
        !            15: f(u)=u+1;
        !            16: print(f(5));kill(f);
        !            17: f=12
        !            18: g(u)=if(u,,return(17));u+2
        !            19: g(2)
        !            20: g(0)
        !            21: setrand(10)
        !            22: n=33;until(n==1,print1(n," ");if(n%2,n=3*n+1,n=n/2));print(1)
        !            23: m=5;while(m<20,print1(m," ");m=m+1);print()
        !            24: \\
        !            25: default(realprecision,28)
        !            26: default(seriesprecision,12)
        !            27: print((x-12*y)/(y+13*x));
        !            28: print([1,2;3,4])
        !            29: print1(x+y);print(x+y);
        !            30: print((x-12*y)/(y+13*x));
        !            31: print([1,2;3,4])
        !            32: print1(x+y);print1(" equals ");print(x+y);
        !            33: print1("give a value for s? ");s=input();print(1/(s^2+1))
        !            34: printtex((x+y)^3/(x-y)^2)
        !            35: for(i=1,100,for(j=1,25,if(i+j==32,break(2)));print(i))
        !            36: {
        !            37:   u=v=p=q=1;
        !            38:   for (k=1,400,
        !            39:     w=u+v; u=v; v=w;
        !            40:     p *= w; q=lcm(q,w);
        !            41:     if (k%50==0,
        !            42:       print(k" "log(p)/log(q))
        !            43:     )
        !            44:   );
        !            45: }
        !            46: /* added in 1.924.13 */
        !            47: install(addii,GG)
        !            48: addii(1,2)
        !            49: kill(addii)
        !            50: /* up to here */
        !            51: getheap
        !            52: print("Total time spent: ",gettime);
        !            53: \q