variable = {Subscript[[Rho], 1, 1][t], Subscript[[Rho], 1, 2][t],Subscript[       [Rho], 1, 3][t], Subscript[[Rho], 2, 1][t], Subscript[[Rho], 2, 2][t], Subscript[[Rho], 2, 3][t],Subscript[[Rho], 3, 1][t], Subscript[[Rho], 3, 2][t], Subscript[[Rho], 3, 3][t]};
      init = {Subscript[[Rho], 1, 1][0] == 1, Subscript[[Rho], 1, 2][0] == 0, Subscript[[Rho], 1, 3][0] == 0, Subscript[[Rho], 2, 1][0] == 0, Subscript[[Rho], 2, 2][0] == 0, Subscript[[Rho], 2, 3][0] == 0, Subscript[[Rho], 3, 1][0] == 0, Subscript[[Rho], 3, 2][0] == 0, Subscript[[Rho], 3, 3][0] == 0};

     equations = {Derivative[1][Subscript[[Rho], 1, 1]][t] == [Gamma]/
 2 - [Gamma] Subscript[[Rho], 1, 1][t] -
 I (1/2 [CapitalOmega]Ra Subscript[[Rho], 1, 3][t] -
    1/2 [CapitalOmega]Ra Subscript[[Rho], 3, 1][t]) +
 1/2 [CapitalGamma] Subscript[[Rho], 3, 3][t],