function im = CS4640_create_square(M,N,x_min,y_min,s)
%

im = zeros(M,N);

im(x_min:x_min+s-1,y_min:y_min+s-1) = 1;
