cubique simple : c=0
cubique à  faces centrées : c=1

rayon de chaque sphère :

Sphères bleues et rouges : nb=2
Seulement les sphères bleues : nb=1

A: (0, 0, 0)
B: (1, 0, 0)
C = Point(Circle[B, Distance[A, B], Segment[A, B]])
a, D, E, F, G, H, faceABCD, faceADHE, faceABFE, faceBCGF, faceCDHG, faceEFGH, arêteAD, arêteAB, arêteBC, arêteCD, arêteAE, arêteDH, arêteEH, arêteEF, arêteBF, arêteFG, arêteCG, arêteGH = Cube(A, B, C)
I = Midpoint(B, D)
J = Midpoint(B, G)
K = Midpoint(G, E)
L = Midpoint(H, C)
M = Midpoint(A, H)
N = Midpoint(A, F)
R = 0.14 (slider from 0 to 0.6, step 0.01)
b = Sphere(D, R)
c_1 = Sphere(B, R)
d = Sphere(G, R)
e = Sphere(F, R)
f = Sphere(H, R)
g = Sphere(E, R)
h = Sphere(C, R)
i = Sphere(A, R)
j = Sphere(L, R)
k = Sphere(K, R)
l = Sphere(N, R)
m = Sphere(J, R)
n_1 = Sphere(I, R)
o = Sphere(M, R)
c = 0 (slider from 0 to 1, step 1)
nb = 2 (slider from 1 to 2, step 1)