realprecision = 19 significant digits
   echo = 1 (on)
? intnum(x=0,Pi,sin(x),1)
2.000000000000000017
? sqr(2*intnum(x=0,4,exp(-x^2),1))
3.141592556720305685
? 4*intnum(x=1,10^20,1/(1+x^2),2)
3.141592653589793208
? intnum(x=-0.5,0.5,1/sqrt(1-x^2))
1.047197551196597747
? 2*intnum(x=0,100,sin(x)/x,3)
3.124450933778112629
? \p38
   realprecision = 38 significant digits
? prod(k=1,10,1+1/k!)
3335784368058308553334783/905932868585678438400000
? prod(k=1,10,1+1./k!)
3.6821540356142043935732308433185262946
? Pi^2/6*prodeuler(p=2,10000,1-p^-2)
1.0000098157493066238697591433298145174
? prodinf(n=0,(1+2^-n)/(1+2^(-n+1)))
0.33333333333333333333333333333333333322
? prodinf(n=0,-2^-n/(1+2^(-n+1)),1)
0.33333333333333333333333333333333333322
? solve(x=1,4,sin(x))
3.1415926535897932384626433832795028841
? sum(k=1,10,2^-k)
1023/1024
? sum(k=1,10,2.^-k)
0.99902343750000000000000000000000000000
? 4*sumalt(n=0,(-1)^n/(2*n+1))
3.1415926535897932384626433832795028841
? 4*sumalt(n=0,(-1)^n/(2*n+1),1)
3.1415926535897932384626433832795028842
? sumdiv(8!,x,x)
159120
? suminf(n=1,2.^-n)
0.99999999999999999999999999999999999999
? 6/Pi^2*sumpos(n=1,n^-2)
0.99999999999999999999999999999999999999
? getheap
[19, 190]
? print("Total time spent: ",gettime);
Total time spent: 1646
? \q