Here, we do the calculations of Section 1.9 using matrices.

s={{a,b},{c,d}}

MatrixForm[s]

s1={{A,B/2},{B/2,C}}

MatrixForm[s1]

s.s1.Transpose[s]

Simplify[%]

MatrixForm[%]