Разложить функцию в ряд

double snf(int n, double x)
{
    double sn=0,slag=1,st=x;
    for(int i = 0; i < n; i++)
    {
        slag = st * sin(i*(PI/4));
        st *= x;
        sn += slag;
    }
    return sn;
}
///////////
double f(double x)
{
    return (x*sin(PI/4))/(1-(2*x*cos(PI/4)+x*x));
}
 
/////////
cout<<f(10)<<"="<<snf(40,10);

double snf(int n, double x)
{
    double sn=0,slag=1,st=x;
    for(int i = 1; i <= n; i++)
    {
        slag = slag  * x;
        sn += slag * sin(i*(PI/4));
    }
    return sn;
}

double snf(double x,double eps){/*eps - это погрешность, т.е. ряд, если он сходится, 
с каждой итерацией увеличевается на всё меньшее значение и когда он при какой либо
итерации увеличится меньше чем, например 0.001, это значит мы посчитали ряд с погрешностью
0.001, соответственно вам нужно задать погрешность. Если же ряд расходится то при
каждой следующей итерации ряд будет увеличиваться на всё большую величину, в 
таком случае сумма ряда будет равна бесконечности.*/
    int i=1;
    double x1,x2=-1;
    do{
        x1=x2;
        x2=pow(x,i)*sin(i*PI/4);
    }while(abs(x2-x1)>=eps);
    return x2;
}

sn += (slag *= x) * sin(i*(PI/4));