diff --git a/Packages/Number_types/changes.txt b/Packages/Number_types/changes.txt
index 27e9512eff4..ca895e0838e 100644
--- a/Packages/Number_types/changes.txt
+++ b/Packages/Number_types/changes.txt
@@ -1,3 +1,6 @@
+4.25     (25 June 01)
+- operator<<(MP_Float) was missing the approximation when zero.
+
 4.24     (25 June 01)
 - to_interval(Gmpz) is exact when < 2^53.
 
diff --git a/Packages/Number_types/src/MP_Float.C b/Packages/Number_types/src/MP_Float.C
index 888bbcc40da..4f0e9648e87 100644
--- a/Packages/Number_types/src/MP_Float.C
+++ b/Packages/Number_types/src/MP_Float.C
@@ -269,7 +269,7 @@ operator<< (std::ostream & os, const MP_Float &b)
 {
   // Binary format would be nice and not hard to have too (useful ?).
   if (b.is_zero())
-    return os << 0;
+    return os << 0 << " [ double approx == " << 0.0 << " ]";
 
   MP_Float::const_iterator i;
   int exp = b.min_exp() * log_limb;