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

216 MATLAB programs
fout = fopen(‘SH2-FD.ezp’, ‘w’);
fprintf (fout, ‘/et g “%s”\n’, tit);
fprintf (fout, ‘/et x “%s”\n’, titx);
fprintf (fout, ‘/sm off\n’);
fprintf (fout, ‘/sd off\n’);
fprintf (fout, ‘/og on\n’);
fprintf (fout, ‘%15.5f %15.5f %15.5f\n’, [W;real(G);imag(G)]);
fclose (fout);
f = 1/(2*pi*mu);
T1 = [0 1 1];
T2 = [1+dt:dt:tmax];
Uyy = f./sqrt(T2.ˆ2-1)+real(f./sqrt((T2.ˆ2-rˆ2)));
Uyy = [0 0 inf Uyy];
Uyy(3) = Uyy(4); % avoid singularity & discontinuity in plot
T = [T1, T2];
plot (T, Uyy);
titx = ‘Dimensionless time, t*Cs/r’;
xlabel(titx);
title(tit);
grid on;
pause;
% Create EasyPlot file
fout = fopen(‘SH2-TD.ezp’, ‘w’);
fprintf (fout, ‘/et g “%s”\n’, tit);
fprintf (fout, ‘/et x “%s”\n’, titx);
fprintf (fout, ‘/sm off\n’);
fprintf (fout, ‘/sd off\n’);
fprintf (fout, ‘/og on\n’);
fprintf (fout, ‘%15.5f %15.5f\n’, [T;Uyy]);
fclose (fout);
close;