4.2.1.- S=1/2 con ganancias - Schrödinger

Splus={{0,1},{0,0}};
L= Splus;(*operador de Linblad de ganancias*)
Ldg=ConjugateTranspose[L]; (*Ldagger*)
H2={{(Subscript[w, 1]),0},{0,(Subscript[w, 2])}}; (*Hamiltoniano 2-dimensional*)
Ro={{(Subscript[ro, 11]),(Subscript[ro, 12])},{(Subscript[ro, 21]),(Subscript[ro, 22])}}; 

(*Calcular el Right Hand Side de la ecuación de Lindblad con ganancias para los cuatro componentes de Ro*)
RHSg11=Simplify[(-I (H2.Ro-Ro.H2)-(g/2) ((Ldg .L .Ro)[[1,1]]+(Ro.Ldg .L)[[1,1]] -2 (L. Ro.Ldg)))[[1,1]] ,g∈Reals]; 
RHSg12=Simplify[(-I (H2.Ro-Ro.H2)-(g/2) ((Ldg .L .Ro)[[1,2]]+(Ro.Ldg .L)[[1,2]] -2 (L. Ro.Ldg)))[[1,2]] ,g∈Reals];
RHSg21=Simplify[(-I (H2.Ro-Ro.H2)-(g/2) ((Ldg .L .Ro)[[2,1]]+(Ro.Ldg .L)[[2,1]] -2 (L. Ro.Ldg)))[[2,1]] ,g∈Reals];
RHSg22=Simplify[(-I (H2.Ro-Ro.H2)-(g/2) ((Ldg .L .Ro)[[2,2]]+(Ro.Ldg .L)[[2,2]] -2 (L. Ro.Ldg)))[[2,2]] ,g∈Reals];

(*Left Hand Side = Right Hand Side. Las ecuaciones que surgen para cada componente de Ro de la ecuación de Lindblad*)
Print[ Subscript["Ro'",11],"(t)=" ,RHSg11] 
Print[ Subscript["Ro'",12],"(t)=" ,RHSg12]
Print[ Subscript["Ro'",21],"(t)=" ,RHSg21]
Print[ Subscript["Ro'",22],"(t)=" ,RHSg22]

(*Resolución de las acuaciones diferenciales acopladas obtenidas en el paso anterior*)

DSolve[{ro11'[t]==g ro22[t], ro12'[t]==-1/2ro12[t](g+2I (Subscript[w, 1])-2I (Subscript[w, 2])),ro21'[t]==-1/2ro21[t](g-2I (Subscript[w, 1])+2I (Subscript[w, 2])),ro22'[t]==-g ro22[t]}, {ro11[t],ro12[t],ro21[t],ro22[t]} ,t]