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

function [ ] = SH Plate(x, z, z0) 233
title(tit);
titx = ‘Time’;
xlabel(titx);
pause;
EasyPlot(‘SH plate TD.ezp’, tit, titx, T, U, ‘r’);
% Frequency domain solution via method of images
ni = 20; % Number of images
fac = 0.25/mu/i;
fmax = N*ffund;
nf = 600;
df = fmax/nf;
F = [df:df:fmax];
W = 2*pi*F;
H = zeros(size(F));
for j=-ni:ni
t1 = sqrt(x2+(dz-j*h2)ˆ2)/cs;
t2 = sqrt(x2+(sz-j*h2)ˆ2)/cs;
H = H + besselh(0,2,t1*W) + besselh(0,2,t2*W);
end
H = fac*H;
plot (F, real(H));
hold on;
plot (F, imag(H),‘r’);
tit = sprintf(...
‘Green function at x=%5.3f, z=%5.3f due to SH line load in plate at z=%5.3f’,...
x, z, z0);
titx = ‘Frequency (rad/s)’;
grid on;
title(tit);
xlabel(titx);
pause;
hold off;
EasyPlot(‘SH plate FD1.ezp’, tit, titx, F, H, ‘c’);
% Frequency domain solution via normal modes
nm = 20;
% number of normal modes
fac = 1/mu/i;
tit = sprintf(...
‘Green function at x=%5.3f, z=%5.3f due to SH line load in plate at z=%5.3f’,...
x, z, z0);
x = abs(x/h);
z = z/h;
z0 = z0/h;
W = W*h*(1-i*xi)/cs;
H = 0.5*exp(-i*x*W)./W;
W2 = W.ˆ2;
for j=1:nm
j1 = pi*j;