read `Z355.txt`:

L:=Z355:


#xi(s): the left side of the FE
xi:=proc(s):
s*(s-1)*Pi^(-s/2)*GAMMA(s/2)*Zeta(s):
end:

CheckD:=proc(s,N) local B,gu,i:

B:=evalf(1+gamma/2-log(2*sqrt(Pi))):

gu:=exp(-B*s):

gu:=evalf(gu*mul( (1-s/(1/2+I*L[i]))*exp (s/(1/2+I*L[i])),i=1..N))  :
gu:=evalf(gu*mul( (1-s/(1/2-I*L[i]))*exp (s/(1/2-I*L[i])),i=1..N))  :

evalf(gu/xi(s)):
end: