Нужна программа на C++ для подсчёта ряда

//============================================================================
// Name        : series_ln.cpp
// Author      : UranFlex @gmail.com
// Version     : 0.1 alpha
// Copyright   : UranFlex 2013
// License     : Free for all
// Description : Рассчитать сумму ряда для ф-ции ln( ( 1 + x ) / ( 1 - x ) ) с заданной точностью
//               C++, Ansi-style
//============================================================================
 
#include <iostream>
#include <cmath>
 
using namespace std;
 
// функция вычисляет для x с точностью eps значение суммы ряда ln( ( 1 + x ) / ( 1 - x ) )
// и записывает вычисленное значение в sum
// если ( | x | >= 1 или eps == 0 ), то функция вернет false,
// если функции переданы правильные данные, то она возвратит true
bool CalcSumLn( const double x, const double eps, double &sum );
 
int main() {
 
    // запрашиваем данные у пользователя
    double x, eps;
    cout << "Программа вычисляет сумму ряда для функции  ln( ( 1 + x ) / ( 1 - x ) )" << endl;
    cout << "Введите значение аргумента, | x | < 1: ";
    cin >> x;
    cout << "Введите точность вычислений, | epsilon | > 0: ";
    cin >> eps;
 
    // пробуем вычислить сумму ряда
    double sum;
    if ( !CalcSumLn( x, eps, sum ) ) { // eсли сумму нельзя посчитать,
        cerr << "Введены неверные данные!" << endl; // то выводим сообщение об ошибке
        return 1; // выходим из программы
    }
 
    // а если все нормально - выводим результаты
    cout << "Для заданного x: " << endl;
    cout << "Сумма ряда с заданной точностью = " << sum << endl;
    cout << "Значение, вычисленное с помощью библиотечной функции " << log( ( 1 + x ) / ( 1 - x ) ) << endl;
 
    return 0;
}
 
bool CalcSumLn( const double x, const double eps, double &sum ) {
 
    if ( fabs( x ) >= 1 || eps == 0 ) // если функция получила неверные данные
        return false; // возвращаем false
 
    double epsAbs = fabs( eps ); // на случай, если точность вычислений передана отрицательной
    double curr = x; // первый элемент ряда
    sum = curr; // сумма пока равна первому элементу
 
    unsigned counter = 3; // следующее число для вычисления очередного элемента ряда ( для второго это 3)
 
    while ( fabs( curr ) >= epsAbs ) { // пока очередной элемент ряда по модулю больше точности вычислений
        curr *= x * x * ( counter - 2 ) / counter; // вычисляем очередной элемент ряда
        sum += curr; // прибавляем его к общей сумме
        counter += 2; // увеличиваем следующее число для вычисления очередного элемента ряда
    }
 
    sum *= 2; // вычисляем конечную сумму ряда
 
    return true; // так как все прошло гладко вернем true
}