The biconditional is $(a_n+b_n\Delta)(c_n+d_n\Delta)y_n = (c_n+d_n\Delta)(a_n+b_n\Delta)y_n$ iff $d_n/b_n = \Delta c_n / \Delta a_n = A$, where $A$ is a constant.
I showed the if part, but I am finding it difficult to prove the only if part; I got so far to the following:
$$(a_nd_n - c_nb_n)\Delta y_n = -b_n[c_{n+1}\Delta y_n +y_n \Delta c_n + d_{n+1}\Delta\Delta y_n + (\Delta y_n)\Delta d_n]+d_n[a_{n+1}\Delta y_n +(\Delta a_n)y_n+b_{n+1}\Delta\Delta y_n +(\Delta b_n)\Delta y_n$$
how from the last identity can I show that necessarily there exists a constant $A$ s.t. $d_n/b_n = \Delta c_n / \Delta a_n = A$?
Thanks.
Edit: just in case you don't understnad: $\Delta y_n = y_{n+1} - y_n$.
Here is a slightly different approach which might be helpful. We start with a few notational conventions which help to simplify calculations.
Now we are ready to go.
Comment:
In (1) we use $\Delta=\operatorname{E}-1$.
In (2) we observe $(a-b)(c-d)$ cancels out.
We observe from (4) the operators commute if and only if the coefficient functions of $y_{n+1}$ and $y_{n+2}$ are equal on both sides. Comparing the coefficients of $y_{n+2}$ we see \begin{align*} b_nd_{n+1}=d_nb_{n+1}\\ \end{align*} from which condition $\frac{d_n}{b_n}=\frac{d_{n+1}}{b_{n+1}}=\text{const.}$ follows.