We report an implementation of the linear response function for singlet excited states for the coupled cluster models CCS, CC2 and CCSD. The implementation is based on the derivation of excited state response functions as derivatives of excited state quasienergy Lagrangians. Secular divergencies are explicitly eliminated and response equations and response functions therefore are numerical stable in the static limit. Calculations are performed for the polarizabilities of pyrimidine and s-tetrazine in their lowest singlet excited states. We find that the results for the excited state polarizabilities are sensitive to the accuracy of the excitation energies and that a qualitative correct description of the dispersion is first obtained at a correlated level.