function f=stablepdf(x,a)
%symmetric stable law with parameter ``a''  see Bratley et al p. 154
%the pdf is sum(f).
k=1:100;
i=2*k-1;
f=x.^(i-1);
f=(-1).^(k-1).*f.*gamma(i/a)/gamma(i);



