I would like to suggest that the vnl_real_polynomial class documentation be altered to mention that in the coefficient array coefficients[i] is the coefficient on X^(degree-i).  I recently had a bug in my code resulting from assuming that coefficients[i] is the coefficient on X^(i).

--Eric