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

[wormix]

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

#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]

*протирает глаза*
http://programmersforum.ru/showthread.php?t=223802
WTF?!