ぱらぱらめくる『曲面と多様体』3. 多様体論

• 3. 多様体論
• 付録。多次元トーラス

df<-10

Rs<-runif(df-1)
as<-runif(df-1)
bs<-runif(df-1)
T<-5
t<-seq(from=0,to=1,length=1000)*2*pi*T

X<-matrix(0,length(t),df)

X[,1]<-Rs[1]*cos(as[1]*t+bs[1])
X[,2]<-Rs[2]*sin(as[1]*t+bs[1])

for(i in 2:(df-1)){
	X[,i+1]<-Rs[i]*sin(as[i]*t+bs[i])
	X[,1:i]<-X[,1:i]+Rs[i]*cos(as[i]*t+bs[i])*X[,1:i]/sqrt(sum(X[,1:i]^2))
}

library(rgl)
plot3d(X)