cl <- helix(500, 1000, 36, delta=pi/6, n.smooth=1e3) ; str(cl)
require(rgl)
open3d()
spheres3d(cl$smooth, radius=1,col=2)
## ellipsoids are oriented following the helix
sizes <- equal_sizes(40, 20, 20,NROW(cl$positions))
rgl.ellipsoids(cl$positions, sizes, cl$angles, col="gold")
Run the code above in your browser using DataLab