Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes