deSolve

ロトカ=ヴォルテラの向きの分離

my.Moving.Frame.2 <- function(X,k=length(X[1,])){ n <- length(X[,1]) # No. points d <- length(X[1,]) # Dimension # inter-point distance L <- rep(0,n-1) for(i in 1:(n-1)){ L[i] <- sqrt(sum((X[i,]-X[i+1,])^2)) } diff.X <- list() #diff.X[[1]] …