Вычисление значения функции, заданной рядом Тейлора (рекурсия)

#include <iostream>
#include <cmath>
#include <stdlib.h>
#include <conio.h>
#include <clocale>
#include <iomanip>
 
using namespace std;
 
double mysin_x(double x, double eps, long &n)
{
    double s, an;
    n = 0;
    an = 1;
    s = 0;
    while (fabs(an)>eps)
    {
        s += an;
        an *= -x*x /(2.*n +2.0)/(2.0*n+3.0);
        n++;
    }
    return s;
}
 
int main()
{
    setlocale(LC_ALL, "Russian");
    double xn, xk, dx, eps=0.0000001;
    long n = 0;
    cout << "Программа для вычисления значения функции\nsin(x)/x с помощью ряда Тейлора\n";
 
    cout<<"Начальное значение xn = ";
    cin>>xn;
    cout<<"Конечное значение xk = ";
    cin>>xk;
    cout<<"Шаг dx = ";
    cin>>dx;
 
    cout << "--------------------------------------------------------------------------\n";
    cout << "|  Текущее значение x  |  Текущее значение суммы  |   Счётчик итерации   |\n";
    cout << "--------------------------------------------------------------------------\n";
 
    while(xn <= xk)
    {
        cout<< "|" << setw(15)<< setprecision(8) << xn << "       |    "
            << setw(15) << setprecision(8) << mysin_x(xn, eps, n) <<"       | "
            << setw(15) << n << "      |" << endl;
        xn = xn + dx;
    }
    cout << "--------------------------------------------------------------------------\n";
    _getch();
    return 0;
}

double mysin_x(double x, double& an,double eps, long n)
{
    if(n==1)
        an=1.;
    else       
        an *= -x*x /(2.*n +2.0)/(2.0*n+3.0);
    if(fabs(an)<eps)
       return an;
    return an+mysin_x(x, an,eps,n+1);
}

double an;
int n=1;
double x=0.1;
double y=mysin_x(x, an, 1e-6, n);