Разложение в ряд Тейлора

@v0l0d1ka · Регистрация: 14.12.2010

Есть функция log(x);
Нужно ее разложить в ряд и найти сумму для заданной точность е=0.0001
Считаю х от 0 до 6.0, для того чтобы построить график функции. Но при этом код не работает.
Только при х=2.0 получаю результат.
Вот код моей ф-ии

C++

double myf(double x, double e) {
    
        double  stx, lnt;
        int i, fakt, z;
 
        stx = x - 1; //x в степени
        fakt = 1;    //значение знаменателя
        lnt = 0;      // первый элемент в разложении
        i = 1; //{счётчик
        z = 1; //знак элемента
 
        while (fabs(stx/fakt) >= e)
        {
                    lnt=lnt+z*stx/fakt;
                      i++;
                      stx=stx*(x-1);
                      fakt=i;
                      z=z*(-1);  //{смена знака}
        }
        return lnt;
}

Можно ли как-то оптимизировать программу чтобы получить результат для х=6.0 ?
Заранее спасибо

@zss · 12.09.2016, 08:06

Вычисление функций разложением в ряд Тейлора

C++

double my_ln (double x,double eps)
{
    const double lnln=M_LN2; //ln2
    int k=0;
    while(x>2.0)
    { 
        x/=2.0;
        k++;
    }
    x-=1.;
    double s=0;
    int n=1;
    double an=x;
    while (fabs(an)>eps)
    {
        s+=an;
        ++n;
        an*=-x*(n-1)/n;
    }
    s+=k*lnln;
    return s;
}

Опции темы