3.3.- S=9/2 - Heisenberg

L110S={{0,0,0,0,0,0,0,0,0,1},{0,0,0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0,0,0},{1,0,0,0,0,0,0,0,0,0}}; (*Matriz de Gell-Mann generalizada "1,10 simétrica"*)
MatrixForm[L110S]


H10={{(Subscript[w, 1]),0,0,0,0,0,0,0,0,0},{0,(Subscript[w, 2]),0,0,0,0,0,0,0,0},{0,0,(Subscript[w, 3]),0,0,0,0,0,0,0},{0,0,0,(Subscript[w, 4]),0,0,0,0,0,0},{0,0,0,0,(Subscript[w, 5]),0,0,0,0,0},{0,0,0,0,0,(Subscript[w, 6]),0,0,0,0},{0,0,0,0,0,0,(Subscript[w, 7]),0,0,0},{0,0,0,0,0,0,0,(Subscript[w, 8]),0,0},{0,0,0,0,0,0,0,0,(Subscript[w, 9]),0},{0,0,0,0,0,0,0,0,0,(Subscript[w, 10])}}; (*Hamiltoniano 10x10*)
MatrixForm[H10]

U10={{Exp[-I (Subscript[w, 1])t],0,0,0,0,0,0,0,0,0},{0,Exp[-I (Subscript[w, 2])t],0,0,0,0,0,0,0,0},{0,0,Exp[-I (Subscript[w, 3])t],0,0,0,0,0,0,0},{0,0,0,Exp[-I (Subscript[w, 4])t],0,0,0,0,0,0},{0,0,0,0,Exp[-I (Subscript[w, 5])t],0,0,0,0,0},{0,0,0,0,0,Exp[-I (Subscript[w, 6])t],0,0,0,0},{0,0,0,0,0,0,Exp[-I (Subscript[w, 7])t],0,0,0},{0,0,0,0,0,0,0,Exp[-I (Subscript[w, 8])t],0,0},{0,0,0,0,0,0,0,0,Exp[-I (Subscript[w, 9])t],0},{0,0,0,0,0,0,0,0,0,Exp[-I (Subscript[w, 10])t]}};
U10dg={{Exp[I (Subscript[w, 1])t],0,0,0,0,0,0,0,0,0},{0,Exp[I (Subscript[w, 2])t],0,0,0,0,0,0,0,0},{0,0,Exp[I (Subscript[w, 3])t],0,0,0,0,0,0,0},{0,0,0,Exp[I (Subscript[w, 4])t],0,0,0,0,0,0},{0,0,0,0,Exp[I (Subscript[w, 5])t],0,0,0,0,0},{0,0,0,0,0,Exp[I (Subscript[w, 6])t],0,0,0,0},{0,0,0,0,0,0,Exp[I (Subscript[w, 7])t],0,0,0},{0,0,0,0,0,0,0,Exp[I (Subscript[w, 8])t],0,0},{0,0,0,0,0,0,0,0,Exp[I (Subscript[w, 9])t],0},{0,0,0,0,0,0,0,0,0,Exp[I (Subscript[w, 10])t]}};
Ro10 = 1/Z {{Exp[-b (Subscript[w, 1])],0,0,0,0,0,0,0,0,0},{0,Exp[-b (Subscript[w, 2])],0,0,0,0,0,0,0,0},{0,0,Exp[-b (Subscript[w, 3])],0,0,0,0,0,0,0},{0,0,0,Exp[-b (Subscript[w, 4])],0,0,0,0,0,0},{0,0,0,0,Exp[-b (Subscript[w, 5])],0,0,0,0,0},{0,0,0,0,0,Exp[-b (Subscript[w, 6])],0,0,0,0},{0,0,0,0,0,0,Exp[-b (Subscript[w, 7])],0,0,0},{0,0,0,0,0,0,0,Exp[-b (Subscript[w, 8])],0,0},{0,0,0,0,0,0,0,0,Exp[-b (Subscript[w, 9])],0},{0,0,0,0,0,0,0,0,0,Exp[-b (Subscript[w, 10])]}};

Print["<",Subscript[Superscript["A","S"],"1,10"],"(0)",Subscript[Superscript["A","S"],"1,10"],"(t)", "> = ",Tr[Ro10.U10dg.L110S.U10.L110S]]