Is there a canonical/performant way to reduce arrays/matrices by removing the border values?

Question

A motivating issue, implemented in Matlab:

N = 1000;
R = zeros(2*N);

for i=0:N-1
    R = R(2:end-1, 2:end-1);
end

For this code timeit() giv