こちらの続き #k 次元立方格子 k<-2 P<-4 A<-expand.grid(rep(list(1:P),k)) N<-length(A[,1]) M2<-matrix(0,N,N) for(i in 1:N){ for(j in 1:N){ tmp<-A[i,]-A[j,] if(sum(tmp^2)==1){ M2[i,j]<-1 } } } M2 image(M2) # ラプラシアンの作成 d<-apply(M2,1,s…
引用をストックしました
引用するにはまずログインしてください
引用をストックできませんでした。再度お試しください
限定公開記事のため引用できません。