Вычисляю экспоненту по формуле Тейлора:
long double cust_exp(double x) {
long double rememberX = x;
long double returnValue = 1. + x;
/* for (int k = 2; k <= 150; k++) {
x = x * rememberX;
returnValue = returnValue + x/fact(k);
}*/
returnValue = 1;
for (int i = 150; i > 0; --i) {
returnValue = 1. + x * returnValue / i;
}
return returnValue;
}
long double fact(int k) {
long double n = 1;
long double returnValue = 0;
for (int i = 1; i <= k; i++) {
n = n * i;
}
returnValue = n;
return returnValue;
}
Закоменченное это то же самое через факториал, только из-за факториала дольше. Нужно с точностью до 6ти знаков после запятой по сравнению с оригинальной функцией. Но вот при степенях больше 20 начинаются расхождения уже на 6м знаке. Есть ли какой-то другой способ расчёта, более точный?
fact
.