N重振り子 シミュレーション~R版
こちらのシミュレーションのR版を走り書き
# calcAandB() : function to calculate matrices A and BcalcAandB <- function(Thetas, Vthetas, L, M, g){ n <- length(Thetas) A <- B <- matrix(0,n,n) for(i in 1:n){ for(j in 1:n){ k <- max(i,j) A[i,j] <- sum(M[k:n])