в ряд разложить же
если координаты X х_1, ..., х_n, а dX = d_1, ..., d_n
то |X+dX|^2= \sum x_i^2+\sum d_i^2 + 2 \sum x_i d_i
кладем A=|Х|, B= A^{-2}(\sum d_i^2 + 2 \sum x_i d_i)
получаем |X+dX|-|Х|=A*(1+B)^{1/2}
раскладываем в ряд по B, которое величина чрезвычайно
маленькая по условию (порядка 10^{-8}).
2-3 значимых членов уже достаточно
разложение в ряд для (1+B)^{1/2} есть например тут
https://en.wikipedia.org/wiki/Square_rootПривет