Отправлено TimMatTU 21 октября 2006 г. 11:58 В ответ на: Понятно отправлено <font color=gray>TimMatTU</font> 21 октября 2006 г. 11:45

Офигеть! MAPLE подсказал, что, оказывается
> f:=1/sqrt(x);
> df:=diff(1/sqrt(x),x);
-1/(2*x^(3/2))
> x1:=x-(f/df);
3x
x1:=subs(x=x1,x1);
9x
x1:=subs(x=x1,x1);
81x
При этом получаем последовательные присваивания:
x[1]=3*x[0]
x[2]=9*x[1]
x[3]=81*x[2]
Причём 9=3^2, 81=9^2....
Для невысокого количества итераций они бы могли эти коэффициенты в маленькую табличку запихать и использовать 2 регистра для адресации с пост-инкрементом... :) Теперь всё ясно, боль спаб!!

• страничка #5 ==> — SM   (21.10.2006 13:10 213.141.159.26, пустое, ссылка)
  • Метод ясен: затравка из таблицы 1/sqrt(x), например, из 256 или 2^N значений (индекс легко получить сдвигом) + метод Ньютона Рафсона заключающийся в реккурентности: X[0]=~1/sqrt(x),D=arg/2, X[i+1]=X[i]*(3/2-D*X[i]), Ans=x*X[i+1] — TimMatTU   (21.10.2006 13:44 213.228.95.64, пустое)