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

Ряд тейлора, трудности - C++

Восстановить пароль Регистрация
 
banda
1 / 1 / 0
Регистрация: 06.02.2011
Сообщений: 46
09.11.2012, 00:10     Ряд тейлора, трудности #1
Прошу помочь с разложением cos(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
double x,E;
const long N=12;
 
 
 
 
double mpow(double x, long y)
{
    double r = 1;
    while((y--)>0) r*=x;
    return r;
}
 
long fact(long x)
{
    if (x<=1) return 1;
    else return x*fact(x-1);
}
 
 
 
double tcos(double x)
{
    double r = 0;
    for(long n=0;n<N;n++) {
        r+=mpow(-1,n)*mpow(x,2*n)/fact(2*n);
    }
    return r;
}
 
int main() {
 
    cin>>x;
    cout<<tcos(x);
 
 
    system("pause");
 
    return 0;
}
В коде что-то не так, когда я беру x=1 или x=2, то выводит правильно, а когда допустим 3 и больше то выдает чушь

Очень прошу, помочь, уже долго пытаюсь разобраться и решить задачу, но безуспешно

Добавлено через 4 часа 44 минуты
помогите
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
09.11.2012, 00:10     Ряд тейлора, трудности
Посмотрите здесь:

Ряд Тейлора C++
ряд Тейлора C++
C++ Ряд тейлора
Ряд Тейлора C++
C++ Ряд Тейлора
C++ Ряд Тейлора
C++ Ряд тейлора
Ряд Тейлора С++ C++

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
asidorchenko
379 / 205 / 25
Регистрация: 09.04.2012
Сообщений: 635
10.11.2012, 15:21     Ряд тейлора, трудности #2
Ваша программа вроде бы правильно работает.

Найти cos через разложение в ряд Тэйлора
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
#include <stdio.h>
#include <math.h>
 
#define M_PI 3.14159265358979323846
 
float fcos(int n, float xn, float fact)
{
if (n == 5)
 return 0.0;
 
int fct = 2*n;
 
float nfact;
 
if (n == 0)
nfact = 1.0;
else if (n == 1)
nfact = fact * (float)fct;
else
nfact = fact * (float)(fct - 1) * (float)fct;
 
float coef = pow( (float) -1.0, (float) n) / nfact;
float step = 2.0*n;
 
float result = coef * pow(xn, step);
 
return (result+fcos(n+1, xn, nfact));
 
}
 
int main()
{
float i;
for(i=0.0; i<1.0;i+=0.01)
{
printf("%f: %f ", i, fcos(0, i,1.0));
printf("%f\n", cos(i));
}
 
return 0;
}
Yandex
Объявления
10.11.2012, 15:21     Ряд тейлора, трудности
Ответ Создать тему
Опции темы

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