貝の関数を作る

昨日はWikipediaやWolframの貝のパラメタ表示のメモをした 自分なりのそれも作りたい ムール貝風 phis <- seq(from=0,to=pi/8,length=100) rs <- sin(seq(from=0,to=pi,length=length(phis))) t <- seq(from=-50,to=-1,length=300) k <- 0.3 theta <- 0.1 X …

はまぐりの数学 貝のらせん(Wolfram) 貝のらせん(Wikipedia) uu <- seq(from=0,to=2*pi,length=50) vv <- seq(from=-2*pi,to=2*pi,length=100) uv <- expand.grid(uu,vv) x <- 5/4*(1-uv[,2]/(2*pi))*cos(2*uv[,2])*(1+cos(uv[,1]))+cos(2*uv[,2]) y <- 5/…