Отправлено продвигающийся 21 октября 2003 г. 04:46 В ответ на: Ответ: отправлено ВН 20 октября 2003 г. 15:31

1. вопрос возник вот отсюда

Y=4/3*ln(is)+pow*ln(2)

is=C(1-x)=2^N(1-x)

Y=4/3[N*ln(2)+ln(1-x)] + pow*ln(2)

так вот в апликухе на техасе вычисляется именно ln(is) - и результат в Q11 - так как is<8207 и (4/3)*ln(8207)<=12

а мне показалось что перекинуть N "поближе" к pow - и тем самым считать только ln(1-x) - и это уже в Q15
вот в этом был вопрос - впроечем видимо апликуха на техасе кривовата - в очередной раз -
к тому же они хоть и приводят в качестве исходной для расчета формулу Y=4/3*ln(is)+pow*ln(2)
но так потом нигде ни словом о ее полном вычислении - что только про логарифм и экспоненту.

да экспоенту считать так предложено

e^(x) = e^(x_int)*e^(x_fract)

при x_fract < 0.5 это все быстро сходится, иначе медленно
та что если x_fract > 0.5
то применим "хитрость" :)

e^(x) = e^(x_int+1)*e^(x_fract-1)

тут вроде все чисто и у них
x_int - малое - в моем случае 1 или 2
в их случае до 12 вроде

2. да интерес вызывает возможность посчитать is^(4/3) в лоб без изысков. Через тейлора - таким образом только один ряд считать в отличие от 2 при ln & exp.

Вопрос только отчкуда этот диапазон для x: 0.5 - 1 ?

дальше в 4 степень - это всего 4 умножения - думаю можно не сильно волноваться :)

вот тут например
http://www.exponenta.ru/educat/class/courses/ma/theme14/theory.asp#2

есть (1+x)^m

но что то не видно про область сходимости - хотя ряд вида [x^(n)]/n! куда сходиться то? - ну думаю что при x<1 он точно сходиться и вполне успешно. Остальное видимо нам не столь интересно.
Спасибо.

• Ответ: — ВН   (21.10.2003 11:43, 506 байт)