? allocatemem(20*10^6);
? gettime;oo=[1];
? \p96
   realprecision = 96 significant digits
? intcirc(s=1,0.5,zeta(s))-1
-1.0018787399041304420 E-94 - 6.296959716514689584 E-98*I
? intlaplaceinv(x=2,1,1/x)-1
7.37240875280988848267836598625171458524622502643426332739 E-55 + 1.71877297
0049496327338218359953890598325317714312174299845 E-54*I
? m=intnumstep();
? intlaplaceinv(x=2,1,1/x,m+1)-1
1.4981364335015034647 E-95 + 8.258510818309625628 E-95*I
? intlaplaceinv(x=5,1,1/x)-1
-5.985125294502659644 E-85 + 8.085986448700465906 E-85*I
? intlaplaceinv(x=100,1,1/x)-1
-2.0471231058052011760 E-55 - 3.351221763351295415 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.072350107766755152 E-97 + 9.907178318246741303 E-100*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
1.01780270266616498382717300808289765224 E-76
? intnum(x=0,[oo,1],g(x))-Euler
0.E-96
? intnum(x=0,1,1/sqrt(x))-2
-1.9262844386075338042 E-59
? intnum(x=[0,-1/2],1,1/sqrt(x))-2
0.E-96
? intnum(x=0,oo,sin(x)/x)-Pi/2
20.778594771898826028
? intnum(x=0,[oo,1],sin(x)/x)-Pi/2
0.00429640980192847594663632193954407532716512471317646262273262828586471032
064716268232582299149209
? 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.363352709384396655 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.3352846492208801743 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
0.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.7790370147830353644 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.0040645786513746241 E-306 + 0.E-330*I
? sumnum(n=1,2,1/(n^3+n+1),tab,1)-a
1.0040645786513746241 E-306
? c=sumnum(n=1,2,1/(n^2+1),tab,1);
? d=sumpos(n=1,1/(n^2+1));
? d-c
1.9691903647788732240 E-306
? sumnum(n=1,2,n^(-4/3),,1)-zeta(4/3)
-1.5810786927596923975695032330081859885853923189652965380063515132600239871
9095595253532786696377207848456709735455733654166843737571750244058595083139
34284299570579730419329071710335445546831167530844494561104268 E-110
? tab=sumnuminit([2,-3/2]);
? sumnum(n=1,[2,-3/2],1/(n*sqrt(n)),tab,1)-zeta(3/2)
-1.1904145143013527400 E-305
? suminf(n=1,2^(-n))-1
-1.1125369292536006916 E-308
? sumpos(n=1,2^(-n))-1
-2.781342323134001728 E-308
? sumnum(n=1,[2,log(2)],2^(-n),intnumstep()+1,1)-1
-1.9519460423754424134 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.3906711615670008645 E-308
? print("Total time spent: ",gettime);
Total time spent: 17637
? \q
