Get["~/doc/research/first_order_singularities/mma/schofield.wl"]

θ0Cas=(Sqrt[1.16951]);
hCas[θ_]:=(1-((θ^2)/((θ0Cas)^2)))θ(1-0.222389 (θ^2)-0.043547 (θ^4)-0.014809 (θ^6)-0.007168 (θ^8))
h0=a b ρ/.{a->((c2^2)/c4),b->(((((-c4)/(c2^3))))^(1/2)),ρ->2.00881}/.{c2->Ghs[[3]]2!,c4->Ghs[[5]]["Value"]4!}
m0=-Ghs[[3]]2! h0

totCas=13;
ds=invDerivativeList[totCas+1][θ|->((h0 hCas[θ])/(((t[θ]))^(Δ[2])))][θ0Cas];
dMs=m0 NestList[D[#,θ]&,(θ/((((θ^2)-1))^(1/8))),totCas]/.θ->θ0Cas;
datCas=Table[-Sum[dMs[[k+1]]BellY[n,k,ds[[;;n-k+1]]],{k,0,n}]/(n+1)!,{n,0,totCas}]

totCas=13;
ds=invDerivativeList[totCas+1][θ|->((h0 hCas[θ])/(((-t[θ]))^(Δ[2])))][0];
dMs=m0 NestList[D[#,θ]&,(θ/(((1-(θ^2)))^(1/8))),totCas]/.θ->0;
datCasH=Table[-Sum[dMs[[k+1]]BellY[n,k,ds[[;;n-k+1]]],{k,0,n}]/(n+1)!,{n,0,totCas}]

totCas=6;
ds=invDerivativeList[totCas+1][θ|->((t[θ])/(((h0 hCas[θ]))^(1/Δ[2])))][1];
dMs=m0 NestList[D[#,θ]&,(θ/(((h0 hCas[θ]))^15)),totCas]/.θ->1;
datCasM=Table[-Sum[dMs[[k+1]]BellY[n,k,ds[[;;n-k+1]]],{k,0,n}]/(n+1)!,{n,0,totCas}]

2/Δ[]-1/Δ[]-1

2/Δ[]-1

-Φs[[1]]//N

ParametricPlot[Evaluate[{{((t[θ])/((H[2][θ])^(1/Δ[]))),((Abs[t[θ]/((H[2][θ])^(1/Δ[]))])^(β[]))dGdξList[NF[2],H[2]][1,θ][[2]]}//.rules[H[2]]/.S[2],{((t[θ])/((H[9][θ])^(1/Δ[]))),((Abs[t[θ]/((H[9][θ])^(1/Δ[]))])^(β[]))dGdξList[NF[9],H[9]][1,θ][[2]]}//.rules[H[9]]/.S[9],
{((((θ^2)-1))/(((h0 hCas[θ]))^(1/Δ[]))),-m0 θ/(((h0 hCas[θ]))^(β[]/Δ[]))}}],{θ,0,1.4},AspectRatio->1/GoldenRatio,PlotPoints->100]

S[2]

m0/(((h0 hCas[1]))^(β[]/Δ[]))

Φs[[1]]

datCasM/.m0rule

totOurs=17;
ds=invDerivativeList[totOurs][θ|->((h1[θ]/.sol1)/(((t[θ]))^(Δ[2])))][θc/.sol1];
dFRegs=NestList[D[#,θ]&,RFReg[n1][θ]/.sol1,totOurs]/.θ->θc/.sol1;
dF=dFRegs+Table[Piecewise[{{AL m! Gamma[m-1] B^(m-1)/π,m>1},{0,True}}],{m,0,totOurs}];
dMs=Rest[NestList[D[#,θ]&,((fp[θ])/((t[θ])^2))-(1/(8π))Log[((t[θ])^2)],totOurs]]/.θ->θc/.Map[Derivative[#][fp][θc]->dF[[#+1]]&,Range[0,totOurs]]/.sol1;
datOurs=Table[Sum[dMs[[k]]BellY[n,k,ds[[;;n-k+1]]],{k,1,n}]/(n)!,{n,1,totOurs}]

datOurLow2=Chop[dGdξLowList[NF[2],H[2],Join[rules[H[2]],S[2]]][14]]

datOurLow9=Chop[dGdξLowList[NF[9],H[9],Join[rules[H[9]],S[9]]][14]]

datOurHigh2=Chop[dGdξList[NF[2],H[2],Join[rules[H[2]],S[2]]][14,0]]

datOurHigh9=Chop[dGdξList[NF[9],H[9],Join[rules[H[9]],S[9]]][14,0]]

datOurMid2=Chop[dΦdηList[NF[2],H[2],Join[rules[H[2]],S[2]]][14,1]]

datOurMid9=Chop[dΦdηList[NF[9],H[9],Join[rules[H[9]],S[9]]][14,1]]

Export["~/doc/research/first_order_singularities/paper/data/glow_series_ours_9.dat",Table[{i-1,datOurLow9[[i]]},{i,1,Length[datOurLow9]}]]

Export["~/doc/research/first_order_singularities/paper/data/glow_series_caselle.dat",Prepend[Table[{i,datCas[[i]]},{i,1,Length[datCas]}],{0,0}]]

Export["~/doc/research/first_order_singularities/paper/data/ghigh_caselle.dat",Prepend[Table[{i,datCasH[[i]]},{i,1,Length[datCasH]}],{0,0}]]

Export["~/doc/research/first_order_singularities/paper/data/ghigh_series_ours_2.dat",Thread[{Range[0,14],datOurHigh2}]]

Export["~/doc/research/first_order_singularities/paper/data/ghigh_series_ours_9.dat",Thread[{Range[0,14],datOurHigh9}]]

"~/doc/research/first_order_singularities/paper/data/ghigh_series_ours_9.dat"

Export["~/doc/research/first_order_singularities/paper/data/phi_series_ours_2.dat",Thread[{Range[0,14],datOurMid2}]]

Export["~/doc/research/first_order_singularities/paper/data/phi_series_ours_9.dat",Thread[{Range[0,14],datOurMid9}]]

Gls[[1]]

ListLogPlot[{
Abs@datNum,Abs@Prepend[datCas,Gls[[2]]]
}]

Export["~/doc/research/first_order_singularities/paper/data/Glow_series_numeric

datNum

convNum=Table[{i,#[[i]]/i!},{i,2,Length[#]}]&@Abs[Rest@datNum[[All,2]]]

ListPlot[{
Table[{1/i,#[[i]]/#[[i-1]]/i},{i,2,Length[#]}]&@Abs[Rest@datNum[[All,2]]],
Table[{1/i,#[[i]]/#[[i-1]]/i},{i,2,Length[#]}]&@Abs[datCas],
Table[{1/i,#[[i]]/#[[i-1]]/i},{i,2,Length[#]}]&@Rest[datOurs]
},PlotRange->{{0,0.5},All}]

datOurs=Chop[Abs[dFdξHighList[n1,h1][12]]/.sol1];

datNumHigh=Thread[{Range[0,Length[Ghs]-1],Abs[Ghs]((Range[Length[Ghs]]-1)!)}];

datCasHigh=NestList[ddξ[hCas],m0 (θ/(((1-(θ^2)))^(1/8))),11]/.θ->0/.m0rule

ListLogPlot[{
Abs@datNumHigh,Abs@datCasHigh,Thread[{Range[0,Length[datOurs]-1],datOurs}]
}]