Matlab - Calculo variacion puntos de cuadricula

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