負の数の平方根なども取りたいので複素数で考える # 複素数にも対応している my.mean <- function(x,p){ (sum(x^p)/length(x))^(1/p) } # 実数列は複素数列に直す x <- rnorm(1000,0,1) x <- complex(real=x,imaginary=rep(0,length(x))) ps <- seq(from=-10…