  We implement continued-fraction techniques to solve exactly quantum master
equations for a spin with arbitrary S coupled to a (bosonic) thermal bath. The
full spin density matrix is obtained, so that along with relaxation and
thermoactivation, coherent dynamics is included (precession, tunnel, etc.). The
method is applied to study isotropic spins and spins in a bistable anisotropy
potential (superparamagnets). We present examples of static response, the
dynamical susceptibility including the contribution of the different relaxation
modes, and of spin resonance in transverse fields.
