We present a numerical GRB shock model, where the specific characteristics of a shock wave have been taken into account. The model have been used to simulate collision between faster and slower shock wave. We assumed that slower shock wave can be represented by Gaussian function. The model was applied to simulate the peaks in the light curve for tree different GRB events in all four BATSE channels. By fitting those curves we are able to extract the basic parameters for shock waves, which we present in the table.