clf;
theta = atan(5./12.);
minx = 0;
maxx = 13;
miny = 0;
maxy = 13;

for ix = minx:maxx,
 for iy = miny:maxy,
 AX(ix-minx+1,iy-miny+1) = ix;
 AY(ix-minx+1,iy-miny+1) = iy;
 BX(ix-minx+1,iy-miny+1) = AX(ix-minx+1,iy-miny+1) * cos(theta) -  AY(ix-minx+1,iy-miny+1) * sin(theta);
 BY(ix-minx+1,iy-miny+1) = AY(ix-minx+1,iy-miny+1) * cos(theta) +  AX(ix-minx+1,iy-miny+1) * sin(theta);
 end;
end;

hold on;
axis('equal'); 
for ix = minx:maxx,
 for iy = miny:maxy,
  plot (AX(ix-minx+1,iy-miny+1),AY(ix-minx+1,iy-miny+1),'o'); 
  plot (BX(ix-minx+1,iy-miny+1),BY(ix-minx+1,iy-miny+1),'x'); 
 end;
end;