  Analytical solutions to nonlinear differential equations -- where they exist
at all -- can often be very difficult to find. For example, Duffing's equation
for a system with cubic stiffness requires the use of elliptic functions in the
exact solution. A system with general polynomial stiffness would be even more
difficult to solve analytically, if such a solution was even to exist.
Perturbation and series solutions are possible, but become increasingly
demanding as the order of solution increases. This paper aims to revisit,
present and discuss a geometric/algebraic method of determining system response
which lends itself to automation. The method, originally due to Fliess and
co-workers, makes use of the generating series and shuffle product,
mathematical ideas founded in differential geometry and abstract algebra. A
family of nonlinear differential equations with polynomial stiffness is
considered; the process of manipulating a series expansion into the generating
series follows and is shown to provide a recursive schematic, which is amenable
to computer algebra. The inverse Laplace-Borel transform is then applied to
derive a time-domain response. New solutions are presented for systems with
general polynomial stiffness, both for deterministic and Gaussian white-noise
excitation
