Вычисление заданной функции разложением ее в ряд с использованием рекурсии - C++

[Fliiz]

Доброго времени суток,народ..Подскажите в чем ошибка

#include <iostream>
#include <clocale>
#include <cmath>

using namespace std;

double add(double x)
{
    double const eps=1E-15;
    double Rn=1.0;
    double y=0.0;
    for (int i=1 ; fabs(Rn)>eps; i++)
    {
        Rn*=x*(3.0*i-2.0)/(3.0*i);
        y=y-Rn;
    }
    return y;
}

int main()
{
    double const eps=1E-15;
    double x;
    double raz=0;
    cout<<"zadaite |x<1| x= ";
    cin>>x;
    double y=0.0;
    double contr=1.0-(1.0/pow((1.0+x),(1.0/3.0)));
    cout<<"contr= "<<contr;
    y=add(x);
    cout<<"res= "<<y;
    raz=y-contr;
    if(fabs(raz)<=eps)
    {
        cout<< "Rezultat vychislenia funkcii:" <<y<< endl;
    }

[Аватар]

кросспост. в той теме нельзя продолжить?
http://programmersforum.ru/showthread.php?t=282515