function a = atlok(n, m)
    a = zeros(n, m);
    x = ((1 : n) * (n + 1) - n)' * ones(1, m - n + 1);
    y = ones(n, 1) * ((0 : m - n) * n);
    a(x + y) = 1 : numel(x);
end