Вычислить сумму ряда с заданной точностью

#include <iostream>
#include <cmath>
using namespace std;
 
double sumN(int n,double x)
{
    double s=0;
    for(int i=0;i<n;i++)
        s+=(4*(pow(x,(i+1)*4+1))/((i+1)*4+1));
    return s;
}
double sumEps(double eps, double x)
{
    double s1=0,s2=0;
    for(int i=0;s2-s1<eps;i++)
        {
        s1+=(4*(pow(x,(i+1)*4+1))/((i+1)*4+1));
        s2=s1+(4*(pow(x,(i+2)*4+1))/((i+2)*4+1));
        }
    return s2;
}
int main()
{
    double s=0,x=2;
    int t,n;
    for(;;)
    {
    cout<<"To enter n press 1.\nTo enter epsilon press 2.\nTo exit press 0.\n";
    cin>>t;
    switch(t)
    {
    case 1:
        cout<<"Enter n."<<endl;
        cin>>n;
        s=sumN(n,x);
        cout<<s<<endl;
        break;
    case 2:
        double eps;
        cout<<"Enter epsilon."<<endl;
        cin>>eps;
        s=sumEps(eps,x);
        cout<<s<<endl;
        break;
    case 0:
        return 0;
        break;
    }
    }
    return 0;
}