Following the work of M. Fliess [1,2] and earlier work by Hespel [3], the authors show how to write a MACSYMA program that, given an analytic dynamical system q^′(t)= A₀(q)+&Sgr;u_i(t)A_i(q), y(t)=f(q), generates a bilinear system q&barbelow;^′(t)=[M₀+&Sgr;u_i(t)M_i(t)]q&barbelow;(t), y(t)=&mgr;q&barbelow;(t) that approximates the generating series of the original system up to a prescribed order k (and therefore locally describes the solutions of the given system) and is amenable to symbolic computer solution. The problem and the solution are interesting. Unfortunately, the paper contains almost no exposition of the underlying mathematics or explanation of the notation. Hence, to understand the paper it is necessary to read the references.