My original description using lagrangian formulation based upon lagrangians
need to have a velocity dependent dissipation function added to the
lagrangian to take into account resistive as well as rective components.
Unforunately lagrangian formulation is not a tool used by many EEs.