Я пытаюсь численно интегрировать обобщенное нормальное распределение . Я ожидаю, что площадь под ним будет 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)