n = 5;
x = 3;
y = 3;
d = 1.5;
Z = magic(n);
% Matriz de coordenadas para X e Y
X = repmat(1:n, n, 1);
Y = X';
% La distacia de cada punto a (3, 2)
distancia = sqrt((X-x).^2 + (Y-y).^2)
% La suma de todos los puntos separados d
sum(Z(distancia > d))
distancia =
2.8284 2.2361 2.0000 2.2361 2.8284
2.2361 1.4142 1.0000 1.4142 2.2361
2.0000 1.0000 0 1.0000 2.0000
2.2361 1.4142 1.0000 1.4142 2.2361
2.8284 2.2361 2.0000 2.2361 2.8284
ans =
208