@renormalize I looked over the write-up above somewhat carefully from post 64 and I believe (4) is in error.
In nomenclature you should be able to follow, I get the following:
## E_{1r}=E_{2r}/\tau_{12}-\rho_{12} E_{2l}/\tau_{12} ## and
##E_{1l}=\rho_{21}E_{2r}/\tau_{12}+(\rho_{21}^2/\tau_{12}+\tau_{12})E_{2l} ##.
This gives ## (M_{21}/M_{11})^2=\rho_{21}^2 ##, but that is not the energy reflection coefficient of the system. The energy reflection coefficient of the system is ## R=(E_{1l}/E_{1r})^2 ##.
The expression ## R=\rho_{21}^2 ## only holds for a single source incident on an isolated interface. The person who wrote this up seems to have blundered with their formula (4) above. I'm not infallible, but on this one I'm pretty sure I am right, and the "book" is wrong.[Edit: See below=the book did get it right]. Even though the Fresnel coefficient ## \rho_{21} ##, (which has ## R=\rho_{21}^2 ## when there is a single source incident on the isolated interface), remains a good number for multiple sources and multiple layers, this ## R ## is no longer a good number, but instead needs to be determined by what the whole system is doing.
Edit: I looked it over some more and wondered if perhaps I made a mistake, because their ## M ## is for the whole system. and I think I must retract the above, because I see how they get (4). Since the final interface has no left going wave on the right side, they set the right going wave to some arbitrary number, and the left-going wave to zero. They then have ## R=(E_{1l}/E_{1r})^2=(M_{21}/M_{11})^2 ##. My mistake.