# Show 12 random vertices in various ways. Force 4-tuples into
# planes so that quads3d works.
x <- rnorm(12)
y <- rnorm(12)
z <- rnorm(12)
for (i in c(4,8,12))
z[i] <- predict(lm(z ~ x+y,
data=data.frame(x=x[(i-3):i], y=y[(i-3):i], z=z[(i-3):i])))[4]
rgl.open()
points3d(x,y,z)
lines3d(x-3,y+3,z)
segments3d(x-6,y+6,z)
triangles3d(x-9,y+9,z,col='red')
quads3d(x-12,y+12,z,col='green')
text3d(x-15,y+15,z,1:12)
# Add labels
text3d(-3*(0:5)+3, 3*(0:5), 0,
c('points3d','lines3d','segments3d',
'triangles3d', 'quads3d','text3d'),
adj = 0)
rgl.bringtotop()
Run the code above in your browser using DataLab