Abstract
The Integral Implicit Method (IIM) (Simonneau & Crivellari 1988) solves the RT problems in stellar atmospheres model computations. A cubic spline model assures the continuity of the source function and that of the first and second derivatives at each discrete optical depth.