The task to numerically compute reliable and accurate emergent Stokes profiles is of great relevance in solar physics. Aiming at facilitating the comprehension of the advantages and drawbacks of the different formal solvers, this work presents a reference paradigm for their characterization based on the concepts of order of accuracy, stability, and computational cost.