Given the theory developed so far, there are two ways to compute field-dependent NMR shifts, both with advantages and disadvantages. The first option is the computation of finite-field shifts by doing the rotational averaging numerically. To obtain N and D, the two integrals given by Eqs. ( ) and ( ) can be solved numerically by introducing a grid (χ i , θ j )