请问第二个数值积分为啥不行
mathematica吧
全部回复
仅看楼主
level 2
2023年12月22日 12点12分 1
level 2
P[x_]:=\[Nu]^3/(N[E^((h \[Nu])/(K (Re[r NIntegrate[r^2,{r,6,x}]/.r->x])^(1/4)))]-1) x^2
Q[p_,\[Nu]_]:=NIntegrate[P[r]/.{h->1,K->5},{r,3,p},{\[Phi],0,2\[Pi]}]
2023年12月22日 12点12分 2
吧务
level 10
2023年12月23日 07点12分 3
不止这个问题。\[Nu],h和K其实都没传进P里去。
2024年01月06日 15点01分
吧务
level 15
除了3楼的链接,另请参看:
https://tieba.baidu.com/p/4886024974
一个姑且能用的版本:
Clear[P, Q]
P[x_?NumericQ, \[Nu]_, h_,
K_] := \[Nu]^3/(N[
E^((h \[Nu])/(K (Re[r NIntegrate[r^2, {r, 6, x}] /. r -> x])^(1/4)))] -
1) x^2
Q[p_, \[Nu]_] := NIntegrate[P[r, \[Nu], 1, 5], {r, 3, p}, {\[Phi], 0, 2 \[Pi]}]
2024年01月06日 15点01分 5
1