Вычисление функции

[XSGRAFF]

Тип проекта - Win32 Project
Задание: Вычислить и вывести на экран в виде таблицы значения функции,заданной с помощью ряда Тейлора,на интервале от хнач до хкон с шагом dx с точностью e.Таблицу снабдить заголовком и шапкой. Каждая строка таблицы должна содержать значение аргумента,значение функции и количество просуммированных членов ряда.

Мой вариант программы:

#include <iostream.h>
#include <iomanip.h>
#include <math.h>
int main()
{
double s,sp,x,xn,xk,xd,e;
int n;
cout<<"vvedi xn,xk,xd,e"<<endl;
cin >> xn>>xk>>xd>>e;
if (xn>=xk) cout<<"ochibka"<<endl;
else
{
cout<<" _______________________________________"<<endl;
cout<<"| n | n |   x    |   F  |        Arth(x)   |"<<endl;
cout<<"|---------------------------------------|"<<endl;
//x=xn;
int m=1; 
for (double x=xn; x<=xk; x+=xd)
{n=0;
sp=0;
	do 
	{
		s=1./(pow(x,(2*n+1))*(2*n+1));
		sp+=s;
	 	n++;
		
	}
	while (fabs(s)>e);
		cout<<"|"<<setw(3)<<m<<"|"<<setw(3)<<n<<setw(9)<<x<<"|"<<setw(12)<<sp<<"|"<<setw(12)<<tanh(x)<<endl;
		m++;
}
cout<<" ---------------------------------------"<<endl;
}
return 0;
}

При прогонке должны совпадать значения 3го и 4го столбца,но они не совпадают. Не могу найти ошибку. Прошу помочь найти ошибку и исправить ее. Заранее спасибо.

[Sazary]

Исходя из вышеприведенной формулы, вроде, должно выглядеть как-то так:

int m=1;
double tmp,a;
for (double x=xn; x<=xk; x+=xd)
{
n=1;
sp=0;
m = 1;
tmp = x;
s = 1/tmp;
while (fabs(s)>e)
	{
		sp+=s;
		m += 2;
		tmp *= (x*x);
	 n++;
	 s = 1/(m*tmp);
	}

-----------
А что это за формула? Гиперболический арктангенс? Так это не tanh, а atanh. Но только она работает с числами от -1 до 1.

[XSGRAFF]

Да - это гиперболический арктангенс

[Sazary]

Вот.

#include <iostream.h>
#include <iomanip.h>
#include <math.h>
#include <conio.h>

int main()
{
double s,sp,x,xn,xk,xd,e;
int n;
cout<<"vvedi xn,xk,xd,e"<<endl;
cin >> xn>>xk>>xd>>e;
if (xn>=xk) cout<<"ochibka"<<endl;
else
{
cout<<" _______________________________________"<<endl;
cout<<"| n | n |   x    |   F  |        Arth(x)   |"<<endl;
cout<<"|---------------------------------------|"<<endl;
//x=xn;
int m=1;
double tmp,a;
for (double x=xn; x<=xk; x+=xd)
{
n=1;
sp=0;
m = 1;
tmp = x;
s = 1/tmp;
while (fabs(s)>e)
	{
		sp+=s;
		m += 2;
		tmp *= (x*x);
	 n++;
	 s = 1/(m*tmp);
	}


	cout<<"|"<<setw(3)<<m<<"|"<<setw(3)<<n<<setw(9)<<x<<"|"<<setw(12)<<sp<<"|"<<setw(12)<<0.5*log((x+1)/(x-1))<<endl;
	m++;
}
cout<<" ---------------------------------------"<<endl;
}
getch();
return 0;
}

Функция atanh считает только при x от -1 до 1. Поэтому нужно использовать другую функцию для вычисления гиперболического арктангенса. Например, логарифм.

[XSGRAFF]

спасибо за помощь