? allocatemem(20*10^6);
? gettime;oo=[1];
? \p96
   realprecision = 96 significant digits
? intcirc(s=1,0.5,zeta(s))-1
-1.001878740 E-94 - 1.227730322 E-97*I
? intlaplaceinv(x=2,1,1/x)-1
7.37240875280988848267836598625171458524622502643 E-55 + 1.71877297004949632
7338218359953890598325317714312 E-54*I
? m=intnumstep();
? intlaplaceinv(x=2,1,1/x,m+1)-1
1.498136434 E-95 + 8.25851081 E-95*I
? intlaplaceinv(x=5,1,1/x)-1
-5.985125294502659643 E-85 + 8.085986448700465906 E-85*I
? intlaplaceinv(x=100,1,1/x)-1
-2.047123106 E-55 - 3.351417365 E-53*I
? A=intmellininv(s=2,4,gamma(s)^3);
? tab=intfuncinit(t=[-oo,4.5],[oo,4.5],gamma(2+I*t)^3,1);
? intmellininvshort(2,4,tab)-A
-3.072350108 E-97 + 1.323681855 E-99*I
? f(x)=1/(exp(x)-1)-exp(-x)/x
? F=truncate(f(t+O(t^7)));
? g(x)=if(x>1e-18,f(x),subst(F,t,x))
? intnum(x=0,[oo,1],f(x))-Euler
3.7783741346110721665958525482 E-67
? intnum(x=0,[oo,1],g(x))-Euler
0.E-96
? intnum(x=0,1,1/sqrt(x))-2
-1.4009863823182375326 E-54
? intnum(x=[0,-1/2],1,1/sqrt(x))-2
0.E-96
? intnum(x=0,oo,sin(x)/x)-Pi/2
11.84217614
? intnum(x=0,[oo,1],sin(x)/x)-Pi/2
0.00448267113448250456172244085159836765958193607953130895411895619311247452
425300891424918944797442
? intnum(x=0,[oo,-I],sin(x)/x)-Pi/2
0.E-96
? intnum(x=0,[oo,-2*I],sin(2*x)/x)-Pi/2
0.E-96
? intnum(x=0,1,(1-cos(x))/x^2)+intnum(x=1,oo,1/x^2)-intnum(x=1,[oo,I],cos(x)
/x^2)-Pi/2
-9.36335271 E-97
? intnum(x=0,[oo,1],sin(x)^3*exp(-x))-0.3
0.E-96
? intnum(x=0,[oo,-I],sin(x)^3*exp(-x))-0.3
-1.335284649 E-89
? tab=intnuminit(0,[oo,-I],m+1);
? intnum(x=0,oo,sin(x)^3*exp(-x),tab)-0.3
0.E-96
? intnum(x=0,[oo,-I],x^2*sin(x))+2
5.618011626 E-96
? tab=intnuminit(-1,1);
? intnum(x=-1,1,intnum(y=-sqrt(1-x^2),sqrt(1-x^2),x^2+y^2,tab),tab)-Pi/2
-1.779037015 E-94
? \p308
   realprecision = 308 significant digits
? a=sumpos(n=1,1/(n^3+n+1));
? tab=sumnuminit(2);
? b=sumnum(n=1,2,1/(n^3+n+1),tab);
? a-b
-1.004064578 E-306 + 0.E-320*I
? sumnum(n=1,2,1/(n^3+n+1),tab,1)-a
1.004064578 E-306
? c=sumnum(n=1,2,1/(n^2+1),tab,1);
? d=sumpos(n=1,1/(n^2+1));
? d-c
1.969190365 E-306
? sumnum(n=1,2,n^(-4/3),,1)-zeta(4/3)
-2.4216544207504506114899981581969759071445497499100570634176004383007422693
6954882413244126542269381152110599476350651837206116316883640508880554373752
13154920980196906076673594248661407504915027261660100 E-107
? tab=sumnuminit([2,-3/2]);
? sumnum(n=1,[2,-3/2],1/(n*sqrt(n)),tab,1)-zeta(3/2)
-1.190414514 E-305
? suminf(n=1,2^(-n))-1
-1.112536929 E-308
? sumpos(n=1,2^(-n))-1
-2.781342324 E-308
? sumnum(n=1,[2,log(2)],2^(-n),intnumstep()+1,1)-1
-1.950277237 E-305
? tab=sumnuminit(2,,-1);
? a=sumnumalt(n=1,2,1/(n^3+n+1),tab,1);
? b=sumalt(n=1,(-1)^n/(n^3+n+1));
? a-b
-1.390671162 E-308
? print("Total time spent: ",gettime);
Total time spent: 39618
? \q
