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)
Note that non-ascii characters in the above source code will be escaped (such as \x9f).
download
return to the top page