Формула синуса на Си

[vinternete]

даны 2 числа X i N
найти сумму по формуле
формула синуса :
X − X^3/(3!) + X^5/(5!) − … + (−1)^N·X^(2·N+1)/((2·N+1)!)
значение должно быть ~ sin(x)

#include<stdio.h>
#include<math.h>
main()
{
        int x,i,n,fa;
        double s=1,f=1,re=1,j;
        printf("vvedite N - ");
        scanf("%d",&n);
 
        printf("vvedite X - ");
        scanf("%d",&x);
        for(i=1;i<=n;i++)
        {
                f*=i;
                fa=2*i+1;
                
for(j=1;j<=fa;j++)
re*=j;
 
                s+=(pow(-1,i)*pow(x,fa))/re;
re=1;
        }
        printf("Virajenie ravno= %f\n",s);
printf("sin(x)=%f\n",sin(x));
}

где ошибка?

[KobolD]

А вот зачем тебе f если ты ее нигде не используешь

f*=i;

Попробуй цикл с i от нуля до n делать, а то помоему у тебя первое слагаемое выпадает.

[Mandrivnyk]

#include<stdio.h>
#include<math.h>

float factor(int n);

int main()
{
    float s, x;
    int n;
    printf("vvedite N - ");
    printf("vvedite X - ");
    scanf("%f", &x);
    s = 0.0;
    for (int i = 0; i <= n; i++)
        s += pow(-1, i) * pow(x, 2 * i + 1) / factor(2 * i + 1);
    printf("Virajenie ravno = %f\n", s);
    printf("sin(x) = %f\n", sin(x));

    return 0;
}

float factor(int n)
{
    float f = 1.0;
    for (int i = 2; i <= n; i++)
        f *= i;
    return f;
}

Факториал сделал действительным, чтобы обойти переполнение.

[Volinad]

Начнем с того, что x и n не целые числа, а вещественные. Исправите- получите результат для синуса в радианах. Но, вычисленная по вашей формуле сумма далека от приближенного значения. Проверьте ваш цикл

[vinternete]

Сообщение от Mandrivnyk

#include<stdio.h>
#include<math.h>

float factor(int n);

int main()
{
    float s, x;
    int n;
    printf("vvedite N - ");
    printf("vvedite X - ");
    scanf("%f", &x);
    s = 0.0;
    for (int i = 0; i <= n; i++)
        s += pow(-1, i) * pow(x, 2 * i + 1) / factor(2 * i + 1);
    printf("Virajenie ravno = %f\n", s);
    printf("sin(x) = %f\n", sin(x));

    return 0;
}

float factor(int n)
{
    float f = 1.0;
    for (int i = 2; i <= n; i++)
        f *= i;
    return f;
}

Факториал сделал действительным, чтобы обойти переполнение.

спасибо