lambda_min <- 0.4 
lambda_max <- 5

lambdalist <- seq(lambda_min, lambda_max, length = 1000)
Glist <- rep(0, length(lambdalist))

lalist <- c(1, .5, .125)
falist <- c(0, 1, 5)

for (i in 1:length(lambdalist))
  Glist[i] <- eulot(lambdalist[i], lalist, falist)

plot(lambdalist, Glist, type = "l", lwd = 2,
     xlab = "lambda", ylab = "G(lambda)")

# solve for zero
lambdabar <- uniroot(function(x) eulot(x, lalist, falist), interval = c(lambda_min, lambda_max))$root