logp <-function(x) log(max(x,exp(1))) x <- -2.4 delta<- 0.1 (abs(x))^(4+delta)/(logp(abs(x)))^2 x<- seq(from=-2.8,to=-2,length=3) x s<-c() j<-1 for (i in x){ s[j]<-(abs(i))^(4+delta)/(logp(abs(i)))^2 j=j+1 } s