In this paper, we consider a class of functional equation Q(λ)Y (λ) −P(λ)Z(λ) = η related to sine type functions, where the known P,Q are appropriate entire functions of exponential type. We are concerned with the existence and uniqueness of the solution (Y,Z) under certain circumstances. Furthermore, we modify the Lagrange interpolation to deal with the situation of the interpolation nodes being counted by multiplicities, which is significant to solve the above functional equation.