リーマン面

• 説明は長くなり、この日の記事の組みひもから遠くなるので、この記事は別ブログのこちらへ引っ越し
• 

theta<-seq(from=0,to=2*pi,length.out=100)
r<-seq(from=0,to=1,length.out=10)
xy<-expand.grid(r,theta)
plot(xy[,1],xy[,2])
x<-xy[,1]*cos(xy[,2])
y<-xy[,1]*sin(xy[,2])
plot(x,y)

z1<-sqrt(xy[,1])*cos(xy[,2]/2)
z2<-sqrt(xy[,1])*cos(xy[,2]/2+pi)
xx<-c(x,x)
yy<-c(y,y)
zz<-c(z1,z2)

plot3d(xx,yy,zz, col=rainbow(1000))


X <- par3d("userMatrix") 
movie3d( par3dinterp( userMatrix=list(X,rotate3d(X, pi/2, 1, 0, 0),rotate3d(X, pi/2, 0, 1, 0) )), duration=5 ,movie="test2",dir=".")

• 

################

theta<-seq(from=0,to=2*pi,length.out=100)
r<-seq(from=0,to=1,length.out=10)
xy<-expand.grid(r,theta)
plot(xy[,1],xy[,2])
x<-xy[,1]*cos(xy[,2])
y<-xy[,1]*sin(xy[,2])
plot(x,y)

a<-complex(real=0.1,imaginary=0.2)*4
b<-complex(real=0.2,imaginary=0.3)
c<-complex(real=0.4,imaginary=0.5)

d<-complex(real=x,imaginary=y)
e<-(d-a)*(d-b)*(d-c)
plot(Re(d),Im(d))


z1<-sqrt(e)
z2<--sqrt(e)
xx<-c(Re(d),Re(d))
yy<-c(Im(d),Im(d))
zz<-c(Re(z1),Re(z2))

plot3d(xx,yy,zz, col=rainbow(1000))


X <- par3d("userMatrix") 
movie3d( par3dinterp( userMatrix=list(X,rotate3d(X, pi/2, 1, 0, 0),rotate3d(X, pi/2, 0, 1, 0) )), duration=5 ,movie="test2",dir=".")