A method of symplectic integrations with adaptive time-steps for individual Hamiltonians in the planetary N-body problem