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

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

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

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

