Паскаль.Почему условие не срабатывает?

[Max$$]

Добрый день! Нужна ваша помощь. По л/р на Pascal написал прогу.

 x: real;
  y: real;
begin
  x := -10;
  while x <= 2 do
  begin
    if (x>=-8) and (x<0) then y:=cos(x);
    if x<-8 then y:=cos(x)/cos(x) + (power (x,0.1*x)); 
    if x>=0 then y:=(power(x,1/3))/sin(x) + x*x*x;
    if (x>-10) and (x<>-8.8) and (x<-7.9) then writeln ('x=',x,'  функция не определена') else writeln('x=',x,'  y=',y);
    x := x + 0.3;
  end;
end.

Вопрос такой: почему условие (x<>-8.8) не срабатывает? Когда х=-8.8 все равно выводится ФУНКЦИЯ НЕ ОПРЕДЕЛЕНА.
Мне самому не нравится, как выглядит эта строка, но я не нашел подходящее исключение, чтобы вместо вывода NaN на экран выводить текст. А он выводится в результате вычисления (power (x,0.1*x)) при отрицательном дробном Х.



________
Код нужно оформлять по правилам:
тегом [CODE]..[/СODE]
(это кнопочка на панели форматирования с решёточкой #)
Не забывайте об этом!

Модератор.

[Serge_Bliznykov]

вещественные числа НЕЛЬЗЯ сравнивать на совпадение.
(т.е. нельзя для них использовать операции =, <=, >= )
Это особенности машинного представления вещественных (дробных) чисел.

почему так - смотрите, например, в темах на форуме:
http://www.programmersforum.ru/showp...05&postcount=9
http://www.programmersforum.ru/showthread.php?t=307730
http://programmersforum.ru/showthread.php?p=1281382

вот, попробуйте заменить строчку на такую:

if (x>-10) and (abs(x+8.8)>0.0001) and (x<-7.9) 
   then writeln ('x=',x,'  функция не определена') else writeln('x=',x,'  y=',y);

только мне не понятно - зачем сначала вычислять, а потом оперделять, что функция не определена?
может быть, стоит делать наоборот - проверять область допустимых значений и только если функция определена, тогда вычислять её?

[Max$$]

Спасибо. Условие сработало, но теперь при значении переменной -8,8 программа выводит ответ NaN, хотя калькулятор говорит, что значение есть (0,147521...). Не понимаю, в чём дело.

Значение -8,8 рассчитывается вот этой функцией:
if x<-8 then y:=cos(x)/cos(x) + (power (x,x/10));

Калькулятор её считает нормально, паскаль выводит NaN. Кто прав?

[Aleksandr H.]

offtop

cos(x)/cos(x)

[Max$$]

согласен ) . Была формула такая дана.

if x<-8 then y:= 1 + (power (x,0.1*x));

Паскаль выдает:

x=-8.8 y=NaN

[Max$$]

хотя -8.8^-0.88 считается на калькуляторе вполне: 0,14752147583696697713276877352338

[Аватар]

Возведение в действительную степень не определено для отрицательных чисел. А power в данном случае считает его действительным и вычисляет как Exp(-0.88 * Ln(-8.8)). Чему равен логарифм от -8.8? То-то и оно. А если показатель степени считать рациональным, то нет проблем - корень 100 степени из (-8.8)^88 ну и обратить - все ok, что и посчитал калькулятор. Умней чем паскаль на компьютере )) По идее калькулятор должен сдохнуть, если предложить посчитать -8.8^-0.87, разве что с комплексными умеет работать

[Serge_Bliznykov]

Сообщение от Аватар

По идее калькулятор должен сдохнуть, если предложить посчитать -8.8^-0.87

угу. но мелкософт так не думает.
windows 7. стандартный виндовый калькулятор спокойно считает -8.8^-0.87
причём, он выдаёт одинаковый ответ и для -8.8^-0.87 и для 8.8^-0.87

[Вадим Мошев]

Сообщение от Аватар

По идее калькулятор должен сдохнуть, если предложить посчитать -8.8^-0.87, разве что с комплексными умеет работать

Если не ошибаюсь, операция возведения отрицательного числа в дробную степень не определена в принципе.
А логарифм от отрицательных чисел вполне вычисляем - но только на множестве комплексных чисел.
Если что, ОБЩАЯ ФОРМУЛА вычисления логарифма записывается так: Ln(z) = ln(|z|) + i*arg(z). У положительных действительных чисел аргумент равен нулю, поэтому мнимой части у результата нет. По такой формуле можно вычистить логарифм любого, отличного нуля, комплексного числа. Действительные числа также являются комплексными.

Но ни в одном из известных мне ЯП логарифм отрицательного числа нигде не давал комплексного результата - в лучшем случае - это NaN, в худшем - ошибка Invalid floating point operation.
Встроенный ЯП в MathCAD - исключение - он вычисляет логарифм на множестве комплексных чисел.

[Max$$]

Т.е. нет программного решения? Если есть, то как в программе описать, я не понял, извините...