In this paper, we perform a systematic investigation of the statistics associated with the product of two independent and non-identically distributed κ-μ random variables. More specifically, we develop novel analytical formulations for many of the fundamental statistics of interest, namely, the probability density function, cumulative distribution function, and moment-generating function. Using these new results, closedform expressions are obtained for the higher order moments, amount of fading and channel quality estimation index, while analytical formulations are obtained for the outage probability, average channel capacity, average symbol error probability, and average bit error probability. These general expressions can be reduced to a number of fading scenarios, such as the double Rayleigh, double Rice, double Nakagami-m, κ-μ/Nakagami-m, and Rice/Nakagami-m, which all occur as special cases. Additionally, as a byproduct of the work performed here, formulations for the κ-μ/κ-μ composite fading model can also be deduced. To illustrate the efficacy of the novel expressions proposed here, we provide useful insights into the outage probability of a dualhop system used in body area networks, and demonstrate the suitability of the κ-μ/κ-μ composite fading for characterizing shadowed fading in device-to-device channels.