zeta(2)
Pi^2/6
zeta(-11)
bestappr(%)
-bernfrac(12)/12
zeta(1+x+O(x^5))
Euler

gamma(1/2)
sqrt(Pi)
gamma'(1)
gamma(s)*gamma(1-s)
Pi/sin(Pi*s)

lambda(s)=zeta(s)*Pi^-(s/2)*gamma(s/2);
lambda(I)
lambda(1-I)
lfunlambda(1,I)

zetamult([2,1])
zeta(3)
zetamult([3,5])+zetamult([5,3])+zeta(8)
zeta(3)*zeta(5)
a = zetamult([2,2,2])
b = zetamult([3,3])
bestappr(zeta(6)/a)
bestappr(zeta(6)/b)

c = zetamult([3,2,1])
lindep([a,b,c])
59*a-54*b+9*c
V=zetamultall(6);
V[zetamultconvert([3,2,1],2)]

dilog(1)
dilog(I)
-Pi^2/48+I*Catalan
polylog(3,1)

gamma(1/4 + O(5^10))
algdep(%,4)

zeta(-3+O(5^20))
bestappr(%)
(5^3-1)*bernfrac(4)/4

real(zetahurwitz(2,I))
(-Pi^2/sinh(Pi)^2-1)/2
zetahurwitz(0,1,1)
-log(2*Pi)/2
zetahurwitz(-3+O(5^20),1)

?

?4
?atan

??atan