У меня дошли руки выдрать из ДИСП80 алгоритм вычисления квадратного корня, записать его в понятном для всех виде, и
спросить, откуда ноги растут у коэффициентов, взятых из "Computer Approximations" by Hart et al., NY,. 1968.
Вкратце, берём мантиссу аргумента, считаем −0.204405847M^2+0.890194099M+0.313567057, вставляем полученному числу половинный порядок аргумента, домножаем на sqrt(2), если он был нечётный, и делаем две итерации Ньютона y = (x/y + y)/2
Зачем-то там полином считается с включённым округлением, а итерации - с выключенным, это осталось непонятным.
(Кстати говоря, в Диспаке последовательность команд на пару считываний/записей короче, чем в ОС Дубна).
Оказалось, что квадратичный полином получен не подгонкой по точкам, а с помощью алгоритма имени
Евгения Яковлевича Ремеза. Вот надо же было лезть в американскую книжку за результатами применения отечественного алгоритма!
Leo