4-NOV-2016: Checking old CAFE implementation
of new Rousselier integration algorithm,
new3.f.
The solver is
Slatec
dnls1.f routine.
The damage zone with 10800 FEs at 5.7e-3 s:
8-NOV-2016: Negative beta (damage variable).
This is indicative of errors in the solver:
