Two-stage waveform relaxation method for the initial value problems with non-constant coefficients