Summary: The modified regularized long wave equation is solved numerically by the \(B\)-spline finite element method. Three polynomial invariant conditions are reported. Solitary wave development, motion, and interaction are studied through computer simulation. Earlier estimates of the magnitude of the coefficient of inelasticity are confirmed. The development of a Maxwellian initial condition into solitary waves is shown to mimic that of the RLW equation.