If I have an experiment with discrete probabilities like drawing black and white balls, the probability of drawing both black and white is just a multiplication of probabilities of drawing black and white. Similarly, if I want to work with "OR" type of situations, I add the probabilities.
But how does this work with continuous random variables? Say, I have the following experiment: I have a large bucket of some material and I am trying to determine the melt temperature. So I take small samples and put them on a burner and note the temperature when they melt. Let's assume that after many of such measurements, I get a gaussian with some mean value and variance.
If I then would like to say something about another bucket of the same material, I might ask, what is the probability of a sample melting at temperature between $(T,T+{\rm d}T)$. Then the probability should be $\Pr(\text{melt at }T)\approx f(T)\,{\rm d} T$ or more exactly $\int_T^{T+{\rm d}T}f(x)\,{\rm d} x$ where $f(T)$ is the gaussian.
But what if I want to ask about two samples from taken from the bucket both melting at $\{T,T+{\rm d}T\}$ or one of the samples from the bucket melting between $(T_1,T_1+{\rm d}T)$ and the other from the bucket melting between $(T_2,T_2+{\rm d}T)$? Is it just $f^2(T)\,{\rm d}x^2$ and $f(T_1)f(T_2)\,{\rm d}x^2$ respectively? Or more exactly
$$\int_T^{T+{\rm d}T}f(x)\,{\rm d} x\int_T^{T+{\rm d}T}f(x)\,{\rm d} x$$
or
$$ \int_{T_1}^{T_1+{\rm d}T} f(x){\rm d} x \int_{T_2}^{T_2+{\rm d}T}f(x)\,{\rm d} x $$
But this leads to having new probability distribution functions $f(x)f(y)$ which from what I have read are not necessarily distribution functions.
And if I want to ask about the probability of the melt temperature being between $(T_1, T_1+{\rm d}T)$ or between $(T_2,T_2+{\rm d}T)$, would it be $f(T_1)\,{\rm d}T + f(T_2)\,{\rm d}T$ or more precisely
$$\int_{T_1}^{T_1+{\rm d}T} f(x)\,{\rm d} x + \int_{T_2}^{T_2+{\rm d}T} f(x)\,{\rm d} x\text{?}$$