Time integration method