function lagplot(z,lag)
[n m] = size(z);
plot(z(1:(n-lag)),z((lag+1):n),'.')