default(echo,1);
allocatemem(20 * 10^6);
gettime; oo = [1];
\p96
intcirc(s=1, 0.5, zeta(s)) - 1
intlaplaceinv(x=2, 1, 1/x) - 1
m = intnumstep();
intlaplaceinv(x=2, 1, 1/x, m+1) - 1
intlaplaceinv(x=5, 1, 1/x) - 1
intlaplaceinv(x=100, 1, 1/x) - 1
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
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
intnum(x = 0, [oo,1],  g(x)) - Euler
intnum(x = 0, 1,  1/sqrt(x)) - 2
intnum(x = [0,-1/2], 1,  1/sqrt(x)) - 2
intnum(x = 0, [oo,1], sin(x)/x) - Pi/2
intnum(x = 0, [oo,-I], sin(x)/x) - Pi/2
intnum(x = 0, [oo,-2*I], sin(2*x)/x) - Pi/2
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

intnum(x = 0, [oo, 1], sin(x)^3*exp(-x)) - 0.3
intnum(x = 0, [oo,-I], sin(x)^3*exp(-x)) - 0.3
tab = intnuminit(0,[oo,-I], m+1);
intnum(x = 0, oo, sin(x)^3*exp(-x), tab) - 0.3
intnum(x = 0, [oo, -I], x^2*sin(x)) + 2
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
\\
\p 308
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
sumnum(n=1,2, 1/(n^3+n+1), tab, 1) - a
c = sumnum(n=1,2,1/(n^2+1),tab,1);
d = sumpos(n=1,1/(n^2+1));
d - c
sumnum(n=1,2,n^(-4/3),,1) - zeta(4/3)
tab = sumnuminit([2,-3/2]);
sumnum(n=1,[2,-3/2],1/(n*sqrt(n)),tab,1)-zeta(3/2)
\\
suminf(n=1, 2^(-n)) - 1
sumpos(n=1, 2^(-n)) - 1
sumnum(n=1,[2,log(2)],2^(-n), intnumstep()+1, 1) - 1
\\
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
