Quantum entanglement of a system of two coupled quantum harmonic oscillators with a Hamiltonian H[over ̂]=1/2(1/m_{1}p[over ̂]_{1}^{2}+1/m_{2}p[over ̂]_{2}^{2}+Ax_{1}^{2}+Bx_{2}^{2}+Cx_{1}x_{2}) can be found in many applications of quantum and nonlinear physics, molecular chemistry, and biophysics. Despite this, the quantum entanglement of such a system is still a problem under study. This is primarily due to the fact that the system is multiparametric and the quantum entanglement of such a system is not defined in a simple analytical form. This paper solves this problem and shows that quantum entanglement depends on only one parameter that has a simple physical meaning: the reflection coefficient R∈(0,1). The reflection coefficient R has a simple analytical form and includes all the parameters of the system under consideration. It is shown that for certain values of the coefficient R, the quantum entanglement can be large. The developed theory can be used not only for calculating quantum entanglement, but also for many other applications in physics, chemistry, and biophysics, where coupled harmonic oscillators are considered.