こちら のつづき リミットサイクル リミットサイクルになる条件 library(rgl) A<-0.5 B<-2 X0<-A Y0<-B/A dt<-0.1 T<-30 Nt<-T/dt M<-matrix(0,3,Nt+1) N<-c() f1<-function(X){ y<- A+X[1]^2*X[2]-B*X[1]-X[1] return(y) } f2<-function(X){ y<- B*X[1]-X[1…
引用をストックしました
引用するにはまずログインしてください
引用をストックできませんでした。再度お試しください
限定公開記事のため引用できません。