2012-05-01から1ヶ月間の記事一覧
Npt<-10000 t<-(0:Npt)/1000*pi C1<-1 C2<-0 w<-1 cc<-C1*cos(w*t+C2) xframe<-c(-1,1) yframe<-c(-1,1) for(i in 1:(Npt-1)){ plot(sin(cc),cos(cc),xlim=xframe,ylim=yframe,type="l") par(new=TRUE) plot(sin(cc[i]),cos(cc[i]),col=2,xlim=xframe,ylim=y…
Npt<-10000 t<-(0:Npt)/1000*pi C1<-1 C2<-0 w<-1 cc<-C1*cos(w*t+C2) xframe<-c(-1,1) yframe<-c(-1,1) for(i in 1:(Npt-1)){ plot(sin(cc),cos(cc),xlim=xframe,ylim=yframe,type="l") par(new=TRUE) plot(sin(cc[i]),cos(cc[i]),col=2,xlim=xframe,ylim=y…
