Разложение функции в ряд Тейлора (с++)

[eleonora_zv]

Провести програмну реалізацію розрахунку приблизного значення функції y, розклавши функцію f(x) в ряд Тейлора, при цьому необхідно вивести та застосувати рекурентне відношення для обрахування загальної частини для кожного доданку ряду або його частки. Аргумент функції має змінюватися в діапазоні [Xbegin, Xend] з кроком Xdelta. Параметри Xbegin, Xend, Xdelta та точність (допустима похибка) розрахунків е, яка визначає умову завершення рекурсивного розрахунку, задаються користувачем з клавіатури. Передбачити ситуації, коли значення функції можу бути не визначене для заданих аргументів.Снимок.JPG

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

using namespace std;

double factorial (double p)
{
int n=1;
p=2*n+1;
if (p == 0) 
return 1; 
else 
return p*factorial(p - 1);
}
long double sin_kub(double p, double x, double e, long double y)
{
double fact=factorial(p);
y=((pow(3, p)-3)*pow(x, p)*pow(-1, p-1)/fact);
if (y>e)
return y=(y+sin_kub(p+1, x, e, y))*0.25;
else return y*0.25;
}
void main()
{
setlocale (0,"Russian");
double xbegin, xdelta, xend, e, x;
cout<<"Введите xbegin (начальное значение): ";
cin>>xbegin;
cout<<"Введите xend(конечное значение х): ";
cin>>xend;
while (xend<xbegin)
{
cout<<"Не корректно задано."<<endl;
cout<<"Введите xbegin (начальное значение): ";
cin>>xbegin;
cout<<"Введите xend(конечное значение х): ";
cin>>xend;
}
cout<<"Введите xdelta(шаг): ";
cin>>xdelta;
cout<<"Введите e(точность): ";
cin>>e;
while (e<0.01)
{
cout<<"Ваше \"е\" слишком маленькое. Введите другое: ";
cin>>e;
}
for (double x=xbegin; x<=xend; x+=xdelta)
{
long double f=sin_kub(3, x, e,0);
cout<<"x="<<x<<"      "<<"f(x)="<<f<<endl;
}

getch ();
}

не выводит значения х и у. где ошибка?

[Abstraction]

Провести програмну реалізацію розрахунку приблизного значення функції y, розклавши функцію f(x) в ряд Тейлора, при цьому необхідно (...) не выводит значения х и у. где ошибка?

Well, you could use debugger to find the problem. Actually, it would be much faster then getting answer from some forum on the Internet.
However (as a hint) take note that when you tell computer that

int n=1;
p=2*n+1;

- it does just that. Might as well simply write

int n=1;
p=3;

And, unless three equals zero, you get yourself an 'infinite' recursion, which actually should deplete the stack quickly and drop your application with some sort of message.

[eleonora_zv]

мне не выдает ошибок и я не пойму в чем проблема.
у меня есть рекулентная формула

y=((pow(3, p)-3)*pow(x, p)*pow(-1, p-1)/fact)

где р зависит от n и все время возрастает

[koljsch]

Оффтоп
Abstraction, жжешь)))

[Abstraction]

мне не выдает ошибок и я не пойму в чем проблема.

Повторяю по-русски: в том, что Ваша функция factorial описывается как "зашёл и не вышел". В отладчике Вы бы это увидели сразу - приучитесь им пользоваться.