Static equations for thin inextensible elastic rods, or elastica as they are sometimes called, have been studied since before the time of Euler. In this paper, we examine how to model the dynamic behavior of elastica. We present a fairly high speed, robust numerical scheme that uses (i) a space discretization that uses cubic splines, and (ii) a time discretization that preserves a discrete vers...