AbstractIn 1975 C. F. Chen and C. H. Hsiao established a new procedure to solve initial value problems of systems of linear differential equations with constant coefficients by Walsh polynomials approach. However, they did not deal with the analysis of the proposed numerical solution. In a previous article we study this procedure in case of one equation with the techniques that the theory of dyadic harmonic analysis provides us. In this paper we extend these results through the introduction of a new procedure to solve initial value problems of differential equations with not necessarily constant coefficients.