This works for the figure on the left:
xy = [1:14]*200;
[X, Y] = meshgrid(xy);
xyo = xy([4 11]);
figure(1)
plot(X, Y, '.k')
hold on
plot(X(:,[4 11]), Y(:,[4 11]), 'ok')
plot(X([4 11],:), Y([4 11],:), 'ok')
plot(X([4 11],[4 11]), Y([4 11],[4 11]), 'sr')
hold off
axis square
This will get you started. Make the necessary changes to get the one on the right. See the documentation for the legend function to get the legend textbox.
