サイクロイド

メモ

wiki記事 # 2変量フーリエ n<-2 nf<-10 rs<-runif(nf) rs<-sort(rs,decreasing=TRUE) ps<-runif(nf)*2*pi #ps<-rep(0,nf) Niter<-1000 dt<-1 t<-seq(from=0,to=1,length.out=Niter)*2*pi*10 xy<-matrix(0,Niter,n) for(i in 1:nf){ #xy[,1]<-xy[,1]+rs[i]*c…