Форум программистов, компьютерный форум CyberForum.ru

Необходимо найти ошибку при вычислении ряда Тейлора с заданной точностью - C++

Восстановить пароль Регистрация
 
Equality
0 / 0 / 0
Регистрация: 11.12.2012
Сообщений: 4
26.03.2013, 13:39     Необходимо найти ошибку при вычислении ряда Тейлора с заданной точностью #1
Доброе время суток! Есть задача. Посчитать ряд Тейлора с заданной точностью. Сама функция - cos2x. (Разложение прикрепляю во вложенные файлы к теме).

Проблема заключается в том, что при вычислении в итоге получаю неопределенность. Результат вида -1.#IND. Деления на ноль вроде бы нету, проверил. Тем более, изначально вычисления идут нормальные, потом после какого-то Nого повторения цикла вылазит неопределенность. С чем это связано? Переполнения типа double? Или компилятор не знает какой тип данных брать и дабл суёт в интеджер?
В общем, нужна помощь и совет, как это всё исправить

Написанный код:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
#include <locale.h>
#include <iostream>
#include <cmath>
#include <cstdio>
using namespace std;
 
double n;
double result;
double x;
double eps;
 
 
    double factorial(double N) 
  {
      return !N ? 1 : N * factorial(N - 1);
  }
 
 double power(double Z, int N)
 {
      int i=0;
      if (N==0) return 1.0;
      for(i=0; i<N; i++) Z*=Z;
      return Z;
 }
 
double suma(double x, double eps)
 {
    double Abs=10; double prev, sum=1.0;
    for(n = 1; ; n++)
    {
        sum+=((power(-4.0,n))*(power(x,2.0*n)))/(factorial(2.0*n));
        Abs=fabs(sum);
        if((n>50)||(Abs<eps)) break;
        prev=sum;
        printf("->%lf\n",prev);
    }
    return prev;
 }
 
int main()
{
    setlocale( LC_CTYPE, ".1251" );
    cout.precision(6);
    cout<<"Введите точность: ";
    cin>>eps;
    cout<<"\nВведите x: ";
    cin>>x;
    if(x>-1 && x<1)
    {
        printf("\nОтвет: %lf", suma(x, eps));
        cout<<"\ncmath: "<<cos(2*x)<<endl;
    }
    else
    cout<<"\nОшибка: x может быть в пределах (-1, 1)."<<endl;
}
Изображения
  
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
26.03.2013, 13:39     Необходимо найти ошибку при вычислении ряда Тейлора с заданной точностью
Посмотрите здесь:

C++ Найти сумму ряда с заданной точностью
Найти сумму ряда с заданной точностью c C++
C++ Найти сумму ряда с заданной точностью, если дан общий член ряда
Найти сумму ряда с заданной точностью = 10 C++
C++ Найти сумму ряда с заданной точностью
C++ Найти сумму ряда с заданной точностью
C++ Найти сумму ряда с заданной точностью
C++ Функция бесконечного ряда Тейлора с заданной точностью на С

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

Текущее время: 11:22. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru