In[]:=
deploy
Tue 29 Jun 2021 14:22:30
In[]:=
eps=
-2
10
;​​logf[t_]:=LogGamma2-
1
p
+LogGamma[t+1]-LogGammat-
1
p
+2;​​​​inv[p_]:=Log10t/.FindRootLogGamma2-
1
p
+LogGamma[t+1]-LogGammat-
1
p
+2Log[eps],{t,100},WorkingPrecision100;​​DiscretePlotinv[p],5
-3/4
(p-1)
,{p,11/10,30/10,1/10},FillingAxis,FillingStyle->Directive[Opacity[0.5],Orange],AxesLabel{"p","
log
10
g"},PlotLegendsg,5
-3/4
(p-1)
,ImageSizeLarge,PlotRangeAll
Out[]=
g
5
3/4
(p-1)
In[]:=
vals=Table[{p,SetPrecision[10^inv[p],5]},{p,2,6}];​​TableForm[vals,TableHeadings{None,{"p","g(p)"}}]
Out[]//TableForm=
p
g(p)
2
7853.2
3
856.89
4
413.88
5
288.44
6
232.49