Форум программистов, компьютерный форум, киберфорум
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Sugarless
0 / 0 / 0
Регистрация: 20.12.2015
Сообщений: 7
#1

Вычисление приближенного значения tan(x) через ряд Тейлора/ряд Маклорена - C++

20.12.2015, 21:26. Просмотров 439. Ответов 3
Метки нет (Все метки)

Привет, CF.
Задача: Вычисление приближенного значения tan(x) через ряд Тейлора/ряд Маклорена
Проблема: не совпадают результаты вычислений с гугл-калькулятором. Вразумите в ошибку, пожалуйста.
Формулы:
Числа Бернулли
Ряд Тейлора
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
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
#include <iostream>
#include <math.h>
#include <iomanip>
#include <omp.h>
#define  PI          3.14159265358979323846
double tgx(double); //функция, высчитывающая тангенс
long long  fact(int); //факториал
double vetBern(double); // числа Бернулли
 
int main() {
    int n;
    std::cin >> n;
    double start = omp_get_wtime();
    double Result = tgx((n * PI) / 180);
    double end = omp_get_wtime() - start;
    std::cout << std::fixed << std::setprecision(7) << " result = " << Result << std::endl;
    std::cout << "elapsed time = " << end << std::endl;
 
    return 0;
}
 
double tgx(double x)
{
    double z = 0.0;
    int i;
    if (fabs(x) < PI / 2) {
        for (i = 1; i < 14; i++)
        {
            //это ряд Тейлора z += vetBern(2 * i) *  ((pow(2, 2 * i)) * (pow(2, (2 * i)) - 1)  * pow(x, (2 * i) - 1)) / fact(2 * i);
            z += (vetBern(2 * i) * pow(-4, i) * (1 - pow(4, i)) * pow(x, ((2 * i) - 1))) / fact(2 * i);// ряд Маклорена
            std::cout << z << std::endl;
        }
    }
    else std::cout << " Wrong variable" << std::endl;
    return z;
}
 
long long fact(int  x) { // факториал
    if (x == 0 && x == 1)
        return 1;
    else {
 
        long long result = 1;
        for (int i = 2; i <= x; i++) {
            result *= i;
        }
        return result;
    }
}
 
 
double vetBern(double n)// расчет чисел Бернулли
{
    static_cast<int>(n);
    double bern = 1.0;
    double sum;
    if (n == 0)
    {
        sum = 1.0;
    }
    else
    {
        sum = 0.0;
        for (int k = 1; k <= n;k++)
        {
            sum += ((fact(k + 1) / (fact(n + 1)*fact(k - n))) * vetBern(n - k));
        }
        bern = ((-1) / (n + 1)) * sum;
    }
    return fabs(bern);
}
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
20.12.2015, 21:26
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Вычисление приближенного значения tan(x) через ряд Тейлора/ряд Маклорена (C++):

Вычисление функции через разложение в ряд (Ряд Тейлора) - C++
Привет всем. Задание такого плана: Вычислить и вывести на экран в виде таблицы значения функции, заданной с помощью ряда Тейлора,...

Вычисление значения функции разложением в ряд Тейлора - C++
Дано вещественное число X (|X| &lt; 1) и целое число N &gt; 0. Написать функцию для вычисления суммы Σ = x - x2/2 + x3/3 - ... + (-1)N-1xN/N....

Вычисление значения функции разложением в ряд Тейлора - C++
Помогите пожалуйста написать программу на C++. От неё зависит мой зачёт. Нарушение правил форума п. 5.18 Запрещено размещать задания и...

Вычисление значения функции разложением в ряд Тейлора - C++
Составить программу для вычисления суммы ряда с точностью \varepsilon. Под приближенными значениями суммы ряда с точностью \varepsilon,...

Вычисление значения функции разложением в ряд Тейлора - C++
Всем доброго времени суток, ребят, кто может решить этот пример с пояснениями? я не совсем понимаю сути этой задачи, и уже начинает плыть...

Вычисление значения функции разложением в ряд Тейлора - C++
Рассчитать сумму членов бесконечного ряда с заданной пользователем точностью е для заданного пользователем значения x є (-1;1) :...

3
olya7
477 / 294 / 91
Регистрация: 18.02.2013
Сообщений: 684
20.12.2015, 21:43 #2
long long дает целое число -> целое/целое=целое
0
Sugarless
0 / 0 / 0
Регистрация: 20.12.2015
Сообщений: 7
20.12.2015, 21:56  [ТС] #3
Сделал факториал не int, а double - результат еще отдалился от истины.
Как правильнее было бы типы привести?
0
olya7
477 / 294 / 91
Регистрация: 18.02.2013
Сообщений: 684
20.12.2015, 22:21 #4
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
double vetBern(double n)// расчет чисел Бернулли
{
    static_cast<int>(n);
   
    double *ber=new int n;
     double bern = 1.0;
     ber[0]=1;
    double sum;
    if (n == 0)
    {
        sum = 1.0;
    }
    else
    {
        sum = 0.0;
        for (int k = 1; k < n;k++)
        {
            sum += ((-1) / (n + 1)) * ((fact(k + 1) / (fact(n + 1)*fact(k - n))) * ber[n - k]);
            ber[k]=sum;
            
        }
       // bern = sum;
    }
    delete []ber;
    return sum;
}
1
20.12.2015, 22:21
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
20.12.2015, 22:21
Привет! Вот еще темы с ответами:

Вычисление значения функции разложением в ряд Тейлора - C++
Написать программу, вычисляющую значение функции e^x, используя разложение в ряд Тейлора. Параметр x и точность задаются пользователем....

Вычисление значения функции разложением в ряд Тейлора - C++
Доброго времени суток,помогите пожалуйста написать программу.Ломал голову все не как не получается Вычислить с точностью eps&gt;0: ...

Вычисление значения функции разложением в ряд Тейлора - C++
Добрый вечер всем, работаю с c++ около двух недель, дали задание написать программу с Рядом Тейлора по условию, к сожалению, не могу...

Вычисление значения функции разложением в ряд Тейлора - C++
s=x-x^3/3!+x^5/5!-x^7/7!+⋯+x^25/25! Решите задачу)


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

Или воспользуйтесь поиском по форуму:
4
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru