This paper is concerned with systems of ordinary differential equations with fuzzy parameters. Applying the Zadeh extension principle to the equations, we introduce the notions of fuzzy solutions and of componentwise fuzzy solutions. The fuzzy extension of the solution operator is shown to provide the unique fuzzy solution as well as the maximal componentwise fuzzy solution. A numerical algorithm based on monotonicity properties of membership functions is presented, together with a proof of convergence. In an interplay of interval analysis and possibility theory, these methods allow to process subjective information on the possible fluctuations of parameters in models involving ordinary differential equations. This is demonstrated in two engineering applications: a queueing model for earthwork and a model of oscillations of bell-towers.