When I read an article by Garbe, where have an eigenmatrix $$B M_i + M_i B^T=\lambda_iM_i$$ and have the following eigenmatrices and eigenvalues:
$$ M_0=\begin{bmatrix} 0 & 1 \\ -1 & 0\end{bmatrix},M_1=\begin{bmatrix} \frac{1}{\sqrt{Xg^2-1}} & 0 \\ 0 & -\sqrt{Xg^2-1}\end{bmatrix}, M_{\pm}=\begin{bmatrix} \frac{1}{\sqrt{Xg^2-1}} & \pm1 \\ \pm1 & \sqrt{Xg^2-1}\end{bmatrix}\\ \lambda_0=\lambda_1=0, \lambda_{\pm}=\pm2\omega_0\sqrt{Xg^2-1}$$
I tried to solve the matrix equation, hoping to get the same result as the article, but I didn't have any beneficial results.
Below is my code for the question
B=B = {{0, w0}, {w0(X g^2-1), 0}} M=Array[mM = Array[m, {2, 2}] \[Lambda]=λ = {{\[Lambda]1λ1,\[Lambda]2 λ2}, {\[Lambda]3λ3,\[Lambda]4 λ4}} sol=Solve[Bsol = Solve[B.M+MM + M.B\[Transpose]==\[Lambda]B\[Transpose] == λ M, Flatten[M]] It's not producing the results I want.
Any help is much appreciated.