function h = circle()

% circle creates a circle image 
% of radius 50, with center at
% (125,125)

z = zeros(256,256);
for x = -50: 50;
   for y = -50 : 50;
     r = sqrt(x^2 + y^2);
       if r < 50 
          z(125+x,125+y)=255;
          end
    end
end
h=z;
image(h);