Разложение в ряд. Не правильно считает.

[Sorento2]

Программа высчитывает значение функции и значение этой функции через ряд в данной точке.
Вроде все высчитывается. Но при х>6 ну или где то около этого начинают огромные числа выводится.
Вводил большую-маленькую точность. Уменьшал шаг. Все равно получается не то.
Посмотрите пожалуйста в чем ошибка.

// asd.cpp : main project file.

#include "stdafx.h"
#include <stdio.h>
#include <conio.h>
#include <math.h>

#define pi 3.1415

long faktorial(int n) // ------------------------
{
int i;
long p=1; // -------------------------------
for(i=1;i<=n;i++)
p*=i;
return p;
}
int main(){
double x0 = 0, xk = 6, Dx = 1, y_0, y_1,toch=0.1;
int  k, sdvig = 0;
double max, a;
double j, y;
for(j=x0; j<=xk; j+=0.1)
{
	y=cos(j);
	y_1=0.0;
	k=0.0;
     do   
    {
		if(k>10)
			break;
        y_0=pow(-1.0, k)*pow(j, 2*k)/faktorial(2*k);
        y_1+=y_0;
        k++; 
    }while(fabs(y_1-y)>=toch);
	printf("%f  %f  %i\n", y, y_1, k);
}


getch();
	return 0;
}

[Sciv]

Почитайте про факториал: это очень быстро растущая функция, как раз от 7 и выше она обычно и выходит за пределы типа int. В дельфях есть типы longint И int64 (не знаю, как они в си), у них факториал можно считать до 12.

[Sorento2]

Так даже если я представляю факториал с типом long double(3,4 * 10 в 4932) у меня все равно выходит то что изображено на картинке.

[Sciv]

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

[Abstraction]

Что-то мне кажется, что у Вас точность неправильно считается. Вроде текущий элемент должен быть больше эпсилон, а не разность текущего с предыдущим.

И это неправильно. Нужно оценивать хвост суммы сверху и требовать, чтобы эта оценка была менее eps. Что не вполне тривиально.

Так даже если я представляю факториал с типом long double(3,4 * 10 в 4932) у меня все равно выходит то что изображено на картинке.

Код изменённой функции faktorial приведите.
А вообще-то, следующее слагаемое нужно считать по предыдущему, а не с нуля.
Также напоминаю, что pow() не гарантирует верного результата, если первый аргумент отрицателен (из-за того, что второй аргумент приводится к типу double, так что можете внезапно начать считать (-1)^(7.000000001), что не определено).

[Sorento2]

у меня сейчас в другом проблема уже)
да, точность я считаю не так(просто здесь пробовал).


этот ряд просто у меня используется в приложении Winapi. и если я указываю типом этой функции факториала long double (или даже double) приложение просто повисает и все.

в простом консольном приложении все считается прекрасно. но в winapi проблема.

long double faktorial(int n) // ------------------------
{
int i;
long double p=1; // -------------------------------
for(i=1;i<=n;i++)
p*=i;
return p;
}

[Abstraction]

этот ряд просто у меня используется в приложении Winapi. и если я указываю типом этой функции факториала long double (или даже double) приложение просто повисает и все.

Установите точный момент, когда происходит зависание.

[Sorento2]

Сразу при запуске приложения. То есть окно с приложением выводится, но дальше просто висит и все. Никакие действия далее не воспринимаются.

[Abstraction]

Сразу при запуске приложения. То есть окно с приложением выводится, но дальше просто висит и все. Никакие действия далее не воспринимаются.

То есть, до вызова функции faktorial? Тогда она тут ни при чём.
То есть, никакой Ваш код не отрабатывает, вообще? Убедитесь в этом. Закомментируйте весь код.