Роботу присвячено розробленню методу, алгоритму та програми для розрахунку сейсмічних хвиль на вільній поверхні вертикально-неоднорідного півпростору, збурених локальним джерелом. Для цього використано отримані раніше співвідношення для скалярних потенціалів прямих P-, SV- та SH-хвиль від сили у формі інтегралів Бесселя-Мелліна. Це дало змогу вивести модифікованим матричним методом Томсона-Хаскела формули для повного хвильового поля у вертикально-неоднорідному ізотропному пружному півпросторі, збуреного залежною від часудовільно орієнтованою силою.