
n = 50;

grid = linspace(-2,2,n);


printf(" %i", [n,n,n]);
printf("\n");
      
for z = grid
  for y = grid
    for x = grid
      r = sqrt(x^2 + y^2 + z^2);
      printf(" %g\n", x*(1+0.5/r^3));
    endfor
  endfor
endfor
