N = 8;
rng = -N:(N+2);
sz = size(rng);
xmesh = ones(sz)'*rng;
ymesh = xmesh';
M = sz(2)-3;
xmesh = xmesh(2:M,:);
ymesh = ymesh(2:M,:);

xg = xs + u(1)*xmesh + v(1)*ymesh;
yg = ys + u(2)*xmesh + v(2)*ymesh;
sz = size(xg);
hold on
for j=1:sz(1)
    plot(xg(j,:),yg(j,:),'y','LineWidth',2);
end
for j=1:sz(2)
    plot(xg(:,j),yg(:,j),'y','LineWidth',2);
end
hold off