clear all
close all

screen = zeros(120,120);

centerx = 60;
centery = 45;

r=30;

x = zeros(1,360);
y = zeros(1,360);

a=1;
for i=0:359
    x(i+1)=r.*cos(((i-90)*pi)/180)+centerx;
    y(i+1)=r.*sin(((i-90)*pi)/180)+centery;
    a=a+1;
end


x=double(uint16(x));
y=double(uint16(y));
figure(10)
plot(x+j*y)
axis([0 90 0 75])

figure(11)
for x1=centerx:100
  		for y1=centery:90           
  			dd = ((x1-centerx)*(x1-centerx)+(y1-centery)*(y1-centery)); 
  			d = sqrt(dd);  			
  			if ((d <= 32) && (d>=30))
  				screen(x1,y1) = 1; 
				screen(centerx-(x1-centerx),y1)=1;
				screen(x1,centery-(y1-centery))=1;
				screen(centerx-(x1-centerx),centery-(y1-centery))=1;				 
            end
        end
end

image(screen'*128);