  By introducing an integration factor to the differential one-form of contact
dynamics, equations of motion are derived variationally, and contact Poisson
bracket and contact Lagrangian are formulated. Discrete symplectic integrator,
named hybrid leap-frog method, is found to be a numerical solver of equations
of motion for contact dissipative systems.
