Ряд Тейлора

26.10.2013, 14:35. **Просмотров** 383. **Ответов** 0

Здраствуйте, задача состоит в просчете ряда Тейлора в области от -1 до 1, с вводом шага и точности.
Я уже набросал кое-что, вот код:

Код

#include <iostream>
#include <conio.h>
#include <math.h>
#include <iomanip>
using namespace std;
int fct(int arg)
{
return (arg == 1) ? 1 : arg * fct(arg - 1);
}
int main(){
long double f;
double Xstart=-1,Xend=1,x,dx,e,realF;
int n;
cout<<"Vvedite E"<<endl;
cin>>e;
cout<<"Vvedite shag"<<endl;
cin>>dx;
//------
cout<<"--------------------------\n";
cout<<"   x   |      F     | n |  real function) ";
cout<<endl;
//pow(x,n)*fct(n-3)/(fct(3)*fct(n));
for(x=Xstart;x<=Xend; x+=dx)
 {
  f = 0;
  realF = pow((x-1),-4);
  n = 0;
  while(fabs(realF-f)>e)
   {
    n++;
    f -=pow(x,n)*fct(n+3)/(fct(3)*fct(n));
   }
 cout<<setw(6)<<setprecision(2)<<x<<" | "<<setw(10)<<setprecision(7)<<f<<" | ";
 cout<<setw(3)<<n<<" | "<<setw(10)<<setprecision(7)<<realF<<endl;
 
}
	
	getch();
return 0;
}

Но в выводе отображаются не числа, а 1.#IND, подскажите в чем проблема, и правильно ли составлена программа.
Вот сам ряд Тейлора:

Опции темы