Z = -3:6/512:3-(6/512);
Z = Z(ones(1,length(Z)),:);
Z = complex(Z,Z');
Z = (Z+1)./(Z-1);
f = figure();
frame = 0;
colormap('gray')
axis square
axis off
im = image();
while(true)
    g1 = mod(log(abs(Z))*4,2)<1;
    g2 = mod(angle(Z)+(frame-1)*pi/180, pi/6)<(pi/12);
    frame = frame + 10;
    im.CData = xor(g1, g2);
end