D.77 Solving the NMR equations

To solve the two coupled ordinary differential equations for the spin up and down probabilities, first get rid of the time dependence of the right-hand-side matrix by defining new variables $A$ and $B$ by

$$a=Ae^{{\rm i}\omega t/2}, \quad b=Be^{-{\rm i}\omega t/2}.$$

Then find the eigenvalues and eigenvectors of the now constant matrix. The eigenvalues can be written as $\pm{\rm i}\omega_1/f$, where $f$ is the resonance factor given in the main text. The solution is then

$$\left(\begin{array}{c}A\\ B\end{array}\right) = C_1 \v... ...{{\rm i}\omega_1t/f} + C_2 \vec v_2 e^{-{\rm i}\omega_1t/f}$$

where $\vec v_1$ and $\vec v_2$ are the eigenvectors. To find the constants $C_1$ and $C_2$, apply the initial conditions $A(0)=a(0)=a_0$ and $B(0)=b(0)=b_0$ and clean up as well as possible, using the definition of the resonance factor and the Euler formula.

It's a mess.