light_bulbs <- 100
lifetime <- 1000
lambda_rate <- light_bulbs/lifetime
1/lambda_rate
## [1] 10
\(f_Z(z) = (1/2)e^{-\lambda|z|}\) can be re-written as \(f_Z(z) = \begin{cases} (1/2)e^{-\lambda z}, & \mbox{if } z \ge 0, \\ (1/2)e^{\lambda z}, & \mbox{if }z <0. \end{cases}\)
Since \(X_1\) and \(X_2\) have exponential density, their PDF is
\[ \begin{split} f_Z(z) &= f_{X_1+(-X_2)}(z) \\ &= \int_{-\infty}^{\infty} f_{-X_2}(z-x_1) f_{X_1}(x_1) dx_1 \\ &= \int_{-\infty}^{\infty} f_{X_2}(x_1-z) f_{X_1}(x_1) dx_1 \\ &= \int_{-\infty}^{\infty} \lambda e^{-\lambda(x_1-z)} \lambda e^{-\lambda x_1} dx_1 \\ &= \int_{-\infty}^{\infty} \lambda^2 e^{-\lambda x_1 + \lambda z} e^{-\lambda x_1} dx_1 \\ &= \int_{-\infty}^{\infty} \lambda^2 e^{\lambda z - \lambda x_1 - \lambda x_1} dx_1 \\ &= \int_{-\infty}^{\infty} \lambda^2 e^{\lambda(z-2x_1)} dx_1 \end{split} \]
mu <- 10
variance <- sqrt(100/3)
k_a <- 2/variance
upper_bound_a <- 1/k_a^2
upper_bound_a
## [1] 8.333333
k_b <- 5/variance
upper_bound_b <- 1/k_b^2
upper_bound_b
## [1] 1.333333
k_c <- 9/variance
upper_bound_c <- 1/k_c^2
upper_bound_c
## [1] 0.4115226
k_d <- 20/variance
upper_bound_d <- 1/k_d^2
upper_bound_d
## [1] 0.08333333