Return to program CVS log

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

Imported pari-2.2.1(alpha).

   realprecision = 38 significant digits
   echo = 1 (on)
? alias(ln,log)
? ln(2)
0.69314718055994530941723212145817656807
? for(x=1,5,print(x!))
1
2
6
24
120
? fordiv(10,x,print(x))
1
2
5
10
? forprime(p=1,30,print(p))
2
3
5
7
11
13
17
19
23
29
? forstep(x=0,Pi,Pi/12,print(sin(x)))
0.E-38
0.25881904510252076234889883762404832834
0.49999999999999999999999999999999999999
0.70710678118654752440084436210484903928
0.86602540378443864676372317075293618346
0.96592582628906828674974319972889736763
1.0000000000000000000000000000000000000
0.96592582628906828674974319972889736764
0.86602540378443864676372317075293618348
0.70710678118654752440084436210484903930
0.50000000000000000000000000000000000002
0.25881904510252076234889883762404832838
4.7019774032891500318749461488889827112 E-38
? forvec(x=[[1,3],[-2,2]],print1([x[1],x[2]]," "));print(" ");
[1, -2] [1, -1] [1, 0] [1, 1] [1, 2] [2, -2] [2, -1] [2, 0] [2, 1] [2, 2] [3
, -2] [3, -1] [3, 0] [3, 1] [3, 2]  
? getheap
[3, 28]
? getrand
1
? getstack
0
? if(3<2,print("bof"),print("ok"));
ok
? kill(y);print(x+y);reorder([x,y]);print(x+y);
x + y
x + y
? f(u)=u+1;
? print(f(5));kill(f);
6
? f=12
12
? g(u)=if(u,,return(17));u+2
? g(2)
4
? g(0)
17
? setrand(10)
10
? n=33;until(n==1,print1(n," ");if(n%2,n=3*n+1,n=n/2));print(1)
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
? m=5;while(m<20,print1(m," ");m=m+1);print()
5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 
? default(realprecision,28)
   realprecision = 28 significant digits
? default(seriesprecision,12)
   seriesprecision = 12 significant terms
? print((x-12*y)/(y+13*x));
(x - 12*y)/(13*x + y)
? print([1,2;3,4])
[1, 2; 3, 4]
? print1(x+y);print(x+y);
x + yx + y
? print((x-12*y)/(y+13*x));
(x - 12*y)/(13*x + y)
? print([1,2;3,4])
[1, 2; 3, 4]
? print1(x+y);print1(" equals ");print(x+y);
x + y equals x + y
? print1("give a value for s? ");s=input();print(1/(s^2+1))
give a value for s? printtex((x+y)^3/(x-y)^2)
{{x^{3} + {{3}y}x^{2} + {{3}y^{2}}x + {y^{3}}}\over{x^{2} - {{2}y}x + {y^{2}
}}}
1
? for(i=1,100,for(j=1,25,if(i+j==32,break(2)));print(i))
1
2
3
4
5
6
? 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
og(p)/log(q))));
50 1.561229126903099279206171725
100 1.601335375590875348711141031
150 1.606915548673659127523394774
200 1.618659998991528481508175175
250 1.626284706204746765086080988
300 1.627822776845103001192024532
350 1.632105905172986668189652273
400 1.632424285532931448171405619
? install(addii,GG)
? addii(1,2)
3
? kill(addii)
? getheap
[21, 2993]
? print("Total time spent: ",gettime);
Total time spent: 1138
? \q