  The formula frequently used to describe the c-axis infrared response of the
coupled electron-phonon system of bilayer cuprate superconductors and providing
important insights into the physics of these materials has been originally
obtained at the level of the phenomenological multilayer model. Here we derive
it using diagrammatic perturbation theory.
