psi <- (1+sqrt(5))/2
psi.angle <- (2*pi)/(1+psi)
t <- 0:5
angles <- psi.angle * t
plot(t,angles)
x <- cos(angles)
y <- sin(angles)
plot(x,y)
plot(x,y,type="b")
t <- 0:100
angles <- psi.angle * t
plot(t,angles)
x <- cos(angles)
y <- sin(angles)
plot(x,y)
plot(x,y,type="b")
t <- 0:5
r <- t
angles <- psi.angle * t
plot(t,angles)
x <- cos(angles)
y <- sin(angles)
x2 <- x * r
y2 <- y * r
plot(x2,y2,type="b")
t <- 0:100
r <- t
angles <- psi.angle * t
plot(t,angles)
x <- cos(angles)
y <- sin(angles)
x2 <- x * r
y2 <- y * r
plot(x2,y2,type="b")
library(rgl)
k <- 1
z <- t * k
plot3d(x2,y2,z)