r=-50:50
s=expand.grid(x=r,y=r,z=r)
a=1
for(d in 50:0){s=subset(s,x^2+y^2+z^2<=d^2);a[d+1]=length(s$x)}
write(a,1,1)