Factorization of linear differential operators