## Sample sizes.
N = c(5,10,20,40,80)

## A grid of lambda values.
G = seq(0.1, 5, by=0.05)

for (k in 1:length(N)) {

    ## The sample size.
    n = N[k]

    ## Replications.
    for (j in 1:10) {
        X = rexp(n)
        L = -n*log(G) - sum(X)/G

        plot(G, L, t='l', xlab='Lambda', ylab='log-likelihood', main='')
        title(main=sprintf('Sample size %d', n))

        ii = which(L >= (max(L)-1))
        lines(G[ii], L[ii], col='red')

        scan()
    }
}