Вычислить сумму ряда (приближенное значение функции sin)

# include <stdio.h>
 # include <conio.h>
 # include <math.h>
#pragma argsused
int main(int argc, char* argv[])
{  int j,f,i,x; float eps,sl,s=0.0;
printf ("vedite x, eps");
scanf("%d %f",&x, &eps);
sl=x;
while(fabs(sl)>eps)
{
s=s+sl;
f=1;
for (j=1; j<=2*i+1;j++) f=f*j;
sl=pow(-1,i)*pow(x,2*i+1)/f;
}
 
 
 
        return 0;

/*
 ============================================================================
 Name        : c_sin_taylor.c
 Author      : UranFlex @gmail.com
 Version     : 0.1 alpha
 License     : Free for all
 Copyright   : UranFlex 2013
 Description : Вычисление sin( x ) с помощью разложения в степенной ряд Тейлора
 *             C, Ansi-style
 ============================================================================
 */
 
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
 
// функция возвращает факториал положительного целого числа
unsigned long long Factorial( unsigned long long value );
// функция избавляется от лишних оборотов единичной окружности возвращает только нужное количество радиан
// от угла, переданного как аргумент в радианах
double ClearUnnTurn( double valuePi );
// функция возвращает сумму ряда sin( x ) с заданной точностью eps
// функция принимает аргумент x в радианах
double SumSeriesSinX( double x, const double eps );
 
int main( void ) {
    // запрашиваем исходные данные у пользователя
    double eps, xDegrees;
    printf( "Введите аргумент x в градусах " );
    scanf( "%lf", &xDegrees );
    printf( "Введите точность epsilon " );
    scanf( "%lf", &eps );
 
    double xRadians = xDegrees * M_PI / 180.0; // переводим градусы в радианы
 
    // выводим результаты вычислений
    printf( "Значение sin( %g )\n", xDegrees );
    printf( " - вычисленное с заданной точностью при помощи разложения в степенной ряд = %2.8f\n", SumSeriesSinX( xRadians, eps ) );
    printf( " - вычисленное при помощи библиотечной функции = %2.8f\n", sin( xRadians ) );
 
    return EXIT_SUCCESS;
}
 
unsigned long long Factorial( unsigned long long value ) {
    return !value ? 1 : value * Factorial( value - 1 ); // с помощью рекурсивного вызова ищем факториал числа
}
 
double ClearUnnTurn( double valuePi ) {
 
    int sign = 1; // знак радианной меры
    if ( valuePi < 0 ) // если радианная мера отрицательна
        sign = -1; // запоминаем это
 
    valuePi = fabs( valuePi ); // берем абсолютное значение радианной меры
    int countUnnecessary = valuePi / ( 2.0 * M_PI ); // ищем количество оборотов
 
    if ( countUnnecessary >= 1 ) // если есть лишние обороты
        valuePi -= ( countUnnecessary * 2 * M_PI ); // значит убираем ненужные
 
    return valuePi * sign; // возвращаем результат, учитывая знак
}
 
double SumSeriesSinX( double x, const double eps ) {
 
    x = ClearUnnTurn( x ); // чистим от лишних оборотов
    double sum = x; // сумма изначально равна x, так как при i = 0 будет именно x
    unsigned counter = 1; // счетчик
    int sign = 1; // для знака очередного элемента ряда
    while ( 1 ) {
        sign *= -1; // для знакочередования меняем знак каждую итерацию
        // эо вычисление очередного элемента ряда по формуле
        double currElement = sign * pow( x, ( double ) ( 2 * counter + 1 ) ) / Factorial( 2 * counter + 1 );
        sum += currElement; // к общей сумме прибавляем очередной элемент ряда
 
        if ( fabs( currElement ) <= eps ) // если очередной элемент ряда меньше или равно модулю эпсилон,
            break; // то необходимая точность достигнута - прерываем цикл
 
        ++counter; // а если больше, то продолжаем цикл, увеличив счетчик
    }
 
    return sum; // возвращаем значение функции, вычисленное с помощью разложения в степенной ряд
}

double MySin(double x,double eps)
{
    double an,s;
    int n;
    n=0;
    an=x;
    s=an;
    while(fabs(an)>eps)
    {
        an*=-x*x/(2*n+2)/(2*n+3);
        s+=an;
        n++;
    }
    return s;
}

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
 
double MySin( double x, double eps );
// функция избавляется от лишних оборотов единичной окружности возвращает только нужное количество радиан
// от угла, переданного как аргумент в радианах
double ClearUnnTurn( double valuePi );
 
int main( void ) {
    // запрашиваем исходные данные у пользователя
    double eps, xDegrees;
    printf( "Введите аргумент x в градусах " );
    scanf( "%lf", &xDegrees );
    printf( "Введите точность epsilon " );
    scanf( "%lf", &eps );
 
    double xRadians = xDegrees * M_PI / 180.0; // переводим градусы в радианы
 
    // выводим результаты вычислений
    printf( "Значение sin( %g )\n", xDegrees );
    printf( " - вычисленное с заданной точностью при помощи разложения в степенной ряд = %2.8f\n", MySin( xRadians, eps ) );
    printf( " - вычисленное при помощи библиотечной функции = %2.8f\n", sin( xRadians ) );
 
    return EXIT_SUCCESS;
}
 
double ClearUnnTurn( double valuePi ) {
 
    int sign = 1; // знак радианной меры
    if ( valuePi < 0 ) // если радианная мера отрицательна
        sign = -1; // запоминаем это
 
    valuePi = fabs( valuePi ); // берем абсолютное значение радианной меры
    int countUnnecessary = valuePi / ( 2.0 * M_PI ); // ищем количество оборотов
 
    if ( countUnnecessary >= 1 ) // если есть лишние обороты
        valuePi -= ( countUnnecessary * 2 * M_PI ); // значит убираем ненужные
 
    return valuePi * sign; // возвращаем результат, учитывая знак
}
 
double MySin( double x, double eps ) {
    x = ClearUnnTurn( x ); // тут
    double an, s;
    int n;
    n = 2;
    an = x;
    s = an;
    while ( fabs( an ) > eps ) {
        an *= -x * x / ( n * ( n + 1 ) ); // тут
        s += an;
        n += 2;
    }
    return s;
}