How to create a sphere filled with equidistant points inside it?

I want to create a sphere with equidistant points inside it. The sphere should be placed essentially in a matrix, so that the sphere can be detected. Wherever there is a sphere.. the code should output 1 and viceversa.

Answer 1

[Walter Roberson]

xc = 10; yc = 15; zc = 12; R = 8;
[X, Y, Z] = meshgrid(linspace(0,20,100));
M = double((X-xc).^2 + (Y-yc).^2 + (Z-zc).^2 <= R.^2);
isosurface(X, Y, Z, M, 0)
axis equal

[Rushil Asthana]

Hey, thanks bro. But this is not exactly what i want. The sphere should be filled with equidistant inside it not on the surface. Could you please help me in this?