x = 1;
y = 1;
l = 2;
w = 7;
ox = 10;
oy = 10;

rectangle('Position',[x,y,l,w],'LineWidth',4)
axis([0,10,0,10]);

hold on;

rectangle('Position', [3,7,2,3], 'LineWidth', 3);

for i = 0:.1:ox
    r = 0 + 10 * rand(1,1);
    if ( r>= x && r<= w+x && i>=y && i<=l+y)
        v = 1 + (3-1) * rand(1,1);
        plot(v,r,'*');
        plot_x = v;
        plot_y = r;
        break;
    end
end

k = 1;

for j = 1:2*w
    for i = 1:2*l
    rectangle('Position',[x+(i-1)/2,y + (j-1)/2,l/(2*l),w/(2*w)], 'LineWidth', 2);
end
end

for j = 1:2*w
    for i = 1:2*l
        p(k) = x + (i-1)/2;
        k = k + 1;
        p(k) = x + (i/2);
        k = k + 1;
        p(k) = x + (i-1)/2;
        k = k + 1;
        p(k) = y + (j/2);
        k = k + 1;
        %fprintf('plot_x = %f, plot_y = %f\n',plot_x,plot_y);
    end
end

grid(gca,'minor');
hold on;