Asymptotic statements from the closed forms (no eigenvalue computation)
gamma_3 = kappa_h g_inf(0) = 0.2527310099
  g_inf([4, 0, 0]) * 4 pi kappa_h |x| = 1.019830
  g_inf([8, 0, 0]) * 4 pi kappa_h |x| = 1.004108
  g_inf([16, 0, 0]) * 4 pi kappa_h |x| = 1.000988
  g_inf([6, 6, 6]) * 4 pi kappa_h |x| = 0.998451

(5.15) orders of the limits on the infinite cubic lattice (L -> inf first), anchor at 0, g1 = g2 = 1; c(r) as in mutual_cost.txt
  (a) midpoint (0,40,0), r << D: c(4->8)=+1.0341, c(8->16)=+0.9956
  (c) midpoint at the anchor: c(8->16)=-3.1350, c(16->32)=-3.0272, c(32->64)=-3.0065
  (b) x1=(3,0,0) fixed, x2=(3+r,0,0): 4 pi kappa_h r^2 [E12 - E12(r=inf)] -> -(x1-h).r_hat = -3: r=8: -2.2024, r=16: -2.5327, r=32: -2.7447, r=64: -2.8662

(5.23) a_mark on the infinite cubic lattice (L -> inf first), collinear: anchor 0, body d, d+1, source d+1+R'; prediction a_mark ~ (2 c_phi/Gamma_inf) (R/d)^2 with c_phi=1/kappa_a+gamma_3/kappa_h=1.681302, Gamma_inf=1.934033, R=R'+1/2, d->d+1/2
  R'=3 d=  8: a_mark=+2.354491e-01  prediction=+2.947882e-01  ratio=0.7987
  R'=3 d= 16: a_mark=+6.851490e-02  prediction=+7.823121e-02  ratio=0.8758
  R'=3 d= 32: a_mark=+1.813394e-02  prediction=+2.016421e-02  ratio=0.8993
  R'=3 d= 64: a_mark=+4.636675e-03  prediction=+5.119511e-03  ratio=0.9057
  R'=6 d=  8: a_mark=+6.421039e-01  prediction=+1.016718e+00  ratio=0.6315
  R'=6 d= 16: a_mark=+2.306867e-01  prediction=+2.698179e-01  ratio=0.8550
  R'=6 d= 32: a_mark=+6.571165e-02  prediction=+6.954595e-02  ratio=0.9449
  R'=6 d= 64: a_mark=+1.715213e-02  prediction=+1.765709e-02  ratio=0.9714

(5.24) a_bond on the infinite cubic lattice: pair (d,0,0),(d,2,0) -> translated by e_x, source (d+1+R',1,0)
  R'=3 d=  8: a_bond=+2.539632e-01  prediction=+3.239171e-01  ratio=0.7840
  R'=3 d= 16: a_bond=+7.600590e-02  prediction=+8.727102e-02  ratio=0.8709
  R'=3 d= 32: a_bond=+2.027305e-02  prediction=+2.258620e-02  ratio=0.8976
  R'=3 d= 64: a_bond=+5.193783e-03  prediction=+5.740509e-03  ratio=0.9048
  R'=6 d=  8: a_bond=+6.966915e-01  prediction=+1.117183e+00  ratio=0.6236
  R'=6 d= 16: a_bond=+2.562273e-01  prediction=+3.009960e-01  ratio=0.8513
  R'=6 d= 32: a_bond=+7.351920e-02  prediction=+7.789935e-02  ratio=0.9438
  R'=6 d= 64: a_bond=+1.922696e-02  prediction=+1.979890e-02  ratio=0.9711

(5.25) ring: a_mark for a body between anchor and source (d=10, d+1, source at d+1+3) -> 1 as L -> inf
  L=   32: a_mark=0.528868
  L=   64: a_mark=0.811924
  L=  256: a_mark=0.959032
  L= 1024: a_mark=0.990075
  L= 4096: a_mark=0.997538

(5.13) anchored line (L = inf): w(d) and Sigma(d) from the resolvent (5.9); predictions w(10d)-w(d) -> ln10/(2 pi sqrt(kappa_h)), Sigma(d)/d -> 2 ln2/(pi kappa_h^{3/2})
  d=    10: w=0.809467  Sigma/d=0.489480  Gamma=11.429
  d=   100: w=1.157998  Sigma/d=0.446296  Gamma=101.429  w(d)-w(d/10)=0.348531
  d=  1000: w=1.522451  Sigma/d=0.441776  Gamma=1001.429  w(d)-w(d/10)=0.364453
  d= 10000: w=1.888714  Sigma/d=0.441322  Gamma=10001.429  w(d)-w(d/10)=0.366263
  predictions: ln10/(2 pi)=0.366468, 2 ln2/pi=0.441271