Eduardo Kausel-Fundamental solutions in elastodynamics_ a compendium-Cambridge University Press (2006)

220 MATLAB programs
a = sqrt(a2);
dt = (tmax-a)/np;
% Cs/Cp
% time step
r = sqrt(xˆ2+zˆ2);
% Source-receiver distance
c = abs(z)/r;
% direction cosine w.r.t. z axis
s = abs(x)/r; % ” ” ” x axis
theta = asin(s);
% Source-receiver angle w.r.t. vertical
crang = asin(a);
% Critical angle w.r.t. vertical
fac = cs/(pi*mu*r); % scaling factor for displacements
% t < tp=a (two points suffice)
T0 = [0 a];
Uxx0 = [0, 0];
Uxz0 = [0, 0];
Uzx0 = [0, 0];
Uzz0 = [0, 0];
T = [a+dt:dt:tmax]; % Time vector
T2 = T.ˆ2;
jl = length(T);
if theta > mang
% Displacements on surface
T1 = [a:dt:1];
% interval from tp to ts
T2 = [1+dt:dt:tmax]; % t > ts=1
% a=tp