I don't understand why if you multiply a scaling matrix with rotation matrix that the resulting matrix, when applied to a shape like an ellipse, only gets scaled and does not get rotated.
$$\begin{bmatrix} 3.7 & 0 \\ 0 & 2.1 \end{bmatrix} \times \begin{bmatrix} 0.9510 & -.3090 \\ 0.3090 & 0.9510 \end{bmatrix}$$
If you do the matrix multiplication the opposite way it does do both when the resulting matrix is applied to an ellipse.
The only way you could be "right" is for a rotation with $\theta=n\pi/2$ ($n \in \mathbb{Z})$. If $n$ is odd you need a scaling that "reverses" the shape of the ellipse, in such a way that the major axis becomes the minor axis, and the minor axis becomes the major axis, with the same proportions as the original ellipse.