Power series approximation to solutions of nonlinear systems of differential equations