--- boost/math/tools/polynomial.hpp	2016-08-17 16:58:10.000000000 -0400
+++ boost/math/tools/polynomial_new.hpp	2016-10-15 21:25:15.904928377 -0400
@@ -334,7 +334,7 @@
    }
    T evaluate(T z)const
    {
-      return boost::math::tools::evaluate_polynomial(&m_data[0], z, m_data.size());;
+      return m_data.size() > 0 ? boost::math::tools::evaluate_polynomial(&m_data[0], z, m_data.size()) : 0;
    }
    std::vector<T> chebyshev()const
    {