num_lightbulb <- 100
lifetime_lightbulb <- 1000
burnout <- lifetime_lightbulb/num_lightbulb
burnout
## [1] 10
Z=X1−X2
X2=X1−Z
fZ(z)=∫∞∞fX1(X1)fX2(X1−Z)dx1
fZ(z)=∫∞Zλe−λX1λe−λ(X1−z)dx1
fZ(z)=∫∞Zλ2e−2λX1+λZdx1
fZ(z)=∫∞Z(1/2)λe−λZ
#Solution
upper_bound_func <- function(x){
result <- (100 / (3 * x**2))
return (result)
}
prob_a <- upper_bound_func(2)
prob_b <- upper_bound_func(5)
prob_c <- upper_bound_func(9)
prob_d <- upper_bound_func(20)
prob_a
## [1] 8.333333
prob_b
## [1] 1.333333
prob_c
## [1] 0.4115226
prob_d
## [1] 0.08333333