求助函数的数值积分解
mathematica吧
全部回复
仅看楼主
level 1
CXC0028 楼主
需要绘制P'/P为函数,Pt/P为自变量的函数图。看到函数后无法下手,在此求助各位大佬了!
2024年07月27日 07点07分 1
level 9
积分上限中有x, 积分微元中也有x, 我很疑惑
我所想的是这样的
首先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函数的资料来改善精度
2024年07月28日 14点07分 3
level 7
以Pt/P为自变量?好吧,那就把它当一个整体来看,楼上也写了,不过注意积分上限那里是P/Pt。
不过你这个函数以x为自变量看着挺好玩(),只不过我这精度也不高就是了。
2024年07月29日 10点07分 4
1