> >>> s.poly(s.roots([1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20])).real > array([ 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11., > 12., 13., 14., 15., 16., 17., 18., 19., 20.]) Sure, but I agree with Heiko that tiny zeros should be detected internally and cleaned up automatically. eric