  Using the Sklyanin-Kharchev-Lebedev method of Separation of Variables adapted
to the cyclic Baxter--Bazhanov--Stroganov or $\tau^{(2)}$-model, we derive
factorized formulae for general finite-size Ising model spin matrix elements,
proving a recent conjecture by Bugrij and Lisovyy.
