Re #18490 - doctest again.

docs/source/fitfunctions/CrystalFieldHeatCapacity.rst

@@ -68,6 +68,19 @@ The `x`-axis is given in Kelvin, and the heat capacity (`y`-axis) is in Joules p
         if row['Value'] != 0:
             print "%7s = % 7.5g" % (row['Name'], row['Value'])
 
+.. testcleanup:: ExampleCrystalFieldHeatCapacity
+
+.. testoutput:: ExampleCrystalFieldHeatCapacity
+   :hide:
+   :options: +ELLIPSIS, +NORMALIZE_WHITESPACE
+
+        B20 = ...
+        B22 = ...
+        B40 = ...
+        B42 = ...
+        B44 = ...
+    Cost function value = ...
+
 Output (the numbers you see on your machine may vary):
 
 .. code::

docs/source/fitfunctions/CrystalFieldMagnetisation.rst

@@ -78,6 +78,19 @@ The `x`-axis is given in Tesla, and the magnetisation (`y`-axis) is in bohr magn
         if row['Value'] != 0:
             print "%7s = % 7.5g" % (row['Name'], row['Value'])
 
+.. testcleanup:: ExampleCrystalFieldMagnetisation
+
+.. testoutput:: ExampleCrystalFieldMagnetisation
+   :hide:
+   :options: +ELLIPSIS, +NORMALIZE_WHITESPACE
+
+        B20 =  0...
+        B22 =  0...
+        B40 = -0...
+        B42 = -0...
+        B44 = -0...
+    Cost function value = ...
+
 Output (the numbers you see on your machine may vary):
 
 .. code::

docs/source/fitfunctions/CrystalFieldMoment.rst

@@ -77,6 +77,19 @@ The `x`-axis is the temperature in Kelvin, and the magnetic moment (`y`-axis) is
         if row['Value'] != 0:
             print "%7s = % 7.5g" % (row['Name'], row['Value'])
 
+.. testcleanup:: ExampleCrystalFieldMoment
+
+.. testoutput:: ExampleCrystalFieldMoment
+   :hide:
+   :options: +ELLIPSIS, +NORMALIZE_WHITESPACE
+
+        B20 =  0...
+        B22 =  0...
+        B40 = -0...
+        B42 = -0...
+        B44 = -0...
+    Cost function value = ...
+
 Output (the numbers you see on your machine may vary):
 
 .. code::

docs/source/fitfunctions/CrystalFieldSusceptibility.rst

@@ -81,6 +81,19 @@ The `x`-axis is given in Kelvin, and the susceptibility (`y`-axis) is in cgs uni
         if row['Value'] != 0:
             print "%7s = % 7.5g" % (row['Name'], row['Value'])
 
+.. testcleanup:: ExampleCrystalFieldSusceptibility
+
+.. testoutput:: ExampleCrystalFieldSusceptibility
+   :hide:
+   :options: +ELLIPSIS, +NORMALIZE_WHITESPACE
+
+        B20 =  0...
+        B22 =  0...
+        B40 = -0...
+        B42 = -0...
+        B44 = -0...
+    Cost function value = ...
+
 Output (the numbers you see on your machine may vary):
 
 .. code::