Деление на ноль

[sithif]

Привет,

если в коде вычислять функцию $f(r)=\frac 1 r \cdote (1-\exp(-r^2))$ в нуле то получится, что мы делим на ноль. Если посчитать предел этой функции в нуле, то получим ноль. Вопрос такой, как лучше запрограммировать такой случай?
У меня два предложения есть. Можно сделать условный цикл,

if(r==0)
   f(r)=0
else
   f(r)=...
end if

Другой вариант, это всегда переопределять $ r=r+\delta$ , а дельта выбрать маленьким (например, 1.E-40, получается ноль).
Как думаете, что лучше? Может еще способы есть?

[rendall]

Я делал как вы, т.е. условие на случай неопределённости.

[Yuri Gendelman]

sithif в сообщении #452488 писал(а):

Может еще способы есть?

Например (r >= 0):

if(r < 1.0e-10) // подобрать нужное значение константы
   f(r)=0
else
   f(r)=...
end if

В арифметике с плавающей точкой ставнение $x==x_0$ не рекомендуется. Надежнее использовать что-то вроде $abs (x - x_0) \le \epsilon$

[venco]

А ещё лучше для маленьких $r$ перейти на ряд Тейлора, т.к. выражение $1-\exp(-r^2)$ будет очень неточным.

[sithif]

Спасибо за советы!