Through a detailed analysis of the properties of a system of differential equations, bounds are given for the error affecting the final result of a numerical integration. These bounds appear to be narrower than those obtained with other methods. The key procedure is to consider carefully the linear part of the system and to bound it taking account of all possible errors. No very significant restriction is made on the system.