I'm new to Mathematica and want to solve a system of equations. First I have to do some matrix calculations and then I want to solve the set of equations for r11, r12, ....., r33 (which stands for the density matrix elements). I always get the error message:
"Equations may not give solutions for all "solve" variables."
and I get the solutions for r11 up to r32 in dependence of r33, so mathematica doesn't solve r33. Maybe there is a problem with the parameters? I already searched the web but couldn't find any help. Although I'm pretty sure the calculation is possible, since it is a common excercise in atomic physics. I post my mathematica code below and as a file attached, I hope someone can help me :)
ClearAll;
H={{0,0,OmP},{0,2*d,OmC},{OmP,OmC,2*d}};
r={{r11,r12,r13},{r21,r22,r23},{r31,r32,r33}};
Kommutator=(i/2)*(H.r-r.H);
G={{G1*r33,0,-(1/2)*r13*(G1+G2)},{0,G2*r33,-(1/2)*r23*(
G1+G2)},{-(1/2)*r31*(G1+G2),-(1/2)*r32*(G1+G2),-r33*(G1+G2)}};
Nul={{0,0,0},{0,0,0},{0,0,0}};
Solve[{Kommutator+G==Nul},{r11,r12,r13,r21,r22,r23,r31,r32,r33}]
Attachments: