Возведение в степень на Pascal

[arcman]

накидал свой костыль на эту тему

function pow(x: real; n: real): real;
var
n_int: integer;
begin
     if (ABS(2.0 - n) < 1E-9) then
     begin
          pow:= x * x;
          exit;
     end;

     if (ABS(3.0 - n) < 1E-9) then
     begin
          pow:= x * x * x;
          exit;
     end;

     if (x > 0) then
     begin
         pow:= EXP(LN(x) * n);
         exit;
     end;

     n_int := round(n);
     if (ABS(n_int - n) > 1E-9) then
     begin
          pow:= 0;
          writeln('#ERR# pow()');
          exit;
     end;

     (*if ((n_int MOD 2) = 0) then*)
     if ((n_int AND 1) = 0) then
     begin
          pow:= EXP(LN(ABS(x)) * n);
     end
     else
     begin
          pow:= (-1) * EXP(LN(ABS(x)) * n);
     end;
end;

[ALEKSEEV-1992]

ой-йё зачем всё так усложнять?
Power(x,y);
x:=Число
y:=Степень

[puporev]

ой-йё зачем всё так усложнять?

Если такой умный, набери это в Турбо Паскаль и посмотри что будет, в следующий раз выеживаться не будешь.

[arcman]

угу, в турбо паскаль отлаживал, а писал вот для этого:
http://ru.wikipedia.org/wiki/Structured_Text

так что в контроллерах где не поддерживается Си, всё достаточно плачевно.
с этим ST намучался уже.

[Fellics{новичок}]

Сорри не туда написал

[tmogoreanu]

Вот модуль pow.tpu pow.rar там содержится функция power(a,b).
Т.о
Uses pow;
...
begin
...
y:=power(2,3);
...
end.

[tmogoreanu]

_____________________________

[rainbow]

exp(ln(x)*n)
где x - основание;
n - степень;
x,n должно быть типа real;

[arcman]

Смею напомнить, что требуется реализация ф-и возведения в произвольную степень.

Она должна быть:
1) Правильно работающая
2) Написана на Паскале (максимально совместима/переносима)

[arcman]

я пока оставил так:

function pow(x: real; n: real): real;
var
n_int: integer;
begin
     if (x > 0) then
     begin
         pow:= EXP(LN(x) * n);
         exit;
     end;

     if (ABS(0.0 - x) < 1E-9) then
     begin
          pow:= 0.0;
          exit;
     end;

     n_int := round(n);
     if (ABS(n_int - n) > 1E-9) then
     begin
          pow:= 0;
          writeln('#ERR# pow()');
          exit;
     end;

     (*if ((n_int MOD 2) = 0) then*)
     if ((n_int AND 1) = 0) then
     begin
          pow:= EXP(LN(ABS(x)) * n);
     end
     else
     begin
          pow:= (-1) * EXP(LN(ABS(x)) * n);
     end;
end;

заместо 1E-9 лучше найти способ подставлять точное значение машинного эпсилон, но в моем случае хватает этого.

(простой пример высчитывающий примерное значение эпсилон, в среде TP 7.1, выдает 4.5474735089E-13 в режиме эмуляции сопроцессора и 5.42101086242752E-0020 в противном случае)

Заместо этого

          pow:= 0;
          writeln('#ERR# pow()');

неплохо было бы вернуть NaN, но я пока не знаю как сделать это на Паскаль (мне и не требуется)

в начало можно добавить пару оптимизаций, что бы зазря не гонять тяжелые ф-и ln() и exp()

     if (ABS(0.0 - n) < 1E-9) then
     begin
          pow:= 1.0;
          exit;
     end;

     if (ABS(1.0 - n) < 1E-9) then
     begin
          pow:= x;
          exit;
     end;

     if (ABS(2.0 - n) < 1E-9) then
     begin
          pow:= x * x;
          exit;
     end;

     if (ABS(3.0 - n) < 1E-9) then
     begin
          pow:= x * x * x;
          exit;
     end;