In the previous chapter’s review of studies on the Lamb’s problem, two methods were mentioned. One of these, originally used in Lamb (1904), involves first analyzing the steady-state case, then using Fourier transforms to synthesize the solution for any time function, resulting in a double integral for the final solution. In the first half of the 20th century, handling this double integral was the main challenge, forcing researchers to use various asymptotic analysis techniques to obtain approximate solutions. However, this situation naturally improved with the introduction of computers into the field of seismology. If our goal is simply to compute this double integral, rather than analyze the wave characteristics based on it, numerical calculations and Fast Fourier Transforms (FFT) can achieve this objective.