diff options
Diffstat (limited to 'schofield.wl')
-rw-r--r-- | schofield.wl | 26 |
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)! |