clear all
clf
set(0,'defaultaxesfontsize',22);
p=0.5; N=20;
x=(p^2-1)*p.^([0:N])
x(1)=x(1)+1;
x=(1/p)*x;
rxx=conv(x,fliplr(x));
stem([-N:N],rxx)


