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