function y = test_image(n,a)

y=zeros(n,n);
y(n/4:3*n/4,n/2)=1;

for j=1:n
  b=round((j-n/2)*a);
  k=mod(b:n+b-1,n)+1;
  y(j,:)=y(j,k);
end;
