Coefficient updates.

1 files changed, 19 insertions, 7 deletions

diff --git a/schofield.wl b/schofield.wl
index 98dc607..d3ad6d2 100644
--- a/schofield.wl
+++ b/schofield.wl
@@ -122,8 +122,20 @@ eqMid[F_, h_][m_] := D[
   {θ, m} ] / m! /. θ -> 1
 
 Φs = {
-  -1.19773338379799339,
+  -Gamma[1/3]Gamma[1/5]Gamma[7/15]/(2 π Gamma[2/3]Gamma[4/5]Gamma[8/15])(4 π^2 Gamma[13/16]^2 Gamma[3/4]/(Gamma[3/16]^2 Gamma[1/4]))^(8/15),
   -0.31881012489061,
+(*
+  Sin[π s(pp l - p k)/p]/Sin[π s/p(pp-p)](M π(ξ+1)Gamma[(2 ξ+2)/(3 ξ+6)]/(2^(2/3)Sqrt[3]Gamma[1/3]Gamma[ξ/(3 ξ+6)]))^(2 Δlk) Q[ξ, (ξ+1)l-ξ k] //. {
+    l -> 1, k -> 3, p -> 3, pp -> 4, ξ -> p/(pp-p), Δlk -> ((pp l-p k)^2-(pp-p)^2)/(4 p pp),
+    M -> 2^((ξ+5)/(3 ξ+6)) Sqrt[3] Gamma[1/3]Gamma[ξ/(3 ξ+6)]/π/Gamma[(2 ξ+2)/(3 ξ+6)](
+      π^2 Gamma[(3 ξ+4)/(4 ξ+4)]^2 Gamma[1/2+1/(ξ+1)]/Gamma[ξ/(4 ξ+4)]^2/Gamma[1/2-1/(ξ+1)]
+    )^((ξ+1)/(3 ξ + 6))
+  } /. Q -> Function[{ξ, η}, Exp[
+      NIntegrate[
+      (Sinh[(ξ+2)t]Sinh[t(η-1)]Sinh[t(η+1)]/Sinh[3(ξ+2)t]/Sinh[2(ξ+1)t]/Sinh[ξ t](Cosh[3(ξ+2)t]+Cosh[(ξ+4)t]-Cosh[(3 ξ+4)t]+Cosh[ξ t]+1)-(η^2-1)/(2 ξ(ξ+1))Exp[-4 t]) /t
+      ,{t, 0, Infinity}, WorkingPrecision->50]
+    ]], 
+*)
   Around[0.110886196683, 2.0 10^-12],
   Around[0.01642689465,  1.0 10^-11],
   Around[-2.639978 10^-4, 1.0 10^-10],
@@ -131,8 +143,8 @@ eqMid[F_, h_][m_] := D[
   Around[2.08865 10^-4,   1.0 10^-9],
   Around[-4.4819 10^-5,   1.0 10^-9],
   Around[3.194 10^-7,     1.0 10^-9],
-  4.31 10^-6,
-  -1.99 10^-6
+  Around[4.31 10^-6,      0.01 10^-6],
+  Around[-1.99 10^-6,     0.01 10^-6]
 }
 
 Gls = {
@@ -149,8 +161,8 @@ Gls = {
   Around[-6.79 10^2,   1.0],
   Around[5.34 10^3,    10.],
   Around[-4.66 10^4,   0.01 10^4],
-  4.46 10^5,
-  -4.66 10^6
+  Around[4.46 10^5, 0.01 10^5],
+  Around[-4.66 10^6, 0.01 10^6]
 }
 
 Ghs = {
@@ -166,9 +178,9 @@ Ghs = {
   0,
   Around[-2.5891 10^4,  2.0],
   0,
-  5.02 10^5,
+  Around[5.02 10^5, 0.01 10^5],
   0,
-  -1.04 10^7
+  Around[-1.04 10^7, 0.01 10^7]
 }
 
 dRule[sym_][f_, i_] := Derivative[i[[1]] - 1][sym][0] -> f (i[[1]] - 1)!