c++ and matlab agree to order 1e-6 on initial p, minusDS and DDS