X1 <- function(x) 5*(log(2)-x)^4/(log(2))^5 lX1 <- function(x) 5*(2-exp(x))^4*exp(x) trange <- (0:(log(2)*100))/100 plot(trange, sapply(trange, X1), type="l", xlab="x",ylab="First order statistic", main="Comparison of X(1) and log(X(1))") points(trange, sapply(trange, lX1), type="l")