a=1.0 b=0.5 N=3 tau=1: kappa=8.75693909432999 gamma=-1.0 4kappa/(2+gamma)=35.0277563773199
a=1.0 b=0.5 N=4 tau=1: kappa=2.5 c1=0.8 c2=0.0 k110=c1/kappa=0.32 k807=0.16 k809=1.0
a=1.0 b=0.5 N=5 tau=1: kappa=1.70193433372847 gamma=0.641234989147323 4kappa/(2+gamma)=2.57748264084281
a=1.0 b=0.5 N=5 tau=2: kappa=120.679620832013 gamma=-1.14206515520185 4kappa/(2+gamma)=562.651681832112
a=1.0 b=0.5 N=6 tau=1: kappa=1.42302494707577 c1=1.56161859761401 c2=1.56161859761401 k110=c1/kappa=1.09739368998628 k807=0.823045267489712 k809=1.66666666666667
a=1.0 b=0.8 N=7 tau=1: kappa=1.23766239920224 gamma=1.24499315246412 4kappa/(2+gamma)=1.52562713207873
a=1.0 b=0.8 N=7 tau=3: kappa=71.1728320169094 gamma=-1.46828432899947 4kappa/(2+gamma)=535.420232267997
a=1.0 b=0.8 N=8 tau=1: kappa=1.17554020779209 c1=5.09822521273348 c2=7.19466395791987 k110=c1/kappa=4.33692116946728 k807=3.69853669574837 k809=5.11098855399564
a=1.0 b=0.8 N=8 tau=3: kappa=8.99315377989274 c1=4.53231943702993 c2=-6.05300048917743 k110=c1/kappa=0.503974417424449 k807=0.0837203070967151 k809=6.77103178116885
a=2.0 b=1.2 N=9 tau=2: kappa=1.85121609542666 gamma=0.361039135450947 4kappa/(2+gamma)=3.13627346134283
a=2.0 b=1.2 N=10 tau=3: kappa=3.80129861431242 c1=104.99420563142 c2=-69.2576838347406 k110=c1/kappa=27.6206150277439 k807=9.25543752983702 k809=133.739879790902
a=1.0 b=0.3 N=11 tau=4: kappa=79.4296875225325 gamma=-1.0431611754154 4kappa/(2+gamma)=332.050437259444
a=1.0 b=0.3 N=12 tau=5: kappa=6291.56427272844 c1=0.220857477957895 c2=-0.231264235212301 k110=c1/kappa=3.51037466016565e-5 k807=8.36241638727722e-6 k809=331.016058887005
a=3.0 b=1.0 N=13 tau=1: kappa=1.08993326845643 gamma=1.70708388502909 4kappa/(2+gamma)=1.17605460492338
a=3.0 b=1.0 N=14 tau=3: kappa=2.24146894908829 c1=300.513792517333 c2=152.452745835689 k110=c1/kappa=134.070022535698 k807=84.0386759737215 k809=422.225691348668

worst discrepancies (relative unless obvious):
  Lemma2.2: N*2Delta = 4 tau K                            5.96e-44
  b2 = b cn Delta                                         1.65e-44
  a2 = a cn/dn Delta                                      3.13e-45
  kappa = dn/cn^2                                         1.34e-44
  geometric orbit = P(u0 + l h) or P(u0 - l h)            4.5e-43
  inner_l = sigma^-1 P(u_l+Delta)                         2.11e-42
  outer_l = sigma P(u_l+Delta)                            1.46e-39
  A = F(u)                                                3.02e-43
  A' = kappa F(u+D)                                       7.3e-43
  A'' = F(u+D)/kappa                                      2.41e-43
  A'^ominus = F/2 + G/4                                   2.93e-43
  A^otimes = (F/2+G/4)(u+D)/kappa                         7.03e-43
  ThmB: F(u+D) = F(u) (complex u)                         5.32e-47
  ThmB: G = gamma F (complex u)                           2.69e-48
  Lemma3.3: X(p-e) = -X(p+e), X=F                         4.54e-41
  Lemma3.3: X(p-e) = -X(p+e), X=G                         4.54e-41
  ThmB residue of F at iK' (eps=1e-12)                    2.07e-23
  geom k808 = 4kappa/(2+gamma)                            8.55e-48
  geom k810 = 4kappa/(2+gamma)                            8.47e-48
  ThmA: F(u)F(u+D) const (complex u)                      2.08e-43
  ThmA: F(u)G(u+D) const (complex u)                      2.05e-43
  ThmA: G(u)F(u+D) = F(u)G(u+D)                           2.08e-43
  geom k110 = c1/kappa                                    8.46e-44
  geom k807 = (c1/2+c2/4)/kappa                           8.57e-44
  geom k809 = kappa(c1/2+c2/4)                            8.66e-44
