I have response times in seconds of 3 tasks which follow Gamma distribution: F ~ $\Gamma(k_1, θ_1)$, G ~ $\Gamma(k_2, θ_2)$, H ~ $\Gamma(k_3, θ_3)$. Tasks 1 and 2 run in parallel, but task 3 runs after the first two. From this setup my goal is to calculate the probability of all tasks being executed within x seconds.
Currently my though was to do it like this: FG+H.
Calculate the shape and scale of FG like this: $$k_{new} = k_1 + k_2 - 1$$ $$θ_{new} = θ_1 + θ_2$$ Source: https://math.stackexchange.com/questions/683894/multiplying-two-gamma-distributions-over-the-same-variable
Calculate the shape and scale by adding H to FG based on this: $$k_{sum}=\frac{(\sum_iθ_ik_i)^2}{\sum_iθ_i^2k_i}$$ $$θ_{sum}=\frac{\sum_iθ_ik_i}{k_{sum}}$$ Source (Paul Harrison's answer): Generic sum of Gamma random variables
After that I would have my final shape and scale values calculate the probability. Is my approach correct?