Trapz给出了负面结果：为什么？

Question

我正在尝试以数字方式整合Generalized Normal Distribution。我希望它下面的区域是1.我得到的是-1。为什么会这样？

mu=0;
alpha=0.5;
beta=2;

x=-2:0.01:2;
densityGG = beta/(2*alpha*gamma(1/beta)) * exp(-power(abs(x-mu)/alpha,beta));

plot(x,densityGG);
trapz(densityGG,x)

Answer 1

该功能实际上是

trapz(X, Y)

所以试试

trapz(x, densityGG)

它应该解决问题。

Answer 2

阅读Trapz的帮助文件......一切都在那里

你需要写：     trapz（X，densityGG）