# Define the rate parameter lambda <- 1.4 # Define the bounds lower <- 0.8 upper <- 3.4 # Calculate the cumulative probabilities using pexp() prob <- pexp(upper, rate = lambda) - pexp(lower, rate = lambda) # Round the result to 3 decimal places round(prob, 3)