Residual.m
Sets up interaction tensors under partial ordering in a liquid crystal with the user-supplied order matrix. All adjustable parameters are set during the call to create.m function. Syntax:
spin_system=residual(spin_system)
Note: this function is only applicable to high-field NMR.
Note: the function overwrites the interaction tensors supplied by the user. Relaxation superoperator, if required, must be computed before this function is called.