p=randn(10,1)
pold=p;
abs(roots(p))
s=(-1).^(1:10);
qq=conv(p,p.*s');qq=qq(1:2:end);qq=qq/max(abs(qq));p=qq;plot(abs(p),'-');
hold on
sleep(2);
qq=conv(p,p.*s');qq=qq(1:2:end);qq=qq/max(abs(qq));p=qq;plot(abs(p),'-');
sleep(2)
qq=conv(p,p.*s');qq=qq(1:2:end);qq=qq/max(abs(qq));p=qq;plot(abs(p),'-');
sleep(2);
qq=conv(p,p.*s');qq=qq(1:2:end);qq=qq/max(abs(qq));p=qq;plot(abs(p),'-');
sleep(2);
qq=conv(p,p.*s');qq=qq(1:2:end);qq=qq/max(abs(qq));p=qq;plot(abs(p),'-');
sleep(2);
qq=conv(p,p.*s');qq=qq(1:2:end);qq=qq/max(abs(qq));p=qq;plot(abs(p),'-');
sleep(2);
qq=conv(p,p.*s');qq=qq(1:2:end);qq=qq/max(abs(qq));p=qq;plot(abs(p),'-');
sleep(2);
qq=conv(p,p.*s');qq=qq(1:2:end);qq=qq/max(abs(qq));p=qq;plot(abs(p),'-');
sleep(2);
qq=conv(p,p.*s');qq=qq(1:2:end);qq=qq/max(abs(qq));p=qq;plot(abs(p),'-');
abs(roots(p))
