вычислить функцию путем разложения ее ряда через рекурсию

не получается вычислить функцию с помощью ряда с заданной точностью eps

double f(double a, double b, double &p, int l)
{
    if (l == 0)
        return 1;
    else
        p *= pow(a, 2 * l + 1) * ((1 / fabs(l)) + 1 / (2 * l + 1));
    if (fabs(p) < b)
        return p;
    return p + f(a, b, p, l + 1);

}   
int main()
{
    setlocale(LC_ALL, "Rus");

    double x = 2;
    double eps = 0.0001;
    double s = 1;
     cout << f(x, eps,s,1) << "\n";
     cout << 2 * x - x * pow(eps, pow(-x, 2)) - sin(x);

}