A new method of numerical integration for some transient response problems