Driving anomalous coefficients near 0 at Tc.

Not quite at 0 due to differences in fitting procedure (?).
Key change: solving (1-lambda)^2 = 0 requires a smaller tolerance than
solving (1-lambda) = 0. So solve sqrt(abs(1-lambda)^2) instead.