2011-02-20から1日間の記事一覧

2つの周期の積

書きかけ n<-4 t<-seq(from=0,to=1,length.out=1000)*2*pi*3 nt<-10000 dt<-0.001 p1<-3 p2<-2 a1<-1 a2<-0.4 delta<-runif(1) delta<-1 k<-1 x<-matrix(0,nt,n) X<-x X[,1]<-a1*cos(p1*t)*(a2*cos(p2*(t+delta))+k) X[,2]<-a1*sin(p1*t)*(a2*cos(p2*(t+delt…

トーラスになる偏微分方程式