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

  • Автор темы maryfly009
  • Дата начала
Статус
Закрыто для дальнейших ответов.
M

maryfly009

#1
Здравствуйте уважаемые программисты! Помогите пожалуйста найти ошибку в программе.
Код:
//Подключение библиотеки стандартного ввода и вывода
#include<stdio.h>
//Подключение библиотеки математических функций
#include<math.h>

//Функция вычисления факториала
double Factorial(unsigned n)
{
//Объявление и инициализация переменных функции
double p=1.0;
unsigned i = 2;
//Цикл вычисления значения факториала числа
for(;i<=n;i++) p*=(double)i;
//Возврат полученного значения
return p;
}

int main()
{
//Объявление и инициализация
//переменных программы
double sum=0, i=0, t=1, x, eps;

//Приглашение к вводу x
printf("x = ");
//Ввод значения x
scanf("%lf",&x);

//Приглашение к вводу eps
printf("eps = ");
//Ввод значения eps
scanf("%lf",&eps);

//Цикл вычисления суммы ряда с заданной точностью
do
{
t *= pow((-1), i) * ((pow(x, 2*i+1)) / (Factorial(2 * i + 1)));
sum += t;
i++;
}while( fabs(t) >= eps );

//Вывод значения суммы 
printf("sum = %5.4lf\n", sum);
//Вывод значения библиотечной функции sin 
printf("sin = %5.4lf\n", sin(x));

return 0;
}
при x=1
sum = 0.8319
sin = 0.8415
вроде значения приближенно равны

а при x = 5
sum = -64289.9875
sin = -0.9589
существенная разница

подскажите в чём проблема
 

Вложения

hosm

* so what *
18.05.2009
2 442
6
#2
мне лично это не понравилось: t *=
запросто могу ошибаться)
 
M

maryfly009

#3
спасибо, действительно в этом была ошибка
 
Z

zeus

#4
А я бы сделал бы так
C++:
#include <stdio.h>
#include <math.h>

double calc_sin(double x, double eps)
{
double res = x, t = 0.0;
int i = 1;

while (abs(x-t) >= eps)
{
t = res;
res *= x*x / (2*i) / (2*i + 1);
} // while

return res;
} // calc_sin

int main(void)
{
double x, eps;

printf("x = ");
scanf("%lf",&x); // Ввод значения x

printf("eps = ");
scanf("%lf",&eps); // Ввод значения eps

printf("calc_sin = %5.4lf\n", calc_sin(x, eps)); // Вывод значение синуса
printf("sin = %5.4lf\n", sin(x)); //Вывод значения библиотечной функции sin 

return 0;
} // main
 
Статус
Закрыто для дальнейших ответов.