? allocatemem(20*10^6);
  ***   Warning: new stack size = 20000000 (19.073 Mbytes).
? check(a,b)=my(t=abs((a-b)/b));if(t,ceil(log(t)/log(10)),"-oo");
? oo=[1];
? \p96
   realprecision = 96 significant digits
? check(intcirc(s=1,0.5,zeta(s)),1)
-93
? check(intlaplaceinv(x=2,1,1/x),1)
-53
? m=intnumstep();
? check(intlaplaceinv(x=2,1,1/x,m+1),1)
-94
? check(intlaplaceinv(x=5,1,1/x),1)
-83
? check(intlaplaceinv(x=100,1,1/x),1)
-52
? A=intmellininv(s=2,4,gamma(s)^3);
? tab=intfuncinit(t=[-oo,4.5],[oo,4.5],gamma(2+I*t)^3,1);
? check(intmellininvshort(2,4,tab),A)
-94
? 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));
? check(intnum(x=0,[oo,1],f(x)),Euler)
"-oo"
? check(intnum(x=0,[oo,1],g(x)),Euler)
"-oo"
? check(intnum(x=0,1,1/sqrt(x)),2)
-59
? check(intnum(x=[0,-1/2],1,1/sqrt(x)),2)
"-oo"
? check(intnum(x=0,[oo,1],sin(x)/x),Pi/2)
-2
? check(intnum(x=0,[oo,-I],sin(x)/x),Pi/2)
"-oo"
? check(intnum(x=0,[oo,-2*I],sin(2*x)/x),Pi/2)
"-oo"
? A=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);
? check(A,Pi/2)
-96
? check(intnum(x=0,[oo,1],sin(x)^3*exp(-x)),3/10)
-96
? check(intnum(x=0,[oo,-I],sin(x)^3*exp(-x)),3/10)
-88
? tab=intnuminit(0,[oo,-I],m+1);
? check(intnum(x=0,oo,sin(x)^3*exp(-x),tab),3/10)
-96
? check(intnum(x=0,[oo,-I],x^2*sin(x)),-2)
"-oo"
? tab=intnuminit(-1,1);
? check(intnum(x=-1,1,intnum(y=-sqrt(1-x^2),sqrt(1-x^2),x^2+y^2,tab),tab),Pi/2)
-93
? \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);
? check(a,b)
-305
? check(sumnum(n=1,2,1/(n^3+n+1),tab,1),a)
-305
? c=sumnum(n=1,2,1/(n^2+1),tab,1);
? d=sumpos(n=1,1/(n^2+1));
? check(c,d)
-305
? check(sumnum(n=1,2,n^(-4/3),,1),zeta(4/3))
-110
? tab=sumnuminit([2,-3/2]);
? check(sumnum(n=1,[2,-3/2],1/(n*sqrt(n)),tab,1),zeta(3/2))
-305
? check(suminf(n=1,2^(-n)),1)
-307
? check(sumpos(n=1,2^(-n)),1)
-307
? check(sumnum(n=1,[2,log(2)],2^(-n),intnumstep()+1,1),1)
-304
? 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));
? check(a,b)
-307
? \p96
   realprecision = 96 significant digits
? T=intnuminitgen(t,0,[1],exp(2*sinh(t)));
? check(intnum(x=0,[1],1/(1+x^2),T),Pi/2)
"-oo"
? T=intnuminitgen(t,0,[[1],1],exp(t-exp(-t)));
? check(intnum(x=0,[[1],1],exp(-x),T),1)
"-oo"
? intfourierexp(t=0,[[1],1],1/2,exp(-t^2))
0.07515645001618094506724269142337819200573257288846457257682421833136020654
01315089580450881170413 - 0.408148557374490719018171389392074027664154880556
033855962874575972355801042051564345987971241498*I
? intfouriercos(t=0,[[1],1],1/2,exp(-t^2))
0.07515645001618094506724269142337819200573257288846457257682421833136020654
01315089580450881170413
? intfouriersin(t=0,[[1],1],1/2,exp(-t^2))
0.40814855737449071901817138939207402766415488055603385596287457597235580104
2051564345987971241498
? \p38
   realprecision = 38 significant digits
? intnumromb(x=0,1,sin(x))
0.45969769413186028259906339255702339518
? intnumromb(x=0,1,sin(x),1)
0.45969769413186028259906339255702338970
? intnumromb(x=1,100,exp(-x^2),2)
0.13940279264033098824961630553871957887
? intnumromb(x=0,1,sin(x)/x,3)
0.94608307036718301494135331382317964743
? f(x)=-log(cos(x));
? F=truncate(f(t+O(t^16)));
? g(x)=if(x>1e-2,f(x),subst(F,t,x));
? sumpos(n=1,g(1/n))
0.94536905472633293526609521540827019810
? print("Total time spent: ",gettime);
Total time spent: 20736
