level 1
我所想的是这样的
首先Plot[Exp[-x]*Sin[Gamma[x]]^2/Gamma[x]^2, {x, 0, 5}, PlotRange -> Full]得出Pt/P在正实轴上的值域大概是从0.到0.279
再
ListPlot[
ParallelTable[{ptp,
ptp + Quiet@
NIntegrate[Exp[-x]*Cos[Gamma[x]]^2, {x, 0, Log[ptp]}]}, {ptp,
0.01, 0.279, 0.001}], Frame -> True,
FrameLabel -> {"\!\(\*FractionBox[SubscriptBox[\(P\), \(t\)], \(P\)]\
\)", Rotate["\!\(\*FractionBox[\(P'\), \(P\)]\)", -Pi/2]},
FrameStyle ->
Directive[Black, Thickness[1/200], FontSize -> 20,
FontFamily -> "Times"], FrameTicks -> Automatic,
FrameTicksStyle ->
Directive[Black, Thickness[1/200], FontSize -> 16,
FontFamily -> "Times"], PlotRange -> Full, PlotStyle -> Orange,
Joined -> True]

似乎NIntegrate在这个条件下精度不高, 你可以查阅NIntegrate函数的资料来改善精度

