Pasted

/**
 * Программа находит решение уравнения ln((x - 1) / 3) / 3 = 0
 * методом Ньютона с подвижным полюсом
 */
 
/******** MAIN.CPP ********/
#include <iostream> // cin, cout
#include <iomanip>  // setprecision
 
using std::cout;
using std::endl;
 
/**
 * Процедура реализует метод 'Ньютона с подвижным полюсом'
 * Возвращает значение X в уравнении ln((x - 1) / 3) / 3 = 0
 * С точностью |X - nextX| < accuracy
 */
extern "C" void __fastcall calculate();
 
extern "C" float X     = 1.5f; // Начальное значение
extern "C" float NextX = NULL; // Следующее значение функции
 
extern "C" float accurecy = 1e-6f; // Точность |X - NextX|
 
int main() {
    cout << X << endl;
 
    calculate();
    cout << std::setprecision(10) << NextX << endl;
 
    return 0;
}
 
 
 
/******** FUNCT.ASM ********/
.386
.MODEL FLAT;
 
EXTERN _X: DWORD
EXTERN _NextX: DWORD
EXTERN _accurecy: DWORD
 
.DATA
BUFFER DWORD 0.0
THREE  DWORD 3.0 ; Константа 3
TWO    DWORD 2.0 ; Константа 2
.CODE
 
; Получаем значение функции  f(x) = ln((x - 1) / 3) / 3
; Постфиксная запись         f(x) = x 1 - 3 / ln 3 /
; Входное значение    _X
; Выходное значение   BUFFER
getF PROC
fld1 ; Для логарифма...
 
FLD _X ; Загружаем X
FLD1   ; Загружаем 1
FSUB   ; Отнимаем от X, 1
 
FLD THREE ; Загружаем в стек 3
FDIV      ; Делим первый раз на три
 
; Натуральный логарифм X = ln(X) = ln(2) * log2(X) = fyl2x(fldln2, X)
; URL: http://www.cyberforum.ru/assembler/thread525843.html
FYL2X  ; st(0) = st(1) * log2( st(0) )
FLDLN2 ; st(0) = ln(2)   st(1) = (X - 1) / 3
FMUL
 
FLD THREE ; Загружаем в стек 3
FDIV      ; Делим второй раз на три
 
FSTP BUFFER ; Сохраняем значение в переменную BUF
 
RET
getF ENDP
 
; Получаем значение производной функции   fdx(x) = 1 / (3 * (x - 1))
; Постфиксная запись                      fdx(x) = 1 3 x 1 - * /
; Входное значение    _X
; Выходное значение   BUFFER
getFdx PROC
FLD1
FLD THREE
FLD _X
FLD1
FSUB
FMUL
FDIV
 
FSTP BUFFER
 
RET
getFdx ENDP
 
; Получаем значение      xk+1 = xk + f(xk) * (f(xk) - 2) / 2 * fdx(xk)
; Постфиксная запись     xk+1 = xk f(xk) f(xk) 2 - * 2 / fdx(xk) * +
; Входное значение  (xk)     _X
; Выходное значение (xk+1)   _NextX
next PROC
FLD _X ; Загружаем X
 
call getF  ; Вычисляем f(xk)
FLD BUFFER ; Загружаем f(xk)
FLD BUFFER
 
FLD TWO
FSUB
FMUL
FLD TWO
FDIV
 
call getFdx  ; Вычисляем fdx(xk)
FLD BUFFER   ; Загружаем fdx(xk)
 
FMUL
FADD
 
FSTP _NextX
 
RET
next ENDP
 
@calculate@0 PROC
; Инициализируем сопроцессор
FINIT
call next;
 
@start:
FLD _X                      ;
FLD _NextX                  ;    Сравниваем точность нашего решения, которая вычисляется
FSUB                        ;    по формуле |_X - _NextX|
FABS                        ;    с заданной точностью, которая хранится в accurecy
                            ;
FLD _accurecy               ;
FCOMIP ST(0), ST(1)         ;
                            ;    Если наше реiение точно, то переходим на метку @exit
JNC @exit                   ;    Иначе продолжаем цикл
 
 
MOV EAX, _NextX
call next         ; Получаем xk+1
MOV _X, EAX
call next ; Какого хуя это тут нужно?!
loop @start
 
 
@exit:
RET
@calculate@0 ENDP
 
 
END