Evolution of psi fonction in NMR

Heres the problem that needs to be solved in Mathematica (it's a research
project for a course) :

A proton is influenced by a 3D magnetic field of the form b0=(0,0,b0). After
time t0 another field b1=(a*cos(wt), a*sin(wt),0) is combined with b0.

I have solved schrodinger equation for b0 and found psi0 from wich I
obtained <mu0>. Now i want to find <mu1> for b0+b1. I need to find psi1, I
am not able to solve this time dependant schrodinger's equation with
mathematica, so I obtained the evolution operator and tryed to apply it to
psi0. The problem I have is that parts of the operator have the form
exp[i*a*exp[i*b]*exp[i*c]] wich give the form
exp[i*a*exp[i*b]*exp[i*c]*d*[t-t0]] to some of the elements of psi1. To find
<mu1>, I need to conjugate and transpose psi1 and i dont know how to find
the conjugate of the elements mentioned before.

I'm wondering if I am doing something wrong or if there's a way to conjugate
and find <mu1>.

Thanks


.