Figures from funda_fig_13.for and knopoff_hudson.for
A. Figures

A.1 Exact Solutions and Stress coefficients RLFRHD(I), TLFRHD(I) from laterally-averaged structure at x = 0 .

A.2 Plots for FLFRHD(I), RATIOD(I), RLFRHL(I), TLFRHL(I), FLFRHL(I), RATIOL(I) from laterally-avearged structurs at x = 0 . The serrated plot shows extremely minute differences in the range of \pm 4 x 10^-16

A.3 Plots for RLFRHF(I), TLFRHF(I), FLFRHF(I), RATIOF(I), RLFRHN(I), TLFRHN(I) from normal Alsop's Figure 13 structure.
B. Additional Figures for Paper 1b

B.1 Figure 7 These results come from program funda_fig_13.for

B.2 Figure 8 . These results come from program knopoff_hudson.for . We have a loss of 3 sig. fig. here. During the root bracketing procedure, these were computed using 7 zeros (i.e. 7 sig. fig.)
ISTART = 2900000
DO 5000 ITERAT=ISTART,7400000
C = DBLE(ITERAT)/1000000.00D+00 + 0.00000001D+00
The computation for 8 zeros has yet to finish (tomorrow). Do you want me to use much higher precision, like increase the zeros from to 9 if that still fails to removes the serated nature?? . or do you want me to just smooth the figure out manually?