MIKU

今日の活動について。
骨盤をRで書こうというテーマで葛藤。

t=0:99999
x=sin(0.01*pi*t)
y=cos(0.01*pi*t)
z=0.00005*t
library(rgl)
x1=c(2.5,2.5,-2.5,-2.5,2.5,2.5,-2.5,-2.5)
y1=c(2.5,-2.5,2.5,-2.5,2.5,-2.5,2.5,-2.5)
z1=c(0,0,0,0,5,5,5,5)
plot3d(x1,y1,z1,col='white')
a1=c(cos(3/2*pi+1/3*pi),sin(3/2*pi+1/3*pi),1)
a2=c(cos(3/2*pi-1/3*pi),sin(3/2*pi-1/3*pi),1)
l1=((cos(3/2*pi+1/3*pi)-x)^2+(sin(3/2*pi+1/3*pi)-y)^2+(1-z)^2)
s1<-which(l1<(1/3)^2)
plot3d(x[s1],y[s1],z[s1],col='white',add=T)
l2=((cos(3/2*pi-1/3*pi)-x)^2+(sin(3/2*pi-1/3*pi)-y)^2+(1-z)^2)
s2<-which(l2<(1/3)^2)
plot3d(x[s2],y[s2],z[s2],col='white',add=T)
plot3d(x,y,z,add=T)

次は骨盤上口の形を変えたい。

あとphomで穴の数を特定したい。
そのためにある数の集合から特定の数を引けるようにしたい。