mypow（）トラフ継続分数とtaylorseriesはpow（）より大きい数値を返します

Question

double taylor_log(double x, unsigned int n){ 
double tmp; 
double sum = 0; 
if(x < 1){ 
    int j = 2; 
    x = 1 - x; 
    tmp = x; 
    sum = -x; 
    for(unsigned int i = 1; i < n; i++){ 
     sum -= ((tmp *= x)/j); 
     j++; 
    } 
    return sum; 
} 
else if (x >= 1){ 
    tmp = ((x-1)/x); 
    for(unsigned int i = 1; i <= n; i++){ 
     sum += (tmp/i); 
     tmp *= ((x-1)/x); 
    } 
    return sum; 

これは正しく動作するテイラー級のログです。 この式を使用して、数値の指数関数を得る。 formula for mypow() これは今私の問題は、私は私の機能のために例えば30回の反復を置けば数）（POWよりも高いことであるPOW

double taylor_pow(double x, double y, unsigned int n){ 
double sum = 1.0; 
int fac = 1; 
double exp = y; 
double lna = taylor_log(x, n); 
for(unsigned int i = 1; i <= n; i++){ 
    fac *= i; 
    sum += (exp * lna/fac); 
    exp *= y; 
    lna *= taylor_log(x, n); 
} 
return sum; 

} ための私のコードです。たとえばpow（2,3）= 8、20回の反復での結果は8.0007 ...そしてその成長です。すべての応答のThans。

Answer 1

int facがオーバーフローしています。それをlong doubleに変更しました。

32ビット符号付きintは最大12の値しか保持しませんが、80または128ビットlong doubleは2000などの値を保持できます。