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

#include <stdio.h>
#include <iostream>
#include <conio.h>
#include <math.h>
#include <stdlib.h>
#include <locale.h>
using namespace std;   
 
 
 
long double fact(int g)
{
    if(g < 0) 
        return 0; 
    if (g == 0) 
        return 1; 
    else 
        return g * fact (g - 1); 
}
 
int main () 
 
{
    setlocale (LC_ALL, "rus");
    char kl;
    printf ("Задание\nВычислите частичную сумму ряда\n\n  S=x-x^2/2!+x^3/3!-x^4/4!+....   \n\nВычисления прекратите, когда модуль очередного слагаемого станет меньше 0,001.\nПроверьте, насколько сумма приблизилась к значению sin(x).\n\n");
    do
    {
        
        int c,b,x,p,g;
        float S;
         float a;
        
        a=1;
        p=0;
        b=1;
        S=0;
        g=1;
        
   
         
     
    
         printf ("Введите x: ");
        scanf ("%d", &x);
        printf (" S =x-x^2/2!+x^3/3!-x^4/4!+....\n");
       
        while (a>0.001)
        {
          a=(-1)*(-1^b)*(x^b)/c;
         
            if (a>0.001)
            {
                printf ("\n (-1)*(-1^%d)*(%d^%d)/(%d)=%.4f\n",b, x,b, c, a);
                b++;
                p++;
                c=fact(g+p);
                //a=(((-1)*(-1^b)*(x^b))/c);
                S=S+a;
                system ("pause");
            }
        }
        printf ("\n S = %.4f\n", S);
        printf ("Нажмите Enter для продолжения или Esc для выхода из программы.\n");
        do
        {
            kl = getch();
                if (kl == 13)
                    printf ("\n");
                else
                    if (kl == 27)
                        return 0;
        }
        while (kl != 13);
    }
    while (1);
}

double mysin(double x,double eps)
{
    double s,an;
    int n;
    n=1;
    an=x;
    s=0;
    while(fabs(an)>eps)
    {
         s+=an;
         n++;
         an*=-x*x/(2.*n-1.0)/(2.0*n-2.0); 
    }
    return s;
}