Вычисление функции через ряд

[Sorento2]

спасибо большое! проверю сейчас))

[Sorento2]

В принципе все понятно. Но работать пока не захотел. Наверно, что то сам не то делаю. Пробую 10-ый или 11-ый раз:D

[Stilet]

Бывает... )

[Sorento2]

#include <math.h>
#include <stdio.h>
#include <conio.h>
double factorialch(double a){
	double value=1.0;
	for(double i=1.0;i<a+1.0;i++) value*=i;
	return value;
}


//double fact(double a)
//{
//if(a == 1)
//return a;
//else return fact(a-1)*a;
//} 


int main()
{
      double n, b, xn, xk, sh, t,  x, sum;
      double i;
      //char string[105];
      char Y,N;
      
      a1: while(1)
      {  
      printf("Vvedite x nachal'noe, x konechnoe, sdvig, shag, tochnost'\n");       
      scanf("%lf%lf%lf%lf%lf", &xn,&xk,&b,&sh,&t);
      {
               
            
printf("+---------------------------------------+\n");
printf("¦    X    ¦cos(x+b) ¦ Function¦    N    ¦\n");
printf("¦---------+---------+---------+---------¦\n");
      for(x=xn;x<xk+sh/8;x+=sh)
      {
    n=0.0;
    sum=0.0;

    while(fabs(pow(x,n)*cos(b+n*M_PI/2)/factorialch(n))>t)
    {
    if(n>10000.0)
    break;  
    n++;
    }
    for(i=0.0; i<n; i++)
    sum+=(pow(x,i)*cos(b+3.14*i/2))/factorialch(1.0*i);
       
printf("¦ %6.2f  ¦ %6.3f  ¦ %7.3f  ¦  %g   ¦\n", x, cos(x+b), sum, n);
      }
printf("+---------------------------------------+\n");


a2: printf("Prodolzhit' vichysleniya? (Y/N)\n");




}
if(getche()=='Y') goto a1;
if(getche()=='N') return 0;
else goto a2;
}
return 0;
}

Вроде как работает, но почему то после 35-ого элемента сбивается ряд, и
не совпадают результаты.

[Stilet]

Ну а отладчиком проходил?
Смотрел что происходит на 36-й итерации?

[Sorento2]

Вот как то так.

[Stilet]

Хм... А ты случайно в своем факториале за пределы переменной не выскакиваешь?
Возможно получаешь при вычислении число, которым процессор просто давится, вот он тебе и пишет -1.#IO, мол "отвали от меня со своими плюшками в километр"

[Sorento2]

Я понимаю, что из-за факториала не считает, но не знаю как исправить это.