function [output] = circle(input,center,radius,intensity)

[m n] = size(input);
output = input;

for i=1:m
  for j=1:n

      if( (i-center(1))^2 + (j-center(2))^2 <= radius^2 )
          output(i,j) = intensity;
      end
  end
end

