Нахождение корня n-ой степни

[Lord_Rufus]

как найти корень n-ой степени?
a - число под корнем
n - степень корня

1. Power(a, 1/n)
2. exp(ln(a)/n)

у кого еще какие варианты?

[NeshSoft]

всегда пользовался Power, а зачем нужны еще варианты?

[Lord_Rufus]

neshsoft
У меня такое задание, такое решение не принимается

[ОДИНОЧЕСТВО В СЕТИ]

пользуйся!

function Degree(const DegBasis, DegParam: Extended): Extended; 

asm 

  // Вход: DegBasis, DegParam --> в стеке 

  // Выход: Result --> ST(0) 

  // Примечание: 

  // В случае некорректных данных в ST(0) возвращается DegParam 

 

  XOR     EDX, EDX 

  FLD     DegParam 

  FTST 

  FNSTSW  AX 

  SAHF 

  JNZ     @@DegParam_is_not_0 

  FFREE   ST 

  FLD1 

  JMP     @@exit 

@@DegParam_is_not_0: 

  JC @@DegParam_is_less_than_0 

  JMP     @@cmp_DegBasis_to_0       

@@DegParam_is_less_than_0: 

  OR EDX, 1 

@@cmp_DegBasis_to_0: 

  FLD     DegBasis 

  FTST 

  FNSTSW  AX 

  SAHF 

  JNZ     @@DegBasis_is_not_0 

  TEST    EDX, 1 

  JZ @@DegParam_is_greater_than_0 

  FSTP    ST 

  JMP     @@exit 

@@DegParam_is_greater_than_0: 

  FSTP    ST(1) 

  JMP     @@exit 

@@DegBasis_is_not_0: 

  TEST    BYTE PTR DegBasis + 9, 80H 

  JZ @@DegBasis_is_greater_than_0 

  FABS 

  FXCH 

  OR EDX, 2 

  FLD     ST 

  FRNDINT 

  FSUB    ST, ST(1) 

  FTST 

  FNSTSW  AX 

  SAHF 

  FSTP    ST 

  JZ @@DegParam_is_integer 

  FLD1 

  FDIV    ST, ST(1) 

  FABS 

  FLD1 

  FCOMP     

  FNSTSW  AX 

  SAHF 

  JC @@1_div_Abs_DegParam_greater_or_equal_to_1 

  JZ @@1_div_Abs_DegParam_greater_or_equal_to_1 

  FSTP    ST 

  FSTP    ST(1) 

  JMP     @@exit 

@@1_div_Abs_DegParam_greater_or_equal_to_1: 

  FISTP   QWORD PTR @@Int_64 

  TEST    BYTE PTR @@Int_64, 1 

  JNZ     @@continue 

  FSTP    ST(1) 

  JMP     @@exit 

@@DegParam_is_integer: 

  FLD     ST 

  FISTP   QWORD PTR @@Int_64 

  TEST    BYTE PTR @@Int_64, 1 

  JNZ     @@continue 

  XOR     EDX, EDX 

@@continue: 

  FXCH 

@@DegBasis_is_greater_than_0: 

  FYL2X 

  FLD     ST 

  FRNDINT 

  FSUB    ST(1), ST 

  FXCH    ST(1) 

  F2XM1 

  FLD1 

  FADD 

  FSCALE 

  FSTP    ST(1) 

  TEST    EDX, 2 

  JZ @@exit 

  FCHS 

  JMP     @@exit       

@@Int_64:          

  DQ 0 

@@exit: 

  FWAIT 

 

end;

пример использования

degree( 4,1/2);

вернет 2

[Вавел из ГМТУ]

Автор будет в восторге

[NSvirus]

Я находил по формуле Бинома Ньютона..

[Lord_Rufus]

ОДИНОЧЕСТВО В СЕТИ
Слишком сложный

[ОДИНОЧЕСТВО В СЕТИ]

вот с коментами правда чуть другой

function New_Power(x, y: FLOATTYPE): FLOATTYPE; //abs(y*log2(x))<990

asm

fld y

fld x

fldz             //Сравним основание степени

fcomip st,st(1)  // с 0 и соответственно установим  флаги процессора

je @Zero

FYL2X            //Стек: ST(0)=t=y*log2(x)

fldz

fcomip st,st(1)  //Флаги выставляем соответственно числу 0-y*log2(x)

ja @Reverse      //Если 0>y*log2(x), то сосчитаем 2**|y*log2(x)|, а после инвертируем

call Core_Exp2

jmp @Exit

@Zero:

fxch st(1)

fstp st         //Освобождаем ненужный регистр

jmp @Exit

@Reverse:

fabs             //Берем абсолютную величин

call Core_Exp2

fld1             //Считаем обратное значение:

fdivrp           //1/(2**|y*log2(x)|)

@Exit:

end;

[Lord_Rufus]

NSvirus
это ряд? дык вроде она для квадратных корней?

[ОДИНОЧЕСТВО В СЕТИ]

по формуле Бинома Ньютона..

он для целых степеней!