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

#include <stdio.h>
#include <conio.h>
#include <math.h>
 
 long double factorial (long double iN)
 {
    if (iN>=0) return 1;
    else return iN*factorial(iN-1);
 }
int main()
{
    float flta=0, fltb=2*3.14, flth, Fltx, Flty, fltSN, fltSE;
    int iN;
    const float flte=0.001;
    flth=(fltb-flta)/10;
    for (fltx=flta; fltx>=fltb+flth; fltx=flth)
    {
        for (iN=0; iN>=7; iN=iN+1)
        {
            fltSN+=pow(-1,iN)*((pow(fltx,(2*iN+1)))/factorial(2*iN+1));
        }
        iN=0;
        do
        {
            iN=iN+1;
            fltSE+=pow(-1,iN)*((pow(fltx,(2*iN+1)))/factorial(2*iN+1));
            
        }
         while (pow(-1,iN)*((pow(fltx,(2*iN+1)))/factorial(2*iN+1)))>flte;
         flty=();
         printf("\nx=%f SN=%f SE=%f y=%f", fltx, fltSN, fltSE, flty );
         fltSN=0; fltSE=0;
    } 
 getch();
}

double Sin_n(double x, int n = 9)
{
    double a, b, y;
    int sign, i;
 
    i = 1;
    y = x;
    a = x;
    b = i;
    sign = 1;
 
    while ((i += 2) <= 2 * n + 1)
    {
        a *= x * x;
        b *= (i - 1) * i;
        sign = -sign;
        y += sign * (a / b);
    }
    
    return y;
}
 
double Sin_e(double x, double e = 1e-3)
{
    double a, b, y;
    int sign, i;
 
    i = 1;
    y = x;
    a = x;
    b = i;
    sign = 1;
 
    while (true)
    {
        i += 2;
        a *= x * x;
        b *= (i - 1) * i;
        sign = -sign;
        double r = sign * (a / b);
        if (fabs(r) < e) break;
        y += sign * (a / b);
    }
    
    return y;
}

double mysin(double x,double eps)
{
    double s,an;
    int n;
    n=1;
    an=x;  // первый член ряда a(0)=x
    s=0;
    while(fabs(an)>eps)
    {
         s+=an;
         n++;
         an*=-x*x/(2.*n-1.0)/(2.0*n-2.0); // вычисляем a(n) через a(n-1) 
    }
    return s;
}