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

232 MATLAB programs
function [] = SH Plate(x, z, z0)
% SH line source at elevation zo above the origin in an elastic plate.
% Plate has unit thickness and material properties
% Origin is at the bottom surface, z is positive up
% Arguments:
% x, z = coordinates of receiver
% z0 = elevation of source
% Default data
mu = 1; % Shear modulus
rho = 1; % mass density
cs = 1; % shear wave velocity
h = 1;
% plate thickness
xi = 0.0001; % material damping
N = 10; % parameter for max. time and max. frequency
if z>h | z0>h
disp(‘Error, source/receiver z must be