Numerical solution of shock response spectra for vibration system based on Matlab is obtained under rectangular, forward & back peak saw tooth pulses. In order to acquire high accuracy of numerical solution, the paper discusses the relation satisfied between Runge-Kutta's integral time step, integral time interval and...