% TEST_APD_PDF.M     Version I.K. 04/09/02
% Comments: 
%   a standard APD random variable X has pdf :
%   p(alpha,lda,x) = [A^(1/lda)/gamma(1+1/lda)]*...
%                   *exp[-(A/alpha^lda)*abs(x)^lda],        if x <= 0,
%                   *exp[-(A/(1-alpha)^lda)*abs(x)^lda],    if x > 0,
%   where A = (2*alpha^lda*(1-alpha)^lda)/(alpha^lda+(1-alpha)^lda).

alpha= 0.5;
lda = 2;
Nprec = 100;
x = (-10:(10/Nprec):10)';
alphax = alpha*ones(size(x));
ldax = lda*ones(size(x));
d = APDpdf(alphax,ldax,x);

plot(x,d)