(* Mathematica script 2danQCDcouplrt0201al01189.m *)

(* Copyright: G. Cvetic (November 2018), Univ. Tecnica Federico Santa Maria
 Permission is granted for anyone to use this script for any purpose 
 on any computer. The script can be freely distributed. *)

(* The script has been checked to work in Mathematica 11.3 *)
(* ************************************************************ *)

Print["2danQCDcouplrt0201al01189: a Mathematica package for analytic running"]
Print["of the strong coupling A (the analog of pQCD coupling a=alpha_s/Pi)"]
Print["in a model with 2 deltas in the spectral function of A,"]
Print["in the c2=-4.9 Lambert renormalization scheme."]
Print["This model:"]
Print["1.) agrees with perturbative QCD at high Q2 up to (Lambda^2/Q2)^5;"]
Print["2.) reproduces correctly the tau semihadronic decay physics."]
Print["The input parameters are such that r_{tau}^{(D=0)}=0.201 ('rt0201')"]
Print["in the (leading-beta0) + (beyond-leading-beta0) resummation (rtauLBbLB)"]
Print["and alpha_s(MZ^2;MSbar)=0.1189 ('al01189')."]
Print["For additional details, see arXiv:1408.6868"]
Print["The software has been checked to work in: Mathematica 11.3"]
Print["The software can be freely distributed."];
Print["Copyright: G.Cvetic, September 2018, Univ. Tecnica Federico Santa Maria"]

(* ************************************************************ *)
(* ************************************************************ *)


(* ::Section:: *)
(*Definitions*)

(* Couplings, functions: *)
LL2::usage= "LL2[Nf]=Lambda_L^2, this is square the Lambert Lambda scale Lambda_L (with Nf flavors), in GeV^2, appearing in the pQCD running coupling a(Q^2) (see: a2d) underlying for the considered Nf-flavor 2danQCD";
alphaMSMZ::usage= "alphaMSMZ=alpha_s(Q^2=MZ^2,MSbar). In this case, alphaMSMZ=0.1189"
c22din::usage= "c22din=c2=beta_2/beta_0 (chosen to be c2ddin=-4.9) is the value of the c2 coefficient (for Nf=3,4,5,6) in the Lambert renormalization scheme used in the considered 2danQCD [c2=-4.9; c3=c2^2/c1. c4=c2^3/c1^2, etc., and cj=beta_j/beta_0; normalization convention: beta_0=(11 -2*Nf/3)/4; beta_1=(102-38*Nf/3)/16; etc.]"
s0::usage= "s0[Nf]=s_0 is a dimensionless parameter for Nf flavor 2danQCD (Nf=3,4,5,6); it is equal for Nf=3,4,5,6; the pQCD-onset scale M_0^2 in Nf=3 2danQCD is then: M_0^2(Nf=3)=s0[3]*LL2[3] (in GeV^2), and similarly for other cases of Nf"
s1s0::usage= "s1s0[Nf]=s_1 is a dimensionless parameter for Nf flavor 2danQCD (Nf=3,4,5,6); the location of the first delta function of the discontinuity function rho_1(sigma)=Im A(Q^2=-sigma - I epsilon) of the A-coupling of the Nf-flavor 2danQCD is then at sigma=M_1^2=s1s0[Nf]*LL2[Nf] (in GeV^2)"
s2s0::usage= "s2s0[Nf]=s_2 is a dimensionless parameter for Nf flavor 2danQCD (Nf=3,4,5,6); the location of the second delta function of the discontinuity function rho_1(sigma)=Im A(Q^2=-sigma - I epsilon) of the A-coupling of the Nf-flavor 2danQCD is then at sigma=M_2^2=s2s0[Nf]*LL2[Nf] (in GeV^2)"
f1ss0::usage= "f1ss0[Nf]=f_1^2 is a dimensionless strength parameter at the first delta in rho_1(sigma) discontinuity function of A(Q^2), in the Nf-flavor 2danQCD; the first delta function in rho_1(sigma) is Pi*f1ss0[Nf]*LL2[Nf]*delta(sigma - M1^2[Nf]), where M1^2[Nf]= s1s0[Mf]*LL2[Nf] (in GeV^2)"
f2ss0::usage= "f2ss0[Nf]=f_2^2 is a dimensionless strength parameter at the second delta in rho_1(sigma) discontinuity function of A(Q^2), in the Nf-flavor 2danQCD; the second delta function in rho_1(sigma) is Pi*f2ss0[Nf]*LL2[Nf]*delta(sigma - M2^2[Nf]), where M2^2[Nf]= s2s0[Mf]*LL2[Nf] (in GeV^2)"  
a2d::usage= "a2d[Nf,u,fi]= alpha_s(Q^2=u*LL2[Nf]*Exp[I fi)/Pi =a(Q^2), where Q^2=|Q^2|*Exp[I fi] and |Q^2|=u*LL2[Nf] (see: LL2); a2d is the underlying pQCD coupling in the Lambert renormalization scheme [c2=c22din=-4.9; c3=c2^2/c1. c4=c2^3/c1^2, etc., and cj=beta_j/beta_0; normalization conventions: beta_0=(11 -2*Nf/3)/4; beta_1=(102-38*Nf/3)/16; etc.], the pQCD coupling underlying the Nf-flavor 2danQCD coupling A (see: A1).";
apert::usage= "apert[Nf,Q2,fi]= alpha_s(Q^2=Q2*Exp[I fi)/Pi =a(Q^2), where Q^2=|Q^2|*Exp[I fi] and |Q^2|=Q2; apert is the underlying pQCD coupling in the Lambert renormalization scheme [c2=c22din=-4.9; c3=c2^2/c1. c4=c2^3/c1^2, etc., and cj=beta_j/beta_0; normalization conventions: beta_0=(11 -2*Nf/3)/4; beta_1=(102-38*Nf/3)/16; etc.], the pQCD coupling underlying the Nf-flavor 2danQCD coupling A (see: A1); note, apert[Nf,Q2,fi]=a2d[Nf,Q2/LL2[Nf],fi].";
tr12d::usage= "tr12d[Nf,s] is the discontinuity function rho_1^{(pt)}(sigma)=Im a(Q^2=-sigma - I epsilon) of the underlying pQCD coupling a(Q^2) (see: a2d and apert) in the considered Nf-flavor 2danQCD, where s=sigma/LL2[Nf] (see: LL2)"
A1::usage= "A1[Nf,Q2,fi] = A(Q^2) = A_1(Q^2), where Q^2=Q2 Exp[I fi] (Q2=|Q^2|, in GeV^2) and A(Q^2) is the 2danQCD analog of the underlying perturbative coupling a(Q^2) = a2d[Nf,Q2/LL2[Nf],fi] (see: a2d); at high |Q^2|> 1 GeV^2 the two couplings practically coincide";
tA2::usage= "tA2[Nf,Q2,fi]={tilde A}_2(Q^2) for Q^2=Q2 Exp[I fi] (Q2=|Q^2|, in GeV^2); this is the first logarithmic derivative of A_1 (see: A1): {tilde A}_2(Q^2) = (-1)/(beta_0 1!) (d/d ln Q^2) A_1(Q^2)";
tA3::usage= "tA3[Nf,Q2,fi]={tilde A}_3(Q^2) for Q^2=Q2 Exp[I fi]; this is the second logarithmic derivative of A_1 (see: A1): {tilde A}_3(Q^2) = (-1)^2/(beta_0^2 2!) (d/d ln Q^2)^2 A_1(Q^2)";
tA4::usage= "tA4[Nf,Q2,fi]={tilde A}_4(Q^2) for Q^2=Q2 Exp[I fi]; this is the third logarithmic derivative of A_1 (see: A1): {tilde A}_4(Q^2) = (-1)^3/(beta_0^3 3!) (d/d ln Q^2)^3 A_1(Q^2)";
tA2d::usage= "tA2d[Nf,nu,Q2,fi]= {tilde A}_{nu}(Q^2), with Q^2=Q2 Exp[I fi]  (Q2=|Q^2|); {tilde A}_{nu} can have nu noninteger; {tilde A}_{nu} is the generalization (to noninteger nu) of the logarithmic derivative (-1)^{nu-1}/(beta_0^{nu-1} (nu-1)! (d/d ln Q^2)^(nu-1) A_1(Q^2)";
A2dNX::usage= "A2dNX[Nf,M,nu,Q2,fi] = A_{nu+M}(Q^2) is the 2danQCD analog of the power of the pQCD coupling a2d[Nf,Q2/LL2[Nf],fi]^(nu+M)=a(Q^2)^(nu+M) (see: a2d), for NX=N1,N2,N3,N4; here, -1 =< nu =< 1 (in general noninteger), and M is integer, M=0,1,..,X; where Q^2=Q2 Exp[I fi] and Q2=|Q^2|=u*LL2[nq] (in GeV^2). Note that for X=1 (N1), A2dN1[...] are the couplings (A_{nu}, A_{nu+1}) in the NLO sense, when in the TPS of a considered observable we only have two terms known (~a^nu and ~a^(nu+1)); for X=2 (N2), A2dN2[...] are the couplings (A_{nu}, A_{nu+1}, A_{nu+2}) in the N^2LO sense, i.e.,  when in the truncated perturbation series of a considered observable we have three terms known (~a^nu, ~a^(nu+1) and ~a^(nu+2)); etc.; Nf indicates the chosen number of active flavors of the 2danQCD (Nf=3,4,5,6), usually Nf=3."
AYNX::usage= "AYNX[Nf,Q2,fi] = A_{Y}(Q^2) is the 2danQCD analog of the power of the pQCD coupling a2d[Nf,Q2/LL2[Nf],fi]^Y=a(Q^2)^Y, for NX=N1,N2,N3,N4; here, Y is integer, Y=1,..,X+1; where Q^2=Q2*Exp[I fi] and Q2=|Q^2|= u*LL2[nq] (in GeV^2). Note that X=1 (N1) is used when we have in the truncated perturbation series of a considered observable only two terms (~a and ~a^2) known; X=2 (N2) is used when three terms are known (~a, ~a^2 and ~a^3); X=3 (N3) is used when four terms are known (~a, ~a^2, ~a^3 and ~a^4); X=4 (N4) is used when five terms are known (~a, ~a^2, ~a^3, ~a^4, and ~a^5). Nf indicates the chosen number of active flavors of the 2danQCD (Nf=3,4,5,6), usually Nf=3."
tdj::usage= "tdj={tilde d}_j (j=1,2,3) are the coefficients for the Nf=3 canonical (and massless) Adler function in the Lambert renormalization scheme used in the considered 2danQCD and in the underlying pQCD [c2=-4.9; c3=c2^2/c1. c4=c2^3/c1^2, etc., and cj=beta_j/beta_0; normalizations: beta_0=(11 -2 Nf/3)/4; beta_1=(102-38 Nf/3)/16; etc.]"
tdjbLB::usage= "tdjbLB={tilde d}_j - {tilde d}_j^(LB) are the 'beyond-the-leading-beta0' (bLB) part of {tilde d}_j Nf=3 coefficients of the canonical (and massless) Adler function in the Lambert renormalization scheme used in the considered 2danQCD and in the underlying pQCD (see also: tdj)."
rtauLBbLB::usage= "rtauLBbLB=r_{tau}(D=0) calculated (in the considered Lambert renormalization scheme) in the considered 2danQCD with the 'leading-beta_0 (LB) + beyond-the leading-beta_0 (bLB)' approach with Nf=3. Note that here rtauLBbLB=0.2010"
rtaudBG22::usage= "rtaudBG22=r_{tau}(D=0) calculated in the considered 2danQCD with the diagonal [2/2] Baker-Gammel approach."
rtauTPS::usage= "rtauTPS=r_{tau}(D=0) calculated in the considered 2danQCD with the truncated perturbation series approach."
aMSnf3::usage= "aMSnf3[Q2] is the pQCD coupling a(Q^2=Q2)=alpha_s(Q2)/Pi for positive Q2 (Q2 in GeV^2) in MSbar renormalization scheme, for the number of active flavors Nf=3; it corresponds to alpha_s(MZ^2;MSbar)=0.1189, and is obtained by 4-loop RGE evolution in MSbar scheme from Q2=M_Z^2 downwards, with 3-loop quark threshold conditions imposed at threshold momenta Q2=(2*mqbar)^2, where mqbar=mqbar(Q2=mqbar^2) are the MSbar quark masses of b (for Nf=5->4) and of c (for Nf=4->3); we used mbbar=4.2 GeV and mcbar=1.27 GeV."

Off[NIntegrate::inumr]

s0[3] = 28.3
 
s0[4] = 28.3
 
s0[5] = 28.3
 
s0[6] = 28.3
 
s0 /: s0::usage = "s0[Nf]=s_0 is a dimensionless parameter for Nf flavor \
2danQCD (Nf=3,4,5,6); it is equal for Nf=3,4,5,6; the pQCD-onset scale M_0^2 \
in Nf=3 2danQCD is then: M_0^2(Nf=3)=s0[3]*LL2[3] (in GeV^2), and similarly \
for other cases of Nf"
 
s1s0[3] = 20.748682896120894
 
s1s0[4] = 20.56125283859566
 
s1s0[5] = 20.201356561403315
 
s1s0[6] = 19.281901615295563
 
s1s0 /: s1s0::usage = "s1s0[Nf]=s_1 is a dimensionless parameter for Nf \
flavor 2danQCD (Nf=3,4,5,6); the location of the first delta function of the \
discontinuity function rho_1(sigma)=Im A(Q^2=-sigma - I epsilon) of the \
A-coupling of the Nf-flavor 2danQCD is then at sigma=M_1^2=s1s0[Nf]*LL2[Nf] \
(in GeV^2)"
 
s2s0[3] = 1.3263540157512053
 
s2s0[4] = 0.5691199334573928
 
s2s0[5] = -1.0828094755190396
 
s2s0[6] = -7.113846642305069
 
s2s0 /: s2s0::usage = "s2s0[Nf]=s_2 is a dimensionless parameter for Nf \
flavor 2danQCD (Nf=3,4,5,6); the location of the second delta function of the \
discontinuity function rho_1(sigma)=Im A(Q^2=-sigma - I epsilon) of the \
A-coupling of the Nf-flavor 2danQCD is then at sigma=M_2^2=s2s0[Nf]*LL2[Nf] \
(in GeV^2)"
 
f1ss0[3] = 0.3117123419801963
 
f1ss0[4] = 0.400635450969464
 
f1ss0[5] = 0.5730950435870608
 
f1ss0[6] = 0.9828910659783962
 
f1ss0 /: f1ss0::usage = "f1ss0[Nf]=f_1^2 is a dimensionless strength \
parameter at the first delta in rho_1(sigma) discontinuity function of \
A(Q^2), in the Nf-flavor 2danQCD; the first delta function in rho_1(sigma) is \
Pi*f1ss0[Nf]*LL2[Nf]*delta(sigma - M1^2[Nf]), where M1^2[Nf]= \
s1s0[Mf]*LL2[Nf] (in GeV^2)"
 
f2ss0[3] = 0.5997353984783986
 
f2ss0[4] = 0.8408539911987749
 
f2ss0[5] = 1.4344402833120198
 
f2ss0[6] = 4.121439967426389
 
f2ss0 /: f2ss0::usage = "f2ss0[Nf]=f_2^2 is a dimensionless strength \
parameter at the second delta in rho_1(sigma) discontinuity function of \
A(Q^2), in the Nf-flavor 2danQCD; the second delta function in rho_1(sigma) \
is Pi*f2ss0[Nf]*LL2[Nf]*delta(sigma - M2^2[Nf]), where M2^2[Nf]= \
s2s0[Mf]*LL2[Nf] (in GeV^2)"
 
M0[3] = 1.3884978587575423
 
M0[4] = 0.9898699722471362
 
M0[5] = 0.5244038264593068
 
M0[6] = 0.12758729827894974
 
LLamb[3] = 0.2610069048390281
 
LLamb[4] = 0.18607367380494796
 
LLamb[5] = 0.09857632747980138
 
LLamb[6] = 0.023983591771874337
 
LL2[3] = 0.06812460437364946
 
LL2[4] = 0.03462341208327018
 
LL2[5] = 0.009717292339405046
 
LL2[6] = 0.0005752126742799184
 
LL2 /: LL2::usage = "LL2[Nf]=Lambda_L^2, this is square the Lambert Lambda \
scale Lambda_L (with Nf flavors), in GeV^2, appearing in the pQCD running \
coupling a(Q^2) (see: a2d) underlying for the considered Nf-flavor 2danQCD"
 
be0[nq_] := (11 - 2*(nq/3))/4
 
c1[nq_] := (102 - 38*(nq/3))/(11 - 2*(nq/3))/4
 
c22d[nq_] := c22din
 
c22din = -4.9
 
c22din /: c22din::usage = "c22din=c2=beta_2/beta_0 (chosen to be c2ddin=-4.9) \
is the value of the c2 coefficient (for Nf=3,4,5,6) in the Lambert \
renormalization scheme used in the considered 2danQCD [c2=-4.9; c3=c2^2/c1. \
c4=c2^3/c1^2, etc., and cj=beta_j/beta_0; normalization convention: \
beta_0=(11 -2*Nf/3)/4; beta_1=(102-38*Nf/3)/16; etc.]"
 
c32d[nq_] := c22d[nq]^2/c1[nq]
 
c42d[nq_] := c22d[nq]^3/c1[nq]^2
 
c52d[nq_] := c22d[nq]^4/c1[nq]^3
 
c2MS[nq_] := be2MS[nq]/be0[nq]
 
be2MS[nq_] := (2857/2 - 5033*(nq/18) + 325*(nq^2/54))/64
 
c3MS[nq_] := be3MS[nq]/be0[nq]
 
be3MS[nq_] := (1/256)*((149753/6 + 3564*Zeta[3]) - 
      (1078361/162 + (6508/27)*Zeta[3])*nq + (50065/162 + (6472/81)*Zeta[3])*
       nq^2 + (1093/729)*nq^3)
 
btP11[x_, nq_] := betaP11[x, nq]
 
betaP11[x_, nq_] := (-be0[nq])*x^2*((1 + (c1[nq] - c22d[nq]/c1[nq])*x)/
      (1 - (c22d[nq]/c1[nq])*x))
 
betaMS4l[x_, nq_] := (-be0[nq])*x^2*(1 + c1[nq]*x + c2MS[nq]*x^2 + 
      c3MS[nq]*x^3)
 
kappa2d = 2.
 
alphaMSMZ = 0.1189
 
alphaMSMZ /: alphaMSMZ::usage = "alphaMSMZ=alpha_s(Q^2=MZ^2,MSbar). In this \
case, alphaMSMZ=0.1185, or 0.1181, or 0.1189"
 
apert[nq_, Q2_, fi_] := a2d[nq, Q2/LL2[nq], fi]
 
a2d[nq_, u_, fi_] := If[fi >= 0, (-(1/c1[nq]))*
      (1/(1 - c22d[nq]/c1[nq]^2 + ProductLog[-1, zu2d[nq, u, fi]])), 
     (-(1/c1[nq]))*(1/(1 - c22d[nq]/c1[nq]^2 + ProductLog[1, 
         zu2d[nq, u, fi]]))]
 
a2d /: a2d::usage = "a2d[Nf,u,fi]= alpha_s(Q^2=u*LL2[Nf]*Exp[I fi)/Pi \
=a(Q^2), where Q^2=|Q^2|*Exp[I fi] and |Q^2|=u*LL2[Nf] (see: LL2); a2d is the \
underlying pQCD coupling in the Lambert renormalization scheme \
[c2=c22din=-4.9; c3=c2^2/c1. c4=c2^3/c1^2, etc., and cj=beta_j/beta_0; \
normalization conventions: beta_0=(11 -2*Nf/3)/4; beta_1=(102-38*Nf/3)/16; \
etc.], the pQCD coupling underlying the Nf-flavor 2danQCD coupling A (see: \
A1)."
 
zu2d[nq_, u_, fi_] := ((-(1/c1[nq]/E))*Exp[(-I)*(be0[nq]/c1[nq])*fi])/
     u^(be0[nq]/c1[nq])
 
aMSnf3[Qf2_] := First[y[Log[Qf2/Q4^2]] /. sol1a]
 
Q4 = 2.54
 
sol1a = {{y -> InterpolatingFunction[{{-2.8215360227571686, 20.}}, 
        {5, 7, 1, {656}, {4}, 0, 0, 0, 0, Automatic, {}, {}, False}, 
        {{-2.8215360227571686, -2.8215360227571264, -2.8215360227570843, 
         -2.8215360227570416, -2.8215360227569994, -2.8215360227569484, 
         -2.8215360227568973, -2.8215360227568462, -2.8215360227567956, 
         -2.821536022756734, -2.8215360227566726, -2.8215360227566113, 
         -2.82153602275655, -2.821536022756476, -2.821536022756402, 
         -2.821536022756328, -2.8215360227562543, -2.8215360227561654, 
         -2.821536022756076, -2.8215360227559874, -2.821536022755898, 
         -2.821536022755791, -2.821536022755684, -2.821536022755577, 
         -2.82153602275547, -2.8215360227553408, -2.821536022755212, 
         -2.8215360227550828, -2.821536022754954, -2.8215360227547985, 
         -2.8215360227546435, -2.821536022754488, -2.821536022754333, 
         -2.821536022754146, -2.821536022753959, -2.8215360227537722, 
         -2.8215360227535853, -2.82153602275336, -2.821536022753135, 
         -2.82153602275291, -2.821536022752685, -2.8215360227524138, 
         -2.821536022752143, -2.821536022751872, -2.821536022751601, 
         -2.8215360227512747, -2.8215360227509487, -2.8215360227506223, 
         -2.821536022750296, -2.8215360227499033, -2.8215360227495103, 
         -2.8215360227491177, -2.8215360227487247, -2.8215360227482518, 
         -2.821536022747779, -2.821536022747306, -2.8215360227468333, 
         -2.8215360227462636, -2.8215360227456943, -2.821536022745125, 
         -2.821536022744555, -2.8215360227438695, -2.821536022743184, 
         -2.8215360227424986, -2.821536022741813, -2.8215360227409874, 
         -2.821536022740162, -2.8215360227393362, -2.8215360227385102, 
         -2.8215360227375164, -2.8215360227365225, -2.821536022735528, 
         -2.8215360227345343, -2.821536022733337, -2.82153602273214, 
         -2.821536022730943, -2.821536022729746, -2.821536022728305, 
         -2.8215360227268635, -2.821536022725422, -2.821536022723981, 
         -2.821536022722245, -2.8215360227205095, -2.821536022718774, 
         -2.8215360227170385, -2.821536022714948, -2.8215360227128583, 
         -2.8215360227107684, -2.8215360227086785, -2.8215360227061614, 
         -2.821536022703645, -2.821536022701128, -2.821536022698611, 
         -2.8215360226955806, -2.8215360226925497, -2.8215360226895188, 
         -2.821536022686488, -2.821536022682838, -2.821536022679188, 
         -2.8215360226755384, -2.8215360226718884, -2.821536022667493, 
         -2.821536022663097, -2.8215360226587016, -2.821536022654306, 
         -2.8215360226490125, -2.821536022643719, -2.8215360226384254, 
         -2.8215360226331314, -2.821536022626756, -2.8215360226203807, 
         -2.8215360226140054, -2.82153602260763, -2.8215360225999517, 
         -2.8215360225922734, -2.8215360225845947, -2.8215360225769164, 
         -2.8215360225676687, -2.8215360225584205, -2.8215360225491724, 
         -2.8215360225399246, -2.821536022528786, -2.821536022517647, 
         -2.8215360225065083, -2.8215360224953696, -2.8215360224819532, 
         -2.821536022468537, -2.8215360224551205, -2.8215360224417037, 
         -2.8215360224255432, -2.821536022409383, -2.8215360223932224, 
         -2.821536022377062, -2.821536022357596, -2.821536022338129, 
         -2.821536022318663, -2.8215360222991968, -2.8215360222757475, 
         -2.8215360222522983, -2.821536022228849, -2.8215360222054, 
         -2.8215360221771513, -2.821536022148903, -2.821536022120655, 
         -2.821536022092407, -2.8215360220583765, -2.821536022024346, 
         -2.8215360219903154, -2.821536021956285, -2.8215360219152865, 
         -2.821536021874288, -2.8215360218332903, -2.821536021792292, 
         -2.8215360217428977, -2.8215360216935035, -2.8215360216441088, 
         -2.8215360215947145, -2.8215360215352017, -2.821536021475689, 
         -2.821536021416176, -2.821536021356663, -2.821536021284956, 
         -2.8215360212132485, -2.8215360211415415, -2.821536021069834, 
         -2.8215360209834293, -2.8215360208970246, -2.8215360208106204, 
         -2.8215360207242157, -2.821536020620096, -2.8215360205159765, 
         -2.8215360204118567, -2.8215360203077373, -2.821536020182264, 
         -2.821536020056791, -2.821536019931318, -2.8215360198058446, 
         -2.8215360196546304, -2.821536019503416, -2.821536019352201, 
         -2.821536019200987, -2.82153601901874, -2.821536018836493, 
         -2.8215360186542457, -2.821536018471998, -2.8215360182523366, 
         -2.821536018032675, -2.821536017813013, -2.8215360175933513, 
         -2.821536017328577, -2.821536017063803, -2.8215360167990284, 
         -2.821536016534254, -2.8215360162150813, -2.8215360158959086, 
         -2.8215360155767355, -2.821536015257563, -2.8215360148727884, 
         -2.8215360144880144, -2.82153601410324, -2.8215360137184655, 
         -2.8215360132545726, -2.8215360127906792, -2.8215360123267863, 
         -2.8215360118628934, -2.8215360113035706, -2.8215360107442473, 
         -2.8215360101849245, -2.8215360096256017, -2.821536008951164, 
         -2.821536008276726, -2.821536007602288, -2.8215360069278503, 
         -2.8215360061145383, -2.8215360053012257, -2.8215360044879136, 
         -2.821536003674601, -2.8215360026937333, -2.821536001712865, 
         -2.8215360007319967, -2.821535999751129, -2.821535998568078, 
         -2.821535997385027, -2.821535996201976, -2.8215359950189254, 
         -2.8215359935918802, -2.821535992164835, -2.8215359907377904, 
         -2.8215359893107452, -2.821535987589213, -2.8215359858676803, 
         -2.821535984146148, -2.8215359824246153, -2.821535980347608, 
         -2.8215359782706, -2.8215359761935925, -2.8215359741165846, 
         -2.821535971610427, -2.8215359691042696, -2.8215359665981117, 
         -2.821535964091954, -2.8215359610676303, -2.821535958043307, 
         -2.821535955018983, -2.821535951994659, -2.8215359483445988, 
         -2.8215359446945385, -2.821535941044478, -2.8215359373944175, 
         -2.821535932988604, -2.821535928582791, -2.8215359241769775, 
         -2.821535919771164, -2.8215359144524226, -2.8215359091336807, 
         -2.8215359038149392, -2.8215358984961973, -2.821535892074482, 
         -2.8215358856527666, -2.8215358792310514, -2.821535872809336, 
         -2.8215358650548095, -2.8215358573002827, -2.8215358495457563, 
         -2.8215358417912295, -2.8215358324258704, -2.821535823060511, 
         -2.821535813695152, -2.821535804329793, -2.8215357930172185, 
         -2.8215357817046436, -2.821535770392069, -2.8215357590794947, 
         -2.821535745412613, -2.821535731745731, -2.821535718078849, 
         -2.8215357044119678, -2.8215356878979936, -2.821535671384019, 
         -2.8215356548700448, -2.8215356383560706, -2.8215356183983324, 
         -2.8215355984405943, -2.821535578482856, -2.821535558525118, 
         -2.821535534400964, -2.8215355102768096, -2.821535486152656, 
         -2.8215354620285016, -2.821535432862453, -2.8215354036964047, 
         -2.8215353745303564, -2.8215353453643077, -2.8215353100954315, 
         -2.8215352748265548, -2.8215352395576785, -2.821535204288802, 
         -2.8215351616310294, -2.8215351189732565, -2.8215350763154836, 
         -2.8215350336577107, -2.8215349820515687, -2.8215349304454267, 
         -2.8215348788392847, -2.8215348272331426, -2.8215347647870024, 
         -2.821534702340862, -2.821534639894722, -2.8215345774485816, 
         -2.821534501867117, -2.8215344262856523, -2.8215343507041877, 
         -2.821534275122723, -2.821534183619752, -2.821534092116781, 
         -2.82153400061381, -2.821533909110839, -2.821533798303082, 
         -2.8215336874953247, -2.821533576687568, -2.821533465879811, 
         -2.8215333316573186, -2.821533197434826, -2.821533063212333, 
         -2.821532928989841, -2.8215327663578798, -2.8215326037259185, 
         -2.8215324410939573, -2.8215322784619956, -2.8215320813479927, 
         -2.8215318842339894, -2.8215316871199865, -2.821531490005983, 
         -2.821531251023664, -2.8215310120413446, -2.8215307730590253, 
         -2.821530534076706, -2.821530244237757, -2.821529954398808, 
         -2.8215296645598587, -2.82152937472091, -2.8215290230822583, 
         -2.8215286714436068, -2.821528319804955, -2.8215279681663037, 
         -2.8215275413983516, -2.8215271146303995, -2.821526687862447, 
         -2.821526261094495, -2.821525742952222, -2.8215252248099487, 
         -2.821524706667675, -2.821524188525402, -2.821523559200021, 
         -2.8215229298746394, -2.8215223005492582, -2.8215216712238766, 
         -2.821520906547469, -2.8215201418710616, -2.821519377194654, 
         -2.8215186125182465, -2.8215176829872877, -2.8215167534563292, 
         -2.8215158239253704, -2.821514894394412, -2.8215137639699814, 
         -2.821512633545551, -2.82151150312112, -2.8215103726966895, 
         -2.8215089973290786, -2.8215076219614676, -2.8215062465938563, 
         -2.8215048712262454, -2.8215031970391, -2.821501522851955, 
         -2.82149984866481, -2.821498174477665, -2.8214961355311248, 
         -2.8214940965845843, -2.821492057638044, -2.8214900186915037, 
         -2.821487534224202, -2.8214850497569004, -2.8214825652895987, 
         -2.821480080822297, -2.8214770518476175, -2.8214740228729376, 
         -2.821470993898258, -2.8214679649235785, -2.8214642700248778, 
         -2.8214605751261774, -2.821456880227477, -2.8214531853287768, 
         -2.821448675459863, -2.821444165590949, -2.821439655722035, 
         -2.8214351458531217, -2.8214296379019985, -2.8214241299508753, 
         -2.821418621999752, -2.8214131140486285, -2.821406382861871, 
         -2.821399651675114, -2.821392920488356, -2.821386189301599, 
         -2.821377957790349, -2.8213697262790984, -2.8213614947678485, 
         -2.821353263256598, -2.82134319010891, -2.8213331169612217, 
         -2.8213230438135337, -2.8213129706658457, -2.8213006350686003, 
         -2.821288299471355, -2.82127596387411, -2.8212636282768644, 
         -2.8212485108968384, -2.8212333935168123, -2.8212182761367863, 
         -2.8212031587567603, -2.821184618035943, -2.8211660773151257, 
         -2.8211475365943084, -2.821128995873491, -2.8211062384471353, 
         -2.8210834810207794, -2.8210607235944236, -2.821037966168068, 
         -2.8210100099007858, -2.8209820536335037, -2.8209540973662213, 
         -2.8209261410989392, -2.8208917688358692, -2.8208573965727997, 
         -2.8208230243097296, -2.8207886520466596, -2.820746353669415, 
         -2.82070405529217, -2.8206617569149253, -2.8206194585376805, 
         -2.8205673582675685, -2.820515257997456, -2.8204631577273442, 
         -2.8204110574572323, -2.8203468225212944, -2.820282587585357, 
         -2.820218352649419, -2.8201541177134812, -2.8200748434306657, 
         -2.81999556914785, -2.8199162948650343, -2.8198370205822187, 
         -2.8197390855706255, -2.819641150559032, -2.819543215547439, 
         -2.8194452805358456, -2.8193241640009887, -2.819203047466132, 
         -2.819081930931275, -2.818960814396418, -2.8188108652219794, 
         -2.818660916047541, -2.8185109668731023, -2.8183610176986638, 
         -2.8181751621793407, -2.817989306660018, -2.817803451140695, 
         -2.8176175956213725, -2.817386966961183, -2.8171563383009937, 
         -2.8169257096408047, -2.8166950809806153, -2.816408548881366, 
         -2.8161220167821166, -2.815835484682867, -2.815548952583618, 
         -2.815192525043272, -2.814836097502926, -2.81447966996258, 
         -2.814123242422234, -2.8136793038287986, -2.8132353652353626, 
         -2.812791426641927, -2.812347488048491, -2.8117938272397, 
         -2.8112401664309097, -2.810686505622119, -2.810132844813328, 
         -2.8094414132546826, -2.8087499816960366, -2.808058550137391, 
         -2.8073671185787457, -2.8065024412495916, -2.805637763920438, 
         -2.8047730865912843, -2.8039084092621307, -2.8028255479998405, 
         -2.8017426867375503, -2.80065982547526, -2.79957696421297, 
         -2.7982189035404956, -2.7968608428680213, -2.7955027821955474, 
         -2.794144721523073, -2.792439008274788, -2.790733295026503, 
         -2.789027581778218, -2.7873218685299332, -2.7851762861555778, 
         -2.7830307037812223, -2.780885121406867, -2.7787395390325114, 
         -2.776036531731887, -2.7733335244312625, -2.770630517130638, 
         -2.7679275098300136, -2.764516986012208, -2.7611064621944026, 
         -2.757695938376597, -2.754285414558791, -2.7499754475081395, 
         -2.7456654804574883, -2.7413555134068366, -2.737045546356185, 
         -2.7315903446448613, -2.7261351429335376, -2.7206799412222136, 
         -2.71522473951089, -2.7083090706562003, -2.7013934018015107, 
         -2.694477732946821, -2.6875620640921314, -2.6787810700836507, 
         -2.6700000760751696, -2.661219082066689, -2.652438088058208, 
         -2.6412711217727605, -2.630104155487313, -2.618937189201865, 
         -2.6077702229164177, -2.5935470151184075, -2.579323807320397, 
         -2.5651005995223866, -2.550877391724376, -2.5327352575728703, 
         -2.5145931234213643, -2.4964509892698588, -2.4783088551183527, 
         -2.4551485976506515, -2.4319883401829503, -2.4088280827152486, 
         -2.3856678252475474, -2.3561974659473655, -2.326727106647184, 
         -2.297256747347002, -2.26778638804682, -2.232374919295588, 
         -2.196963450544356, -2.161551981793124, -2.1163945270918867, 
         -2.071237072390649, -2.026079617689412, -1.9809221629881748, 
         -1.9265709358774212, -1.8722197087666677, -1.8178684816559143, 
         -1.7484323726898843, -1.678996263723854, -1.609560154757824, 
         -1.540124045791794, -1.4564910395639807, -1.3728580333361675, 
         -1.2892250271083545, -1.1828450760510472, -1.0764651249937398, 
         -0.9700851739364325, -0.8637052228791251, -0.7573252718218177, 
         -0.6509453207645105, -0.5445653697072032, -0.4578424462888541, 
         -0.3711195228705051, -0.28439659945215606, -0.19767367603380706, 
         -0.11095075261545803, -0.08573467247558121, -0.060518592335704384, 
         -0.03530251219582756, -0.010086432055950731, -0.007564824041963049, 
         -0.005043216027975366, -0.002521608013987683, 
         -0.0012608040069938414, 0., 0.0012624195982784533, 
         0.0025248391965569067, 0.005049678393113813, 0.00757451758967072, 
         0.010099356786227627, 0.03534774875179669, 0.060596140717365764, 
         0.08584453268293483, 0.1110929246485039, 0.20340087437875976, 
         0.29570882410901556, 0.3880167738392714, 0.48032472356952727, 
         0.5726326732997831, 0.7108150493211634, 0.8489974253425436, 
         0.9871798013639238, 1.125362177385304, 1.2635445534066843, 
         1.4017269294280645, 1.5965096181117162, 1.791292306795368, 
         1.9860749954790198, 2.1808576841626715, 2.3756403728463233, 
         2.5704230615299752, 2.818717477054328, 3.0670118925786807, 
         3.3153063081030334, 3.5636007236273857, 3.8118951391517384, 
         4.060189554676091, 4.3794075994587285, 4.698625644241365, 
         5.017843689024002, 5.33706173380664, 5.656279778589276, 
         5.975497823371914, 6.386972718151893, 6.798447612931872, 
         7.209922507711851, 7.621397402491829, 8.032872297271808, 
         8.444347192051788, 8.97686318993129, 9.509379187810794, 
         10.041895185690297, 10.5744111835698, 11.106927181449304, 
         11.639443179328808, 12.331892516381835, 13.024341853434864, 
         13.716791190487891, 14.40924052754092, 15.101689864593947, 
         15.794139201646976, 16.69958625618837, 17.60503331072976, 
         18.51048036527115, 19.255240182635575, 20.}}, 
        {Developer`PackedArrayForm, {0, 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 
         22, 24, 26, 28, 30, 32, 34, 36, 38, 40, 42, 44, 46, 48, 50, 52, 54, 
         56, 58, 60, 62, 64, 66, 68, 70, 72, 74, 76, 78, 80, 82, 84, 86, 88, 
         90, 92, 94, 96, 98, 100, 102, 104, 106, 108, 110, 112, 114, 116, 
         118, 120, 122, 124, 126, 128, 130, 132, 134, 136, 138, 140, 142, 
         144, 146, 148, 150, 152, 154, 156, 158, 160, 162, 164, 166, 168, 
         170, 172, 174, 176, 178, 180, 182, 184, 186, 188, 190, 192, 194, 
         196, 198, 200, 202, 204, 206, 208, 210, 212, 214, 216, 218, 220, 
         222, 224, 226, 228, 230, 232, 234, 236, 238, 240, 242, 244, 246, 
         248, 250, 252, 254, 256, 258, 260, 262, 264, 266, 268, 270, 272, 
         274, 276, 278, 280, 282, 284, 286, 288, 290, 292, 294, 296, 298, 
         300, 302, 304, 306, 308, 310, 312, 314, 316, 318, 320, 322, 324, 
         326, 328, 330, 332, 334, 336, 338, 340, 342, 344, 346, 348, 350, 
         352, 354, 356, 358, 360, 362, 364, 366, 368, 370, 372, 374, 376, 
         378, 380, 382, 384, 386, 388, 390, 392, 394, 396, 398, 400, 402, 
         404, 406, 408, 410, 412, 414, 416, 418, 420, 422, 424, 426, 428, 
         430, 432, 434, 436, 438, 440, 442, 444, 446, 448, 450, 452, 454, 
         456, 458, 460, 462, 464, 466, 468, 470, 472, 474, 476, 478, 480, 
         482, 484, 486, 488, 490, 492, 494, 496, 498, 500, 502, 504, 506, 
         508, 510, 512, 514, 516, 518, 520, 522, 524, 526, 528, 530, 532, 
         534, 536, 538, 540, 542, 544, 546, 548, 550, 552, 554, 556, 558, 
         560, 562, 564, 566, 568, 570, 572, 574, 576, 578, 580, 582, 584, 
         586, 588, 590, 592, 594, 596, 598, 600, 602, 604, 606, 608, 610, 
         612, 614, 616, 618, 620, 622, 624, 626, 628, 630, 632, 634, 636, 
         638, 640, 642, 644, 646, 648, 650, 652, 654, 656, 658, 660, 662, 
         664, 666, 668, 670, 672, 674, 676, 678, 680, 682, 684, 686, 688, 
         690, 692, 694, 696, 698, 700, 702, 704, 706, 708, 710, 712, 714, 
         716, 718, 720, 722, 724, 726, 728, 730, 732, 734, 736, 738, 740, 
         742, 744, 746, 748, 750, 752, 754, 756, 758, 760, 762, 764, 766, 
         768, 770, 772, 774, 776, 778, 780, 782, 784, 786, 788, 790, 792, 
         794, 796, 798, 800, 802, 804, 806, 808, 810, 812, 814, 816, 818, 
         820, 822, 824, 826, 828, 830, 832, 834, 836, 838, 840, 842, 844, 
         846, 848, 850, 852, 854, 856, 858, 860, 862, 864, 866, 868, 870, 
         872, 874, 876, 878, 880, 882, 884, 886, 888, 890, 892, 894, 896, 
         898, 900, 902, 904, 906, 908, 910, 912, 914, 916, 918, 920, 922, 
         924, 926, 928, 930, 932, 934, 936, 938, 940, 942, 944, 946, 948, 
         950, 952, 954, 956, 958, 960, 962, 964, 966, 968, 970, 972, 974, 
         976, 978, 980, 982, 984, 986, 988, 990, 992, 994, 996, 998, 1000, 
         1002, 1004, 1006, 1008, 1010, 1012, 1014, 1016, 1018, 1020, 1022, 
         1024, 1026, 1028, 1030, 1032, 1034, 1036, 1038, 1040, 1042, 1044, 
         1046, 1048, 1050, 1052, 1054, 1056, 1058, 1060, 1062, 1064, 1066, 
         1068, 1070, 1072, 1074, 1076, 1078, 1080, 1082, 1084, 1086, 1088, 
         1090, 1092, 1094, 1096, 1098, 1100, 1102, 1104, 1106, 1108, 1110, 
         1112, 1114, 1116, 1118, 1120, 1122, 1124, 1126, 1128, 1130, 1132, 
         1134, 1136, 1138, 1140, 1142, 1144, 1146, 1148, 1150, 1152, 1154, 
         1156, 1158, 1160, 1162, 1164, 1166, 1168, 1170, 1172, 1174, 1176, 
         1178, 1180, 1182, 1184, 1186, 1188, 1190, 1192, 1194, 1196, 1198, 
         1200, 1202, 1204, 1206, 1208, 1210, 1212, 1214, 1216, 1218, 1220, 
         1222, 1224, 1226, 1228, 1230, 1232, 1234, 1236, 1238, 1240, 1242, 
         1244, 1246, 1248, 1250, 1252, 1254, 1256, 1258, 1260, 1262, 1264, 
         1266, 1268, 1270, 1272, 1274, 1276, 1278, 1280, 1282, 1284, 1286, 
         1288, 1290, 1292, 1294, 1296, 1298, 1300, 1302, 1304, 1306, 1308, 
         1310, 1312}, {282.45145147841544, -8.496632240158703*^13, 
         278.96300943521095, -7.984817053515383*^13, 275.67984701468936, 
         -7.525943059489636*^13, 272.58119524001035, -7.112448119860981*^13, 
         269.6491911635791, -6.7381233805337914*^13, 266.3187217243788, 
         -6.332222092696874*^13, 263.1842425781003, -5.968307859448284*^13, 
         260.2259239425064, -5.64038304492598*^13, 257.4267107374313, 
         -5.343523354202546*^13, 254.24705975950712, -5.0216203566049805*^13, 
         251.25453043237204, -4.733016659644955*^13, 248.4301896349061, 
         -4.472955470413336*^13, 245.75775356833134, -4.237531301992638*^13, 
         242.7220909006826, -3.982245408906089*^13, 239.86508287750507, 
         -3.7533685631113625*^13, 237.16865243153657, -3.54712811043098*^13, 
         234.61725211622917, -3.3604265974835297*^13, 231.71905692759327, 
         -3.1579729134773387*^13, 228.9914328238352, -2.9764636416692508*^13, 
         226.41712002329564, -2.8129066670267426*^13, 223.98127408435602, 
         -2.664845399115006*^13, 221.21432085379635, -2.5042912682107125*^13, 
         218.61022123002152, -2.360347584912672*^13, 216.15249570442202, 
         -2.2306412929214285*^13, 213.8269710032111, -2.1132240778194824*^13, 
         211.18531614759365, -1.9858989644714996*^13, 208.69914692322274, 
         -1.8717470759313297*^13, 206.35272875819769, 
         -1.7688863403840617*^13, 204.1325291512936, -1.6757715618787691*^13, 
         201.61049825596788, -1.5747990188849662*^13, 199.23691869090206, 
         -1.4842736868916844*^13, 196.99676708484532, 
         -1.4027028998409172*^13, 194.87712269038067, 
         -1.3288611777440186*^13, 192.4692983202082, -1.2487876444553482*^13, 
         190.2032095425581, -1.1769993685627346*^13, 188.06451197859593, 
         -1.1123125254648959*^13, 186.04086893014826, 
         -1.0537552031643387*^13, 183.74207899079676, -9.902556017171355*^12, 
         181.5786130426176, -9.333266587657334*^12, 179.5367749495796, 
         -8.820294940296998*^12, 177.60478561285186, -8.355933005496946*^12, 
         175.41009224089427, -7.852375422569231*^12, 173.34460162153334, 
         -7.40092637991647*^12, 171.39523650450738, -6.994139738537267*^12, 
         169.55075013828161, -6.625902388884094*^12, 167.45543909284984, 
         -6.226580394341158*^12, 165.48348677444488, -5.868582737650249*^12, 
         163.62240679104247, -5.546004009102594*^12, 161.86146062995618, 
         -5.253996053701447*^12, 159.86103116587836, -4.937336145313517*^12, 
         157.97838106196681, -4.653447995053743*^12, 156.20158802620142, 
         -4.3976481430403955*^12, 154.52039876810957, -4.166091658960419*^12, 
         152.61055396680382, -3.914984849628354*^12, 150.8131618318194, 
         -3.6898670822192646*^12, 149.1168386247149, -3.4870240193414946*^12, 
         147.51179430295855, -3.3034063761296616*^12, 145.68843184229303, 
         -3.104284376426687*^12, 143.97243658415536, -2.925772345445131*^12, 
         142.35293895419565, -2.7649245489296826*^12, 140.82059117871432, 
         -2.6193227159522925*^12, 139.07979451999196, -2.461425261949376*^12, 
         137.44150990534197, -2.319871898707579*^12, 135.89535863955822, 
         -2.1923264621709722*^12, 134.4324151887569, -2.076871126257755*^12, 
         132.7704451630649, -1.951664905838787*^12, 131.2063518983672, 
         -1.8394197416662747*^12, 129.73022534817102, 
         -1.7382829104906682*^12, 128.33354309645654, 
         -1.6467334772663308*^12, 126.74682987215097, -1.547451289507677*^12, 
         125.25356804351368, -1.458447277557881*^12, 123.84429498909685, 
         -1.378252125834234*^12, 122.51087315478034, -1.3056596625659094*^12, 
         120.99600856879508, -1.2269349338628435*^12, 119.57037042421655, 
         -1.1563607664272932*^12, 118.22492326214446, -1.092771916437992*^12, 
         116.9518969611649, -1.0352118407612426*^12, 115.50562719735656, 
         -9.727887680102795*^11, 114.14455025561583, -9.168289998062477*^11, 
         112.86003849473244, -8.664084433413071*^11, 111.64467258613236, 
         -8.207685267779744*^11, 110.26389118555379, -7.712721998741117*^11, 
         108.96445165749147, -7.269011476508004*^11, 107.7381157096539, 
         -6.869225265272455*^11, 106.57779892002407, -6.507347674661151*^11, 
         105.25954010722735, -6.114888607227428*^11, 104.01894661450132, 
         -5.763071521465916*^11, 102.84815186610768, -5.446084709783214*^11, 
         101.7403911797555, -5.159157597745363*^11, 100.48182349196166, 
         -4.8479794963011176*^11, 99.29741107082432, -4.5690289664826306*^11, 
         98.17964222332564, -4.317697155491338*^11, 97.12205752694437, 
         -4.090200925626994*^11, 95.92047773121197, -3.8434739528882587*^11, 
         94.78970210738028, -3.6223019296444653*^11, 93.72255777368557, 
         -3.423029931710611*^11, 92.71287674318529, -3.2426579006406665*^11, 
         91.56570402931631, -3.0470362840115533*^11, 90.48613615234994, 
         -2.871678283983928*^11, 89.46732369815263, -2.7136854880749347*^11, 
         88.50337691723819, -2.570678848369745*^11, 87.40814735365157, 
         -2.415579734169768*^11, 86.37746817883666, -2.276548297933897*^11, 
         85.40479879754731, -2.1512860734925452*^11, 84.48451509734821, 
         -2.0379064932356854*^11, 83.43887633802277, -1.914938047045229*^11, 
         82.45487184456161, -1.8047100057349698*^11, 81.52625585527625, 
         -1.7053996828775305*^11, 80.6476578650252, -1.6155110596240787*^11, 
         79.64936409604249, -1.5180189285317856*^11, 78.7099205307601, 
         -1.4306287890916068*^11, 77.82336288908529, -1.3518952630131267*^11, 
         76.9845627882181, -1.2806320434173161*^11, 76.0314699032179, 
         -1.2033396864406992*^11, 75.1345692397172, -1.1340571750614922*^11, 
         74.28816525194144, -1.071638345673508*^11, 73.48736071461185, 
         -1.015142496208284*^11, 72.57742170900265, -9.538659573948763*^10, 
         71.72113731269627, -8.989403047949861*^10, 70.91306854426298, 
         -8.494566574483876*^10, 70.14853886770534, -8.046690786464983*^10, 
         69.27979944099502, -7.560907969109283*^10, 68.46229193004514, 
         -7.125479529877806*^10, 67.6908222989539, -6.7331983367694534*^10, 
         66.96092470685424, -6.3781495028903496*^10, 66.1315190643491, 
         -5.993043880948725*^10, 65.35103235824036, -5.647862118063323*^10, 
         64.61450440553344, -5.3368892595557884*^10, 63.917670518959646, 
         -5.0554351716185005*^10, 63.12581736383351, -4.750149183681879*^10, 
         62.38067491110794, -4.476516554482143*^10, 61.67750624044534, 
         -4.230005620027971*^10, 61.01223870876794, -4.0068974497157875*^10, 
         60.25623741527385, -3.764893373130982*^10, 59.544838591777285, 
         -3.547984961938397*^10, 58.873518469964914, -3.3525788874265392*^10, 
         58.23838775409589, -3.175726353905014*^10, 57.51661471422692, 
         -2.9838918337434914*^10, 56.83743138459693, -2.8119534551299973*^10, 
         56.19651749618162, -2.657061970020342*^10, 55.59015879759351, 
         -2.516879269375137*^10, 54.90106393329078, -2.3648180209043*^10, 
         54.252637168251546, -2.228530336110931*^10, 53.64075251716802, 
         -2.1057570011055195*^10, 53.06186284605754, -1.994644151361085*^10, 
         52.403966279213165, -1.8741135075470245*^10, 51.784903221395325, 
         -1.7660877907836666*^10, 51.20073317280455, -1.6687756614980076*^10, 
         50.648068548919476, -1.580707232329336*^10, 50.01995742229803, 
         -1.4851721404610598*^10, 49.42892829409609, -1.3995504438437057*^10, 
         48.87121775027878, -1.3224218046377867*^10, 48.343590530590525, 
         -1.2526208024864435*^10, 47.743915972659224, 
         -1.1769001647734568*^10, 47.179651219539025, 
         -1.1090382707099224*^10, 46.6472019236353, -1.0479089651334984*^10, 
         46.14347825098892, -9.925882452883198*^9, 45.570952478244706, 
         -9.32574487414211*^9, 45.03224004145994, -8.787905102775045*^9, 
         44.52390800334827, -8.303434252350867*^9, 44.043005370685655, 
         -7.86500651242382*^9, 43.496398921218834, -7.389373594107122*^9, 
         42.982081634046516, -6.963124055992195*^9, 42.496774672804186, 
         -6.579178627749503*^9, 42.03765959770259, -6.231729271569706*^9, 
         41.51579868995468, -5.854785319362643*^9, 41.02477179176796, 
         -5.516986405826074*^9, 40.56144718935548, -5.21271983593024*^9, 
         40.12313299381593, -4.937380986196355*^9, 39.62489700515291, 
         -4.63866072310504*^9, 39.15610576825237, -4.370968898730579*^9, 
         38.713768029613924, -4.129855659235696*^9, 38.29531272053774, 
         -3.9116705302238693*^9, 37.819631779714484, -3.6749506058327446*^9, 
         37.37206924337016, -3.462824717149132*^9, 36.94976795773816, 
         -3.271765007098938*^9, 36.55027220318086, -3.098877502336554*^9, 
         36.096124892118176, -2.9112975884192266*^9, 35.668829699480206, 
         -2.743210990366929*^9, 35.26565749875751, -2.5918208740251737*^9, 
         34.8842626960309, -2.454833049809865*^9, 34.45067385557915, 
         -2.3061989827376313*^9, 34.04272818846543, -2.1730148636485972*^9, 
         33.657818797992256, -2.053063453672143*^9, 33.29370522918982, 
         -1.9445260860292857*^9, 32.87974386414816, -1.826757358532486*^9, 
         32.49027147125756, -1.7212334680307164*^9, 32.12279784876258, 
         -1.6261967021623278*^9, 31.7751829197343, -1.5402053921628172*^9, 
         31.379960198956077, -1.4468971016391506*^9, 31.008124513564262, 
         -1.3632930828843517*^9, 30.657297072562642, -1.2879999475382733*^9, 
         30.325433631789178, -1.2198749072218823*^9, 29.94810097868041, 
         -1.14595045659146*^9, 29.59310332142335, -1.0797165326714578*^9, 
         29.25816823493227, -1.0200685862602824*^9, 28.941342968407263, 
         -9.661007731083119*^8, 28.58109023798873, -9.075365283790141*^8, 
         28.24216810111436, -8.550666928504332*^8, 27.922405682257118, 
         -8.078157065221021*^8, 27.619937581121793, -7.650655634686381*^8, 
         27.275991319652572, -7.186725639503722*^8, 26.95241672897862, 
         -6.771089145298691*^8, 26.647139884930333, -6.39680591775751*^8, 
         26.35837878250757, -6.058184955835439*^8, 26.03000056586923, 
         -5.690693989324061*^8, 25.72107851687321, -5.3614702014687306*^8, 
         25.42963127277802, -5.0650124976747525*^8, 25.1539564478195, 
         -4.796810273204707*^8, 24.840441295229795, -4.5057291889470553*^8, 
         24.545508260034225, -4.244969134795739*^8, 24.267264349832104, 
         -4.0101698726337254*^8, 24.00408319307536, -3.79775648809582*^8, 
         23.704758052535, -3.567212853670207*^8, 23.423180554413257, 
         -3.360692938765378*^8, 23.157542075400375, -3.1747408786190176*^8, 
         22.90628881643836, -3.006523451655933*^8, 22.62051111887847, 
         -2.8239394178219265*^8, 22.351684371414954, -2.6603886941461608*^8, 
         22.098080499821705, -2.5131323060743785*^8, 21.858214991705605, 
         -2.379924931315438*^8, 21.585371270551942, -2.235333823520545*^8, 
         21.328717878375492, -2.1058214580758926*^8, 21.086603643062503, 
         -1.9892171178353673*^8, 20.857610201901384, -1.883741461446733*^8, 
         20.597114775176017, -1.7692460592428452*^8, 20.352083493551813, 
         -1.6666958048058486*^8, 20.120938605257606, -1.574370408104767*^8, 
         19.902324902382826, -1.4908598615904802*^8, 19.653618614453794, 
         -1.4002030030524498*^8, 19.41968316501745, -1.3190083501119374*^8, 
         19.1990108985967, -1.2459126308705689*^8, 18.99030690287018, 
         -1.1797985595032515*^8, 18.75285592324293, -1.1080227907349429*^8, 
         18.529513863418273, -1.0437419326807511*^8, 18.318839990745598, 
         -9.858756842812167*^7, 18.119596958807715, -9.335387528350407*^7, 
         17.89289163524063, -8.767165656590629*^7, 17.679663278791043, 
         -8.258306328202239*^7, 17.478535049929942, -7.800249268053776*^7, 
         17.288324562121332, -7.385980113197328*^7, 17.071878325740027, 
         -6.936180365352508*^7, 16.868305712256863, -6.5333948284909815*^7, 
         16.676290882814, -6.170840682599598*^7, 16.494703922797548, 
         -5.842960583399161*^7, 16.28805224276569, -5.486936139106817*^7, 
         16.093698153800183, -5.168143579281365*^7, 15.910384056304498, 
         -4.881208732578283*^7, 15.737030132346701, -4.621728859406152*^7, 
         15.53972951787781, -4.339956835218587*^7, 15.354176537629446, 
         -4.087667604524172*^7, 15.179169194964794, -3.860603317703274*^7, 
         15.01367550100995, -3.6552761538195014*^7, 14.825302548534843, 
         -3.432293345642799*^7, 14.648152167041712, -3.2326552456840012*^7, 
         14.481075445995913, -3.0529886265628584*^7, 14.32308606070176, 
         -2.8905306234053876*^7, 14.143236544254185, -2.7140901312792927*^7, 
         13.974108301255473, -2.556132656234929*^7, 13.81460310445333, 
         -2.4139859534715977*^7, 13.663778226544254, -2.285461565416314*^7, 
         13.492066229280013, -2.1458644857330173*^7, 13.33059689679583, 
         -2.0208999617367394*^7, 13.178320391177294, -1.908451187858908*^7, 
         13.034335609380408, -1.806784805917523*^7, 12.870392694420506, 
         -1.6963506743983544*^7, 12.716235496342783, -1.5974997452051753*^7, 
         12.570860376573126, -1.5085553061585872*^7, 12.43340597260846, 
         -1.4281445524131233*^7, 12.27688039144668, -1.340791750743907*^7, 
         12.129704258483697, -1.2626074087782426*^7, 11.990918043730197, 
         -1.1922634975717256*^7, 11.859698326872657, -1.1286728454328835*^7, 
         11.710254263608945, -1.0595861536403865*^7, 11.569743121949132, 
         -9.977558659001913*^6, 11.437247484483422, -9.421301980769807*^6, 
         11.311980156240008, -8.918482460724983*^6, 11.169297006029767, 
         -8.372154548040991*^6, 11.035149098208725, -7.8832517454277165*^6, 
         10.908659222403204, -7.443444642623127*^6, 10.78907476995016, 
         -7.045916493663923*^6, 10.65284645011622, -6.613948973736972*^6, 
         10.524773686625114, -6.227421146426228*^6, 10.404017656963704, 
         -5.879737190347025*^6, 10.289858774038132, -5.565499968481949*^6, 
         10.159793066360974, -5.22400472504169*^6, 10.037520406599338, 
         -4.918461217770462*^6, 9.922238623384231, -4.643646821397413*^6, 
         9.813259657382742, -4.3952887415152965*^6, 9.689077580183806, 
         -4.1253587687536557*^6, 9.572342441451285, -3.883870325267463*^6, 
         9.462287062972026, -3.6666884401478665*^6, 9.358253487087833, 
         -3.4704307450841046*^6, 9.239688695676486, -3.2571030025001387*^6, 
         9.128240388862674, -3.066272651151012*^6, 9.023174798765872, 
         -2.8946661496078107*^6, 8.923862707973715, -2.7396062524877912*^6, 
         8.810660922285399, -2.5710402110348814*^6, 8.70426011314745, 
         -2.420267405943712*^6, 8.603958411953428, -2.2846965049765464*^6, 
         8.5091540418403, -2.1622086017251657*^6, 8.401072499937898, 
         -2.0290356687448798*^6, 8.299490694721507, -1.9099335300358443*^6, 
         8.203737214317933, -1.802851088543098*^6, 8.113236481656545, 
         -1.7061115124201663*^6, 8.01004341793803, -1.6009198259187012*^6, 
         7.913062472251744, -1.5068536598904568*^6, 7.8216513123345655, 
         -1.4222896908945886*^6, 7.7352593764243665, -1.3459010772273245*^6, 
         7.636733523457512, -1.2628275328141132*^6, 7.544145173382063, 
         -1.1885494202140977*^6, 7.456879758880363, -1.1217821113395055*^6, 
         7.374410602742649, -1.0614758264793735*^6, 7.280340715612538, 
         -995883.0432396638, 7.191946129998724, -937242.520008805, 
         7.108638788497774, -884537.7719035208, 7.029914818993794, 
         -836938.3284868392, 6.940099221149689, -785158.8906931984, 
         6.855708574150302, -738873.8951987358, 6.776180132411311, 
         -697279.1252361523, 6.701031798429835, -659717.6880538708, 
         6.615277948053653, -618851.669079048, 6.534710010969046, 
         -582327.2180855599, 6.458789409676827, -549508.1383717253, 
         6.387054837568974, -519874.90515720827, 6.305178913530764, 
         -487629.58771182975, 6.228260665088349, -458814.2489359635, 
         6.155784590996484, -432925.6744297614, 6.087309236472004, 
         -409553.0394027399, 6.009135742991215, -384116.0502842458, 
         5.9357019972287235, -361388.343949397, 5.866514531918033, 
         -340971.9459500476, 5.8011508476603675, -322542.03646914987, 
         5.726512236789491, -302480.93960203545, 5.6564052877121815, 
         -284559.4371961015, 5.590357572180719, -268462.8306775181, 
         5.527964690442346, -253934.34036358277, 5.456701001597706, 
         -238117.17838124087, 5.389770283877482, -223989.3628557012, 
         5.326720198269055, -211302.08332706813, 5.267163627789569, 
         -199852.38109921684, 5.19912214351126, -187384.80165123136, 
         5.1352239084599605, -176250.78002494163, 5.0750357662084085, 
         -166253.63023388275, 5.018187102772995, -157232.96128977204, 
         4.953222019998764, -147408.47214057553, 4.892219026140667, 
         -138636.44487595296, 4.834763281891302, -130761.42417023177, 
         4.780499931924022, -123656.68166750633, 4.718472048062758, 
         -115917.29782312008, 4.660233265633988, -109008.33365251542, 
         4.605386236307643, -102806.9430233325, 4.553591152899075, 
         -97213.01414481852, 4.494367565995175, -91118.12993093781, 
         4.438767894704514, -85678.30517649109, 4.386411493077429, 
         -80796.4805782954, 4.336972923854718, -76393.58643378825, 
         4.280426746236961, -71595.35669666606, 4.227346745714622, 
         -67313.72236941963, 4.1773682259705724, -63472.01290739979, 
         4.130179472293515, -60007.791016985255, 4.076189557054684, 
         -56231.66299925661, 4.025515189375428, -52862.823104605006, 
         3.977806904052404, -49840.71893961557, 3.9327660909901665, 
         -47116.05781094189, 3.8812167707068896, -44145.378307335675, 
         3.8328391544312828, -41495.71632947774, 3.7872983222374077, 
         -39119.259434718406, 3.744308178828494, -36977.101251114655, 
         3.695089015968304, -34640.94779717492, 3.6489041911804865, 
         -32557.73910790324, 3.6054326751858556, -30689.730773875774, 
         3.564400324513699, -29006.221720297526, 3.5174058729622732, 
         -27169.78405829915, 3.473314576771692, -25532.59096914116, 
         3.4318186724744826, -24064.853076419935, 3.392655431084794, 
         -22742.350288306625, 3.347785008298948, -21299.328896386458, 
         3.305692460338878, -20013.20094542437, 3.2660826931632414, 
         -18860.461341105383, 3.228703879403107, -17822.007953766013, 
         3.185861348657633, -16688.60768958493, 3.1456770460816346, 
         -15678.70940106607, 3.1078679778380423, -14773.77141765261, 
         3.0721927286739423, -13958.732325139079, 3.0312862909490415, 
         -13068.917325040227, 2.992923812497057, -12276.284566093258, 
         2.956833856399434, -11566.21173154146, 2.9227849523291125, 
         -10926.826885226888, 2.8837269473751728, -10228.572586549442, 
         2.8471037660733507, -9606.761083772135, 2.8126550099016443, 
         -9049.864273103565, 2.78015870757067, -8548.527176609156, 
         2.742865423700157, -8000.860521951379, 2.707902727771148, 
         -7513.2995196587235, 2.675020764780304, -7076.758278640269, 
         2.644006636926133, -6683.867559497119, 2.608398129121896, 
         -6254.530122113322, 2.5750206507183075, -5872.433273565985, 
         2.543634417927447, -5530.418243770949, 2.514035200303934, 
         -5222.682011058031, 2.4800351162390677, -4886.285346230245, 
         2.4481709676189713, -4587.0018757218095, 2.418212591121064, 
         -4319.193100709329, 2.3899640360630823, -4078.290980532077, 
         2.357499449638371, -3814.860858591648, 2.3270799664124446, 
         -3580.574517634345, 2.298484613355232, -3370.9924441794865, 
         2.2715253496532513, -3182.520107184917, 2.2405266016867422, 
         -2976.3479232030427, 2.211486192792782, -2793.050634828395, 
         2.1841919297058072, -2629.134053549956, 2.158463328482664, 
         -2481.770807477907, 2.128863874193799, -2320.507580718693, 
         2.10113987826235, -2177.1899882255184, 2.0750875354130565, 
         -2049.0689444788163, 2.050533581703233, -1933.9214283173947, 
         2.0222698446524383, -1807.8633446491024, 1.9958023924446302, 
         -1695.8765998058263, 1.9709354339228378, -1595.79920886665, 
         1.9475026036700793, -1505.8840813885163, 1.9205138358191, 
         -1407.4092402119422, 1.8952457184228424, -1319.9619493078196, 
         1.871510117658496, -1241.8426713652652, 1.84914725985359, 
         -1171.6789798995965, 1.8233754074500974, -1094.8034750333957, 
         1.799251949955621, -1026.565298749979, 1.7765960704066655, 
         -965.6290456562226, 1.7552542941166835, -910.9171219354831, 
         1.73064386909915, -850.945280726778, 1.7076128094696692, 
         -797.7346741956919, 1.6859872902136082, -750.2365081496821, 
         1.6656198562522477, -707.6051078147693, 1.6421178128970393, 
         -660.8540053225382, 1.6201291857747295, -619.3923199360073, 
         1.5994868317605395, -582.3967630859553, 1.580049048771157, 
         -549.2040141047489, 1.5576046653291904, -512.7865623626776, 
         1.536610690540817, -480.504473133061, 1.516906367280954, 
         -451.7118125026158, 1.4983554920260485, -425.8885788150615, 
         1.4769202571039546, -397.5427969159645, 1.4568752326394168, 
         -372.42797615192666, 1.4380657651316902, -350.03760726310463, 
         1.4203609067935514, -329.96427644418026, 1.3998884102665945, 
         -307.9189591642728, 1.380748609534747, -288.3962533618845, 
         1.3627926852339154, -270.99920639869845, 1.3458947135603503, 
         -255.4088077999752, 1.3263405418356629, -238.27787383673657, 
         1.3080641150326353, -223.11508591186882, 1.2909221907220416, 
         -209.60953717686454, 1.2747936478812383, -197.5115955983196, 
         1.2561152833517266, -184.21102615263024, 1.2386621628004053, 
         -172.44486233564143, 1.2222963752388054, -161.9697436576409, 
         1.2069013912733844, -152.59045679398756, 1.1890581158689668, 
         -142.27302212603482, 1.172389925251829, -133.15091339231412, 
         1.1567640055317638, -125.03376627369344, 1.142068217361302, 
         -117.76902816736762, 1.1250210201301047, -109.77301568949625, 
         1.1091009875673876, -102.707431773902, 1.0941801791544385, 
         -96.42346924404701, 1.0801506531114429, -90.80199903675712, 
         1.0638621418631256, -84.61095843406247, 1.0486550168844029, 
         -79.14354923492874, 1.034405997399879, -74.28353243808365, 
         1.0210111553749437, -69.93794926634563, 1.0054454724887933, 
         -65.1491029093061, 0.9909174470192612, -60.92256757355095, 
         0.9773082539043854, -57.16761218376784, 0.96451780325372, 
         -53.81175576315927, 0.9496405458756196, -50.11122190245086, 
         0.9357591794795131, -46.847251948610925, 0.9227591397999767, 
         -43.949080214293375, 0.910544007288971, -41.36023838284169, 
         0.8963221522940956, -38.503608355317034, 0.8830563020702726, 
         -35.985602943355346, 0.8706359668724631, -33.7510731017746, 
         0.8589682370811885, -31.75606275408284, 0.8453700713686039, 
         -29.553183886282365, 0.8326898270832036, -27.612712645331257, 
         0.8208209109036132, -25.89170647577992, 0.8096737697104921, 
         -24.355979055931943, 0.7966688266438022, -22.659039675034034, 
         0.7845454519325713, -21.16524704771487, 0.7732007783090914, 
         -19.841192562643425, 0.7625484623155333, -18.660312345400722, 
         0.7501074654446555, -17.3545181549919, 0.7385133462369411, 
         -16.2058354033467, 0.7276668003884688, -15.188296252967318, 
         0.7174845534631473, -14.281279096079162, 0.7055793690651228, 
         -13.27756209230903, 0.6944879707819858, -12.395228724720482, 
         0.6841144610205052, -11.614113612335407, 0.6743784995367187, 
         -10.918223657302281, 0.6629821000219543, -10.14754439751313, 
         0.6523679356043748, -9.470547550014798, 0.6424433655418428, 
         -8.87158877236088, 0.633130854373885, -8.338276983569422, 
         0.6222172952436301, -7.747175200740486, 0.6120559066960103, 
         -7.22829954250543, 0.6025571609270961, -6.769525166485725, 
         0.5936462028795294, -6.361261582615215, 0.5831906166957457, 
         -5.908384252409152, 0.5734585735883341, -5.511131099725472, 
         0.564363520271078, -5.160114591749533, 0.5558331623380495, 
         -4.847919626750362, 0.5458117740910314, -4.50131389351275, 
         0.5364866933598208, -4.19750021088834, 0.5277742079815028, 
         -3.929217963398089, 0.5196044686684915, -3.6907401210857187, 
         0.5099946379888483, -3.4257443300681083, 0.5010552303832478, 
         -3.193633169424348, 0.4927052459608637, -2.988797014533809, 
         0.48487716881372217, -2.806817402513201, 0.47565746561387834, 
         -2.6044195142629425, 0.4670836152003931, -2.427265557863772, 
         0.4590772051428713, -2.271026144188963, 0.45157294452719016, 
         -2.132296106452402, 0.44272326574478604, -1.97785707299817, 
         0.43449614983380636, -1.8427761672619203, 0.4268156505253795, 
         -1.7237156773755917, 0.41961859640684007, -1.6180542216364873, 
         0.41112033233491213, -1.5003161284584323, 0.4032225898078003, 
         -1.397407493776708, 0.3958517703847881, -1.3067579158961402, 
         0.38894671908202744, -1.2263515822385669, 0.38078297792036697, 
         -1.136667801379758, 0.37319893377185276, -1.0583331350350518, 
         0.36612322010392045, -0.9893702127204655, 0.3594965972374029, 
         -0.9282305520113068, 0.35165249575013213, -0.8599688343875249, 
         0.34436844839583963, -0.8003847036345119, 0.3375752065388094, 
         -0.7479585915367347, 0.3312153462095447, -0.7015022347839971, 
         0.32367837508053565, -0.6495824131098095, 0.3166829524683756, 
         -0.6042918925544586, 0.31016183528432967, -0.5644640972586279, 
         0.3040593170143305, -0.5291880880692875, 0.2968198159553969, 
         -0.48972459000503404, 0.2901044274199774, -0.455321664443215, 
         0.28384780435925533, -0.42508460496125655, 0.27799586177884056, 
         -0.39831573164444173, 0.271047880469677, -0.36834244186400394, 
         0.26460749500287034, -0.3422291097846831, 0.2586111583183354, 
         -0.31929028385893876, 0.2530063216949103, -0.29899216338064993, 
         0.24635110985804165, -0.27625715622492447, 0.24018726869722812, 
         -0.25646133651827524, 0.23445305007623368, -0.23908084503560403, 
         0.22909743029780916, -0.22370817841070817, 0.22276350225740021, 
         -0.20654999242204575, 0.2169014379707668, -0.19161081599417318, 
         0.2114519169280712, -0.1784952515689809, 0.20636591553717418, 
         -0.1668957225320185, 0.20067806012731856, -0.15462398600687022, 
         0.19539716770023383, -0.14386282080342452, 0.19047447953905652, 
         -0.13435571624872938, 0.1846528592036232, -0.12373583290573131, 
         0.1792784058723679, -0.11450400531053491, 0.17429430319325642, 
         -0.10641118136186915, 0.16965366386523534, -0.09926399265095906, 
         0.1644687239462504, -0.09170223870144609, 0.15966833301024577, 
         -0.08508322608646261, 0.15520581905068526, -0.0792457291939122, 
         0.1499351670658275, -0.07272640279637052, 0.1450861229396252, 
         -0.06707245125362109, 0.14060426804734383, -0.06212733886358272, 
         0.13644477123461973, -0.05776974915554428, 0.13181024428782814, 
         -0.05316731742134924, 0.12753539787729393, -0.04914941042626166, 
         0.12357574402070427, -0.04561516391568519, 0.11893687271799168, 
         -0.04169548002912524, 0.1146858602606766, -0.03830504098065827, 
         0.1107717524716245, -0.03534751460236705, 0.10715262728782338, 
         -0.03274828630116326, 0.10379361024577506, -0.030448699951611417, 
         0.1006654053553251, -0.028402045602943137, 0.09774316641601988, 
         -0.02657067321447996, 0.09549847009150314, -0.02521552081322901, 
         0.09336649013253233, -0.02396909366730999, 0.09133832639751482, 
         -0.02281949257997841, 0.08940604014224415, -0.02175645998127155, 
         0.08756251191927335, -0.02077110319959103, 0.08704220099777898, 
         -0.020498036560251554, 0.0865287031172983, -0.020230702572241913, 
         0.08602187584731155, -0.019968933842493885, 0.08552158147883539, 
         -0.019712569629453676, 0.08547190606066216, -0.019687224389028313, 
         0.0854222944872667, -0.019661931483143887, 0.08537274662687912, 
         -0.01963669075530499, 0.08534798860849073, -0.019624085874176557, 
         0.08532324648238697, -0.019611493997452246, 0.0852985044379997, 
         -0.019598907067836846, 0.08527377826722557, -0.019586333108144567, 
         0.08522435763402807, -0.019561215983216057, 0.0851750003522388, 
         -0.019536150579045466, 0.08512570629146586, -0.019511136741386995, 
         0.08463621518412831, -0.019263802500671345, 0.08415290497985352, 
         -0.01902146475259384, 0.08367565181031315, -0.018783983481722406, 
         0.08320433477352857, -0.018551223550957528, 0.08152987374754905, 
         -0.017738396144949608, 0.0799278318851587, -0.016981121991666598, 
         0.07839333249301508, -0.016274240963070945, 0.07692193575248744, 
         -0.015613191242703505, 0.07550960128172565, -0.014993931951874182, 
         0.07349781053892515, -0.01413730484678876, 0.07159908978376042, 
         -0.013355900881102547, 0.06980366314985931, -0.012640837783794262, 
         0.06810290408990195, -0.01198454614030221, 0.06648915145971344, 
         -0.011380534935447308, 0.06495557718019349, -0.010823206400949014, 
         0.06291834736164034, -0.010107641797643636, 0.061013343731975776, 
         -0.009463775488138546, 0.05922756009195457, -0.00888205390656935, 
         0.05754969163872266, -0.008354498561757734, 0.055969859876067035, 
         -0.007874398487171407, 0.054479389785653025, -0.007436072439761302, 
         0.05269699578637608, -0.006930312980349328, 0.05103360756436605, 
         -0.006476204847926129, 0.049477293262552034, -0.006066777462580613, 
         0.04801770180141572, -0.005696219065811872, 0.046645805493871685, 
         -0.00535964928444798, 0.04535369138056924, -0.005052943450678048, 
         0.04379859841303442, -0.004696924194683291, 0.042351038506508426, 
         -0.0043782665108120875, 0.04099993402706103, -0.004091818950227696, 
         0.03973570670753986, -0.0038333013964969346, 0.038550029265603386, 
         -0.0035991315677878027, 0.037435624730912144, 
         -0.0033862915366638068, 0.03609405249169606, -0.003139314295071171, 
         0.034848480239930146, -0.002918983477625944, 0.03368874251030616, 
         -0.0027215405873334277, 0.03260608119893014, -0.002543875056777954, 
         0.03159290798453862, -0.002383393100959694, 0.030642613124449776, 
         -0.002237917292390627, 0.029496775868956576, -0.0020689604049467917, 
         0.0284357822083976, -0.0019187634588311538, 0.027450388640699402, 
         -0.0017846149252045264, 0.026532662757793875, 
         -0.0016642807101293398, 0.025675757241465445, 
         -0.0015559056799121518, 0.02487372879787971, -0.0014579386789922533, 
         0.02390444912146014, -0.0013439927401557425, 0.023009434822055946, 
         -0.001243079153335708, 0.022180354276982926, -0.001153264813947712, 
         0.02141008828791234, -0.0010729649188295344, 0.02069251628804376, 
         -0.0010008696012024803, 0.020022345943830525, 
         -0.0009358884337533501, 0.019209971733568313, 
         -0.0008601510205493184, 0.018462015672630675, 
         -0.0007933444971549473, 0.017771029747261684, 
         -0.0007341075160092441, 0.01724087512774926, -0.000690267561886129, 
         0.016741907516291822, -0.0006502787171511716}}, {Automatic}]}}
 
A1[nq_, Q2_, fi_] := (1/Pi)*NIntegrate[tr12d[nq, Exp[w]]/
        (1 + (Q2/LL2[nq])*Exp[I*fi]*Exp[-w]), {w, lns0[nq], lnsm}, 
       MaxRecursion -> 30, WorkingPrecision -> 16] + 
     f1ss0[nq]/((Q2/LL2[nq])*Exp[I*fi] + s1s0[nq]) + 
     f2ss0[nq]/((Q2/LL2[nq])*Exp[I*fi] + s2s0[nq]) + J0[nq] - 
     (Q2/LL2[nq])*Exp[I*fi]*J1[nq] + (Q2/LL2[nq])^2*Exp[I*2*fi]*J2[nq]
 
A1 /: A1::usage = "A1[Nf,Q2,fi] = A(Q^2) = A_1(Q^2), where Q^2=Q2 Exp[I fi] \
(Q2=|Q^2|, in GeV^2) and A(Q^2) is the 2danQCD analog of the underlying \
perturbative coupling a(Q^2) = a2d[Nf,Q2/LL2[Nf],fi] (see: a2d); at high \
|Q^2|> 1 GeV^2 the two couplings practically coincide"
 
tr12d[nq_, s_] := If[s != 0, If[s > 0, Im[a2d[nq, s, -Pi]], 
      Im[(1/c1[nq])*(1/(1 - c22d[nq]/c1[nq]^2 + ProductLog[-1, 
           zu2d[nq, -s, 0]]))]], 0]
 
tr12d /: tr12d::usage = "tr12d[Nf,s] is the discontinuity function \
rho_1^{(pt)}(sigma)=Im a(Q^2=-sigma - I epsilon) of the underlying pQCD \
coupling a(Q^2) (see: a2d) in the considered Nf-flavor 2danQCD, where \
s=sigma/LL2[Nf] (see: LL2)"
 
lns0[3] = 3.342861804649192
 
lns0[4] = 3.342861804649192
 
lns0[5] = 3.342861804649192
 
lns0[6] = 3.342861804649192
 
lnsm = 92.10340371976183
 
J0[3] = 0.00466999455549764309108360437157522531`16.
 
J0[4] = 0.00507753887074683231808795530064138381`16.
 
J0[5] = 0.00557750058581599984275926991902285832`16.
 
J0[6] = 0.00622202997429364348515952541839318206`16.
 
J1[3] = 4.84346028881497324758771242290045`16.*^-45
 
J1[4] = 5.29827956855580615903950025463651`16.*^-45
 
J1[5] = 5.87569345738507775148605975808154`16.*^-45
 
J1[6] = 6.66537630861582027402462127582015`16.*^-45
 
J2[3] = 2.446809702456800300219863101055`16.*^-85
 
J2[4] = 2.676734109529465657890833417244`16.*^-85
 
J2[5] = 2.968727412713396780222296111734`16.*^-85
 
J2[6] = 3.368285035422341381188279148526`16.*^-85
 
tA2[nq_, Q2_, fi_] := (1/be0[nq]/Pi)*((Q2/LL2[nq])*Exp[I*fi])*
     (NIntegrate[Exp[-w]*(tr12d[nq, Exp[w]]/
         (1 + (Q2/LL2[nq])*Exp[I*fi]*Exp[-w])^2), {w, lns0[nq], lnsm}, 
       MaxRecursion -> 30, WorkingPrecision -> 16] + 
      (f1ss0[nq]*Pi)/((Q2/LL2[nq])*Exp[I*fi] + s1s0[nq])^2 + 
      (f2ss0[nq]*Pi)/((Q2/LL2[nq])*Exp[I*fi] + s2s0[nq])^2 + J1[nq] - 
      2*(Q2/LL2[nq])*Exp[I*fi]*J2[nq] + 3*(Q2/LL2[nq])^2*Exp[I*2*fi]*J3[nq])
 
tA2 /: tA2::usage = "tA2[Nf,Q2,fi]={tilde A}_2(Q^2) for Q^2=Q2 Exp[I fi] \
(Q2=|Q^2|, in GeV^2); this is the first logarithmic derivative of A_1 (see: \
A1): {tilde A}_2(Q^2) = (-1)/(beta_0 1!) (d/d ln Q^2) A_1(Q^2)"
 
J3[3] = 1.63689500729422241956655214488`16.*^-125
 
J3[4] = 1.79074905870213668885243599517`16.*^-125
 
J3[5] = 1.98615801578900018756979502896`16.*^-125
 
J3[6] = 2.25360240588900587649899107709`16.*^-125
 
tA3[nq_, Q2_, fi_] := (1/2/be0[nq]^2/Pi)*((Q2/LL2[nq])*Exp[I*fi])*
     (NIntegrate[Exp[-w]*(tr12d[nq, Exp[w]]/
         (1 + (Q2/LL2[nq])*Exp[I*fi]*Exp[-w])^2), {w, lns0[nq], lnsm}, 
       MaxRecursion -> 30, WorkingPrecision -> 16] - 
      2*NIntegrate[Exp[-w]*(tr12d[nq, Exp[w]]/
          (1 + (Q2/LL2[nq])*Exp[I*fi]*Exp[-w])^3), {w, lns0[nq], lnsm}, 
        MaxRecursion -> 30, WorkingPrecision -> 16] + 
      (f1ss0[nq]*Pi)*((-s1s0[nq] + (Q2/LL2[nq])*Exp[I*fi])/
        (s1s0[nq] + (Q2/LL2[nq])*Exp[I*fi])^3) + (f2ss0[nq]*Pi)*
       ((-s2s0[nq] + (Q2/LL2[nq])*Exp[I*fi])/
        (s2s0[nq] + (Q2/LL2[nq])*Exp[I*fi])^3))
 
tA3 /: tA3::usage = "tA3[Nf,Q2,fi]={tilde A}_3(Q^2) for Q^2=Q2 Exp[I fi]; \
this is the second logarithmic derivative of A_1 (see: A1): {tilde A}_3(Q^2) \
= (-1)^2/(beta_0^2 2!) (d/d ln Q^2)^2 A_1(Q^2)"
 
tA4[nq_, Q2_, fi_] := (1/6/be0[nq]^3/Pi)*((Q2/LL2[nq])*Exp[I*fi])*
     (NIntegrate[Exp[-w]*(1 - 4*Exp[-w]*(Q2/LL2[nq])*Exp[I*fi] + 
         Exp[-2*w]*(Q2/LL2[nq])^2*Exp[2*I*fi])*(tr12d[nq, Exp[w]]/
         (1 + (Q2/LL2[nq])*Exp[I*fi]*Exp[-w])^4), {w, lns0[nq], lnsm}, 
       MaxRecursion -> 30, WorkingPrecision -> 16] + 
      (f1ss0[nq]*Pi)*((s1s0[nq]^2 - 4*s1s0[nq]*(Q2/LL2[nq])*Exp[I*fi] + 
         (Q2/LL2[nq])^2*Exp[2*I*fi])/(s1s0[nq] + (Q2/LL2[nq])*Exp[I*fi])^4) + 
      (f2ss0[nq]*Pi)*((s2s0[nq]^2 - 4*s2s0[nq]*(Q2/LL2[nq])*Exp[I*fi] + 
         (Q2/LL2[nq])^2*Exp[2*I*fi])/(s2s0[nq] + (Q2/LL2[nq])*Exp[I*fi])^4))
 
tA4 /: tA4::usage = "tA4[Nf,Q2,fi]={tilde A}_4(Q^2) for Q^2=Q2 Exp[I fi]; \
this is the third logarithmic derivative of A_1 (see: A1): {tilde A}_4(Q^2) = \
(-1)^3/(beta_0^3 3!) (d/d ln Q^2)^3 A_1(Q^2)"
 
tA[nq_, nupp_, Q2_, fi_] := tA2d[nq, nupp, Q2, fi]
 
tA2d[nq_, nupp_, Q2_, fi_] := If[Abs[nupp] <= 2., 
     tAMA1l[nq, nupp, Q2/LL2[nq], fi] + 
      (-(Pi*be0[nq]^(nupp - 1)*Gamma[nupp])^(-1))*
       ((Pi*(f1ss0[nq]/s1s0[nq]))*PolyLog[-nupp + 1, (-s1s0[nq])*Exp[(-I)*fi]*
           (LL2[nq]/Q2)] + (Pi*(f2ss0[nq]/s2s0[nq]))*PolyLog[-nupp + 1, 
          (-s2s0[nq])*Exp[(-I)*fi]*(LL2[nq]/Q2)] - 
        NIntegrate[(tr1MA1l[nq, ssig]/ssig)*PolyLog[-nupp + 1, 
           (-ssig)*Exp[(-I)*fi]*(LL2[nq]/Q2)], {ssig, 0, s0[nq]}] + 
        NIntegrate[((tr12d[nq, ssig] - tr1MA1l[nq, ssig])/ssig)*
          PolyLog[-nupp + 1, (-ssig)*Exp[(-I)*fi]*(LL2[nq]/Q2)], 
         {ssig, s0[nq], 10*s0[nq]}] + NIntegrate[
         (tr12d[nq, Exp[w]] - tr1MA1l[nq, Exp[w]])*PolyLog[-nupp + 1, 
           (-Exp[w])*Exp[(-I)*fi]*(LL2[nq]/Q2)], {w, Log[10*s0[nq]], 10.^9}, 
         MaxRecursion -> 30, WorkingPrecision -> 16]), 
     tAMA1l[nq, nupp, Q2/LL2[nq], fi] + 
      (-(Pi*be0[nq]^(nupp - 1)*Gamma[nupp])^(-1))*
       ((Pi*(f1ss0[nq]/s1s0[nq]))*PolyLog[-nupp + 1, (-s1s0[nq])*Exp[(-I)*fi]*
           (LL2[nq]/Q2)] + (Pi*(f2ss0[nq]/s2s0[nq]))*PolyLog[-nupp + 1, 
          (-s2s0[nq])*Exp[(-I)*fi]*(LL2[nq]/Q2)] - 
        NIntegrate[(tr1MA1l[nq, ssig]/ssig)*PolyLog[-nupp + 1, 
           (-ssig)*Exp[(-I)*fi]*(LL2[nq]/Q2)], {ssig, 0, s0[nq]}] + 
        NIntegrate[((tr12d[nq, ssig] - tr1MA1l[nq, ssig])/ssig)*
          PolyLog[-nupp + 1, (-ssig)*Exp[(-I)*fi]*(LL2[nq]/Q2)], 
         {ssig, s0[nq], 10*s0[nq]}] + NIntegrate[
         (tr12d[nq, Exp[w]] - tr1MA1l[nq, Exp[w]])*PolyLog[-nupp + 1, 
           (-Exp[w])*Exp[(-I)*fi]*(LL2[nq]/Q2)], {w, Log[10*s0[nq]], 10.^6}, 
         MaxRecursion -> 30, WorkingPrecision -> 16])]
 
tA2d /: tA2d::usage = "tA2d[Nf,nu,Q2,fi]= {tilde A}_{nu}(Q^2), with Q^2=Q2 \
Exp[I fi]  (Q2=|Q^2|); {tilde A}_{nu} can have nu noninteger; {tilde A}_{nu} \
is the generalization (to noninteger nu) of the logarithmic derivative \
(-1)^{nu-1}/(beta_0^{nu-1} (nu-1)! (d/d ln Q^2)^(nu-1) A_1(Q^2)"
 
tAMA1l[nq_, nu_, u_, fi_] := If[Abs[u*Exp[I*fi] - 1] > Delt[nu], 
     (1/be0[nq]^nu)*(1/(Log[u] + I*fi)^nu - 
       PolyLog[-nu + 1, (1/u)*Exp[(-I)*fi]]/Gamma[nu]), 
     (1/be0[nq]^nu)*Limt1l[u, nu, fi]]
 
Delt[nu_] := 10.^(-2) + (HeavisideTheta[nu - 3]*(nu - 3)*2)/10.^2
 
Limt1l[u_, nu_, fi_] := -(Zeta[1 - nu]/Gamma[nu]) + 
     ((u*Exp[I*fi] - 1)*Zeta[-nu])/Gamma[nu] - 
     ((u*Exp[I*fi] - 1)^2*(Zeta[-1 - nu] + Zeta[-nu]))/(2*Gamma[nu])
 
tr1MA1l[nq_, s_] := (Pi/be0[nq])*(1/(Log[s]^2 + Pi^2))
 
AN4[nq_, M_, nupp_, Q2_, fi_] := A2dN4[nq, M, nupp, Q2, fi]
 
A2dN4[nq_, 4, nupp_, Q2_, fi_] := tA2d[nq, 4 + nupp, Q2, fi]
 
A2dN4[nq_, 3, nupp_, Q2_, fi_] := tA2d[nq, 3 + nupp, Q2, fi] + 
     tk1f[nq, 3 + nupp]*tA2d[nq, 4 + nupp, Q2, fi]
 
A2dN4[nq_, 2, nupp_, Q2_, fi_] := tA2d[nq, 2 + nupp, Q2, fi] + 
     tk1f[nq, 2 + nupp]*tA2d[nq, 3 + nupp, Q2, fi] + 
     tk2f[nq, 2 + nupp]*tA2d[nq, 4 + nupp, Q2, fi]
 
A2dN4[nq_, 1, nupp_, Q2_, fi_] := tA2d[nq, 1 + nupp, Q2, fi] + 
     tk1f[nq, 1 + nupp]*tA2d[nq, 2 + nupp, Q2, fi] + 
     tk2f[nq, 1 + nupp]*tA2d[nq, 3 + nupp, Q2, fi] + 
     tk3f[nq, 1 + nupp]*tA2d[nq, 4 + nupp, Q2, fi]
 
A2dN4[nq_, 0, nupp_, Q2_, fi_] := tA2d[nq, nupp, Q2, fi] + 
     tk1f[nq, nupp]*tA2d[nq, 1 + nupp, Q2, fi] + 
     tk2f[nq, nupp]*tA2d[nq, 2 + nupp, Q2, fi] + 
     tk3f[nq, nupp]*tA2d[nq, 3 + nupp, Q2, fi] + 
     tk4f[nq, nupp]*tA2d[nq, 4 + nupp, Q2, fi]
 
tk1f[nq_, s_] := -k1f[nq, s]
 
k1f[nq_, s_] := c1[nq]*s*(-1 + HarmonicNumber[s])
 
tk2f[nq_, s_] := k1f[nq, s]*k1f[nq, s + 1] - k2f[nq, s]
 
k2f[nq_, s_] := ((s - 1)*(s/2))*c22d[nq] + 
     c1[nq]^2*((s*((s + 1)/2))*((PolyGamma[s + 1] - PolyGamma[1])^2 + 
         (PolyGamma[1, s + 1] - PolyGamma[1, 1]) + 1) - 
       s^2*(PolyGamma[s + 1] - PolyGamma[1]))
 
tk3f[nq_, s_] := (-k1f[nq, s])*k1f[nq, s + 1]*k1f[nq, s + 2] + 
     k1f[nq, s]*k2f[nq, s + 1] + k1f[nq, s + 2]*k2f[nq, s] - k3f[nq, s]
 
k3f[nq_, s_] := (s*(s - 1)*((s + 4)/12))*c32d[nq] + 
     (s*(s + 1)*((s + 2)/2))*c1[nq]*c22d[nq]*((s/(s + 2))*ES[1, s + 1] + 
       2/(s + 1) + 1/(s + 2) - 11/6) + (s*(s + 1)*((s + 2)/6))*c1[nq]^3*
      (Z3[s + 2] - 3*Z2[s + 1] + 3*ES[1, s] - 1)
 
ES[m_, n_] := ((-1)^(m - 1)/(m - 1)!)*(PolyGamma[m - 1, n + 1] - 
      PolyGamma[m - 1, 1])
 
Z3[n_] := ES[1, n]^3 - 3*ES[1, n]*ES[2, n] + 2*ES[3, n]
 
Z2[n_] := ES[1, n]^2 - ES[2, n]
 
tk4f[nq_, s_] := k1f[nq, s]*k1f[nq, s + 1]*k1f[nq, s + 2]*k1f[nq, s + 3] - 
     k1f[nq, s]*k1f[nq, s + 1]*k2f[nq, s + 2] - k1f[nq, s]*k1f[nq, s + 3]*
      k2f[nq, s + 1] - k1f[nq, s + 2]*k1f[nq, s + 3]*k2f[nq, s] + 
     k2f[nq, s]*k2f[nq, s + 2] + k1f[nq, s]*k3f[nq, s + 1] + 
     k1f[nq, s + 3]*k3f[nq, s] - k4f[nq, s]
 
k4f[nq_, s_] := ((s*(s + 1)*(s + 2)*(s + 3))/6)*
     (c42d[nq]*(1/12 - 2/((s + 1)*(s + 2)*(s + 3))) + 
      c22d[nq]^2*(5/12 - (3*s^2 + 6*s + 1)/((s + 1)*(s + 2)*(s + 3))) + 
      (c1[nq]*(c32d[nq]/2))*((1 - 6/((s + 2)*(s + 3)))*Z1[s + 1] + 
        4/(s + 1) - 1/(s + 2) - 1/(s + 3) - 13/6) + (c1[nq]^2*(c22d[nq]/2))*
       ((3*((s + 1)/(s + 3)))*Z2[s + 2] + (12/(s + 2) + 6/(s + 3) - 11)*
         Z1[s + 1] - 6/(s + 1) - 5/(s + 2) - 11/(s + 3) + 38/3) + 
      (c1[nq]^4/4)*(Z4[s + 3] - 4*Z3[s + 2] + 6*Z2[s + 1] - 4*Z1[s] + 1))
 
Z1[n_] := ES[1, n]
 
Z4[n_] := ES[1, n]^4 - 6*ES[1, n]^2*ES[2, n] + 3*ES[2, n]^2 + 
     8*ES[1, n]*ES[3, n] - 6*ES[4, n]
 
AN3[nq_, M_, nupp_, Q2_, fi_] := A2dN3[nq, M, nupp, Q2, fi]
 
A2dN3[nq_, 3, nupp_, Q2_, fi_] := tA2d[nq, 3 + nupp, Q2, fi]
 
A2dN3[nq_, 2, nupp_, Q2_, fi_] := tA2d[nq, 2 + nupp, Q2, fi] + 
     tk1f[nq, 2 + nupp]*tA2d[nq, 3 + nupp, Q2, fi]
 
A2dN3[nq_, 1, nupp_, Q2_, fi_] := tA2d[nq, 1 + nupp, Q2, fi] + 
     tk1f[nq, 1 + nupp]*tA2d[nq, 2 + nupp, Q2, fi] + 
     tk2f[nq, 1 + nupp]*tA2d[nq, 3 + nupp, Q2, fi]
 
A2dN3[nq_, 0, nupp_, Q2_, fi_] := tA2d[nq, nupp, Q2, fi] + 
     tk1f[nq, nupp]*tA2d[nq, 1 + nupp, Q2, fi] + 
     tk2f[nq, nupp]*tA2d[nq, 2 + nupp, Q2, fi] + 
     tk3f[nq, nupp]*tA2d[nq, 3 + nupp, Q2, fi]
 
A5N4[nq_, Q2_, fi_] := tA2d[nq, 5, Q2, fi]
 
A4N4[nq_, Q2_, fi_] := tA4[nq, Q2, fi] + tk1f[nq, 4]*tA2d[nq, 5, Q2, fi]
 
A3N4[nq_, Q2_, fi_] := tA3[nq, Q2, fi] + tk1f[nq, 3]*tA4[nq, Q2, fi] + 
     tk2f[nq, 3]*tA2d[nq, 5, Q2, fi]
 
A2N4[nq_, Q2_, fi_] := tA2[nq, Q2, fi] + tk1f[nq, 2]*tA3[nq, Q2, fi] + 
     tk2f[nq, 2]*tA4[nq, Q2, fi] + tk3f[nq, 2]*tA2d[nq, 5, Q2, fi]
 
A1N4[nq_, Q2_, fi_] := A1[nq, Q2, fi]
 
A4N3[nq_, Q2_, fi_] := tA4[nq, Q2, fi]
 
A3N3[nq_, Q2_, fi_] := tA3[nq, Q2, fi] + tk1f[nq, 3]*tA4[nq, Q2, fi]
 
A2N3[nq_, Q2_, fi_] := tA2[nq, Q2, fi] + tk1f[nq, 2]*tA3[nq, Q2, fi] + 
     tk2f[nq, 2]*tA4[nq, Q2, fi]
 
A1N3[nq_, Q2_, fi_] := A1[nq, Q2, fi]
 
AN2[nq_, M_, nupp_, Q2_, fi_] := A2dN2[nq, M, nupp, Q2, fi]
 
A2dN2[nq_, 2, nupp_, Q2_, fi_] := tA2d[nq, 2 + nupp, Q2, fi]
 
A2dN2[nq_, 1, nupp_, Q2_, fi_] := tA2d[nq, 1 + nupp, Q2, fi] + 
     tk1f[nq, 1 + nupp]*tA2d[nq, 2 + nupp, Q2, fi]
 
A2dN2[nq_, 0, nupp_, Q2_, fi_] := tA2d[nq, nupp, Q2, fi] + 
     tk1f[nq, nupp]*tA2d[nq, 1 + nupp, Q2, fi] + 
     tk2f[nq, nupp]*tA2d[nq, 2 + nupp, Q2, fi]
 
A3N2[nq_, Q2_, fi_] := tA3[nq, Q2, fi]
 
A2N2[nq_, Q2_, fi_] := tA2[nq, Q2, fi] + tk1f[nq, 2]*tA3[nq, Q2, fi]
 
A1N2[nq_, Q2_, fi_] := A1[nq, Q2, fi]
 
AN1[nq_, M_, nupp_, Q2_, fi_] := A2dN1[nq, M, nupp, Q2, fi]
 
A2dN1[nq_, 1, nupp_, Q2_, fi_] := tA2d[nq, 1 + nupp, Q2, fi]
 
A2dN1[nq_, 0, nupp_, Q2_, fi_] := tA2d[nq, nupp, Q2, fi] + 
     tk1f[nq, nupp]*tA2d[nq, 1 + nupp, Q2, fi]
 
A2N1[nq_, Q2_, fi_] := tA2[nq, Q2, fi]
 
A1N1[nq_, Q2_, fi_] := A1[nq, Q2, fi]
 
rtaudBG22 = 0.19517068796961218
 
rtaudBG22 /: rtaudBG22::usage = "rtaudBG22=r_{tau}(D=0) calculated in the \
considered 2danQCD with the diagonal [2/2] Baker-Gammel approach."
 
rtauTPS = 0.17984560907740513
 
rtauTPS /: rtauTPS::usage = "rtauTPS=r_{tau}(D=0) calculated in the \
considered 2danQCD with the truncated perturbation series approach."
 
rtauLBbLB = 0.20102433230210737
 
rtauLBbLB /: rtauLBbLB::usage = "rtauLBbLB=r_{tau}(D=0) calculated (in the \
considered Lambert renormalization scheme) in the considered 2danQCD with the \
'leading-beta_0 (LB) + beyond-the leading-beta_0 (bLB)' approach with Nf=3. \
Note that here rtauLBbLB=0.2010, or 0.2030, or 0.1990"
 
d1bar = 1/12 + (9*(11/2 - 4*Zeta[3]))/4
 
d2bar = -12521/288 + (81*(302/9 - (76*Zeta[3])/3))/16 + 13*Zeta[3] + 
     (3*(401 - 424*Zeta[3] + 200*Zeta[5]))/16
 
d3bar = 144939499/20736 - (5693495*Zeta[3])/864 + (5445*Zeta[3]^2)/8 + 
     9*(1045381/15552 - (40655*Zeta[3])/864 + (5*Zeta[3]^2)/6 - 
       (260*Zeta[5])/27) + 27*(-6131/5832 + (203*Zeta[3])/324 + 
       (5*Zeta[5])/18) + (65945*Zeta[5])/288 - (7315*Zeta[7])/48 + 
     3*(-13044007/10368 + (12205*Zeta[3])/12 - 55*Zeta[3]^2 + 
       (29675*Zeta[5])/432 + (665*Zeta[7])/72)
 
td1bar = 1/12 + (9*(11/2 - 4*Zeta[3]))/4
 
td2bar = -12521/288 - (16*(1/12 + (9*(11/2 - 4*Zeta[3]))/4))/9 + 
     (81*(302/9 - (76*Zeta[3])/3))/16 + 13*Zeta[3] + 
     (3*(401 - 424*Zeta[3] + 200*Zeta[5]))/16
 
td3bar = 144939499/20736 + (8891*(1/12 + (9*(11/2 - 4*Zeta[3]))/4))/2592 - 
     (5693495*Zeta[3])/864 + (5445*Zeta[3]^2)/8 + 
     9*(1045381/15552 - (40655*Zeta[3])/864 + (5*Zeta[3]^2)/6 - 
       (260*Zeta[5])/27) + 27*(-6131/5832 + (203*Zeta[3])/324 + 
       (5*Zeta[5])/18) + (65945*Zeta[5])/288 - 
     (40*(-12521/288 + (81*(302/9 - (76*Zeta[3])/3))/16 + 13*Zeta[3] + 
        (3*(401 - 424*Zeta[3] + 200*Zeta[5]))/16))/9 - (7315*Zeta[7])/48 + 
     3*(-13044007/10368 + (12205*Zeta[3])/12 - 55*Zeta[3]^2 + 
       (29675*Zeta[5])/432 + (665*Zeta[7])/72)
 
d1 = 1.6398212048969851
 
d2 = 15.742079297915776
 
d3 = 83.55174899975434
 
td1 = 1.6398212048969851
 
td2 = 12.826841600321135
 
td3 = 34.57868801343341
 
td1bLB = 0.08333333333333326
 
td2bLB = -2.8843605694609113
 
td3bLB = 9.746656167968794
