The partial fraction decomposition method is explored to establish several interesting trigonometric function identities, which may have applications to the evaluation of classical multiple hypergeometric series, trigonometric approximation and interpolation. 1. Outline and introduction Recently, in an attempt to prove, through the Cauchy residue method, Dougall’s theorem (Dougall [6, 1907], se...