On computation of transient response