box(pos=vec(0,0,0),length=1,height=2,width=3,color=color.red,opacity=0.5)
def axes(Lax):
arrow(pos=vec(-Lax/2,0,0),axis=vec(Lax,0,0),color=color.yellow,fixedwidth=True,shaftwidth=0.05)
arrow(pos=vec(0,-Lax/2,0),axis=vec(0,Lax,0),color=color.green,fixedwidth=True,shaftwidth=0.05)
arrow(pos=vec(0,0,-Lax/2),axis=vec(0,0,Lax),color=color.blue,fixedwidth=True,shaftwidth=0.05)
axes(4)
g = gdots(color=color.magenta)
npts=800
phig=(3-sqrt(5))*pi
for i in range(0,npts):
r=5*sqrt(i/npts)
x=r*cos(i*phig)
y=r*sin(i*phig)
g.plot(x,y)