AbstractFor the modelling of the radiation fields in stellar winds and the resulting forces new efficient algorithms are presented. In the first one, the radiative transfer equation for moving 3D media is solved analytically with the assumption that the source function is known, eg. from the solution of the NLTE rate equations. For a wind with inhomogeneities an a-posteriori error controlled finite element code is described that takes scattering explicitly into account. Finally, we present possibilities for the accurate inclusion of an arbitrary number of spectral lines in a deterministic and in a stochastic way.