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

function [T, Uxx, Uxz, Uzz] = lamb2D(x, z, pois) 223
plot (T,Uxz);
tit = sprintf(...
‘Horizontal displacement due to vertical line load, \\nu=%5.2f’,pois);
title (tit);
grid on;
%axis([0 tmax -0.5 0.5]);
if theta > mang
hold on
plot ([tr tr], [0 wr]);
end
pause;
hold off;
pfac = rho*cs*r;
EasyPlot(‘Lamb2Dxx.ezp’, tit, titx, T, pfac*Uxx, ‘r’);
EasyPlot(‘Lamb2Dzz.ezp’, tit, titx, T, pfac*Uzz, ‘r’);
if z∼=0
% reverse sign & coupling term for source at depth, receiver at surface
EasyPlot(‘Lamb2Dxz.ezp’, tit, titx, T, -pfac*Uzx, ‘r’);
EasyPlot(‘Lamb2Dzx.ezp’, tit, titx, T, -pfac*Uxz, ‘r’);
else
% source+ receiver are at the surface of a lower half-space
EasyPlot(‘Lamb2Dxz.ezp’, tit, titx, T, pfac*Uxz, ‘r’);
EasyPlot(‘Lamb2Dzx.ezp’, tit, titx, T, pfac*Uzx, ‘r’);
end
close;