最近想要复现某篇论文里面关于计算广义高斯分布差分熵的公式：
论文中说式（5）是用mathematica计算得出的，我下了一个mathematica试着自己敲了一下，得到如下结果：
代码如下：
A = s*(Gamma[1/p]/Gamma[3/p])^(1/2)；
Pg = 1/(2*Gamma[1+1/p]*A)*Exp[-Abs[y/A]^p]；
H = Integrate[-Pg*Log[Pg],y]；
我发现好像我写的mathematica并不能计算带绝对值的指数项，请问论文中的式（5）是如何做到的呢？
请各路大神指教，太感谢啦！！！！！

因为你的代码里面没有考虑到p>0和s>0的限制条件
加上假定条件就行了
Assuming[{s > 0, p > 0},Integrate[-Pg*Log[Pg], {y, -\[Infinity], \[Infinity]}]]
你要问我是如何发现的
我用数值积分对比了一下论文的结论猜的
计算比较缓慢，耐心等待