The topic of this paper is formal solutions of linear differential equations with formal power series coefficients. The method proposed for computing these solutions is based on factorization of differential operators. The notion of exponential parts is introduced to give a description of factorization properties and to characterize the formal solutions. The algorithms will be described and their implementation is available.