Clear[p, t, mq, ep, gh]
ff[p_] := 1/(Exp[ep[p]/t] + 1);
D[ff[p], ep[p]]
(*将函数ff[p]对ep[p]求导，可以实现*)
结果为
-(E^((ep[p]/t))/((1 + E^(ep[p]/t))^2 t))；
想实现的是，将D[ff[p], ep[p]]中的p进行积分：
Clear[p, t, mq, ep, gh]
ff[p_] := 1/(Exp[ep[p]/t] + 1);
D[ff[p], ep[p]]
Integrate[D[ff[p], ep[p]], {p, 0, 100}]
这里尚未给出ep[p]的表达式和t的数值。
下面给出ep[p]的表达式和t的数值，
结果就完全不对了：（
Clear[p, t, mq, ep, gh]
ff[p_] := 1/(Exp[ep[p]/t] + 1);
t = 50;
ep[p_] := Sqrt[p^2 + 60^2];
D[ff[p], ep[p]]
NIntegrate[D[ff[p], ep[p]], {p, 0, 100}]
想实现的是，将D[ff[p], ep[p]]中的p进行积分：
直接将D[ff[p], ep[p]]的结果放在NIntegrate命令的被积函数处可以得到结果，但我想用上面那种形式，而不是直接代入D[ff[p], ep[p]]的结果。
请教高手啊！！！ |
回复此楼
+1/177
