4.1.1.- S=1/2 con pérdidas - Schrödinger

Sminus={{0,0},{1,0}};
L= Sminus;(*operador de Linblad de pérdidas*)
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 pérdidas para los cuatro componentes de Ro*)
RHSp11=Simplify[(-I (H2.Ro-Ro.H2)-(p/2) ((Ldg .L .Ro)[[1,1]]+(Ro.Ldg .L)[[1,1]] -2 (L. Ro.Ldg)))[[1,1]] ,p∈Reals]; 
RHSp12=Simplify[(-I (H2.Ro-Ro.H2)-(p/2) ((Ldg .L .Ro)[[1,2]]+(Ro.Ldg .L)[[1,2]] -2 (L. Ro.Ldg)))[[1,2]] ,p∈Reals];
RHSp21=Simplify[(-I (H2.Ro-Ro.H2)-(p/2) ((Ldg .L .Ro)[[2,1]]+(Ro.Ldg .L)[[2,1]] -2 (L. Ro.Ldg)))[[2,1]] ,p∈Reals];
RHSp22=Simplify[(-I (H2.Ro-Ro.H2)-(p/2) ((Ldg .L .Ro)[[2,2]]+(Ro.Ldg .L)[[2,2]] -2 (L. Ro.Ldg)))[[2,2]] ,p∈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)=" ,RHSp11] 
Print[ Subscript["Ro'",12],"(t)=" ,RHSp12]
Print[ Subscript["Ro'",21],"(t)=" ,RHSp21]
Print[ Subscript["Ro'",22],"(t)=" ,RHSp22]

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

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