Generation of a synthetic seismogram

How can we generate a synthetic seismogram?
Procedure to create a synthetic seismogram is as follows:
- Multiply the velocity (calculated from the sonic log) and density logs to generate an acoustic impedance (AI) log. When a density log is not available, the densities can be calculated from the velocities with Gardner's rule: the density is proportional to the ¼ power of the P-wave velocity.

- Calculate from the AI log the reflection coefficients (using Zoeppritz' equation)

- Determine the wavelet from the seismic data

- Convolve the wavelet with the reflection coefficient trace to generate the synthetic trace