Версия для печати темы

Нажмите сюда для просмотра этой темы в обычном формате

Форум «Всё о Паскале» _ Free Pascal, Pascal ABC и другие _ Вычисления. Модуль Math

Автор: klem4 19.10.2006 12:38

Охх ничего с утра опнять не могу. Почему вылетает программа при X = 0 ?

{$mode objfpc}

uses Math, Crt;

var

  x: Double = 0;
  c: array [0..4] of Double = (-1, -1, 0, 2, 1);

function F: Double;
var
  i: Integer;
begin
  result := 0;
  for i := Length( c ) - 1 downto 0 do
   result := result + Power(x, i) * c[i];
end;

begin
  clrscr;
  writeln(F:2:2);
  readln;
end.

Автор: klem4 19.10.2006 12:49

В общем ошибка при Power(0, ...)

Чтож придется делать проверку перед возведением.

Автор: volvo 19.10.2006 13:04

Цитата(rtl.pdf)
17.12.53 power
Synopsis: Return real power.
Declaration: function power(base: float;exponent: float) : float
Visibility: default
Description: power raises base to the power power. This is equivalent to exp(power*ln(base)).


Чему равен Ln(0), можно полюбопытствовать?

Автор: klem4 19.10.2006 13:16

Все ясно, спасибо.

Автор: klem4 19.10.2006 17:00

Как же тогда допускается такое ?!

uses crt, math;
begin

  clrscr;

  writeln(Power(-2, 2):2:1); // выводится 4
  writeln(Power(-2, 3):2:1); // выводится -8

  readln;
end.


Выходит Ln(-2) ?!!!!

:mega_chok: :blink:

Автор: volvo 19.10.2006 18:13

Посмотри на реализацию Power (файл math.pp):

function power(base,exponent : float) : float;

  begin
    if Exponent=0.0 then
      if base <> 0.0 then
        result:=1.0
      else
        InvalidArgument
    else if (base=0.0) and (exponent>0.0) then
      result:=0.0
    else if (abs(exponent)<=maxint) and (frac(exponent)=0.0) then
      result:=intpower(base,trunc(exponent))
    else if base>0.0 then
      result:=exp(exponent * ln (base))
    else
      InvalidArgument;
  end;


(InvalidArgument как раз и бросает Exception)... Видишь, что произойдет, если дробная часть экспоненты ненулевая? Исключение выброшено не будет... А попробуй
  writeln(Power(-2, 2.5) : 2 : 1);
smile.gif