Разложение в ряд Тейлора

[alina999]

Программа не выводит ответ, подправьте пожалуйста код
Вывести на экран таблицу значений функции Y(x) и ее разложения в ряд S(x) с точностью ε (табл. 8.1). Вывести число итераций, необходимое для достижения заданной точности. Вычисление S(x) и Y(x) оформить в виде функций.

Код: (C++)

#include <iostream>
#include <stdlib.h>
#include <string>
#include <string.h>
#include <locale>
#include <math.h>

using namespace std;


double f(double x)
{
return pow(1+x,1/3);
}

int main()
{

setlocale(LC_ALL, "Rus");

double a = -0.9, b = 0.9;
double h = 0.1;
double eps = 0.001;

for (double x = a; x <= b; x += h)
{
double y = f(x);
double sum = pow(1+x,1/3);
double tmp;
int k = 2; while (fabs(sum - y) > eps)
{
k++;
tmp/= pow(3,k)*k*k;
sum += pow(-1,k-1)*(3*k-4)*pow(x,k)*tmp;
}
cout << "x=" << x << " y(x)=" << y << " S(x)=" << sum << " k=" << k << endl;
}
system("pause");
return 0;
}

Добавлено через 3 минуты и 53 секунды:
сумма ряда начинается с 2
s(x)=1 +x/3+ сумма раяда с 2 *(-1)^k-1*(1*2*5*8*...*(3K-4))/((3^k) *k!)
y(x)=3 sqrt(1+x)

[RXL]

А сама?

[Вад]

alina999, что значит, "не выводит"? cout вон там есть. Какие будут идеи, почему до него очередь не доходит?

[alina999]

Точнее не то, что нужно, не могу понять как факториал организовать правильно и вообще как с вашей точки зрения должна выглядеть программа
alina999, что значит, "не выводит"? cout вон там есть. Какие будут идеи, почему до него очередь не доходит?

Добавлено через 1 минуту и 56 секунд:
RXL,  я бы сюда не обратилась, если бы знала сама

[Aether]

Код: (C++)

...

double f(double x)
{
    return pow(1 + x, 1 / 3);
}


int main()
{
    setlocale(LC_ALL, "Rus");

    double a = -0.9, b = 0.9;
    double h = 0.1;
    double eps = 0.001;

    for (double x = a; x <= b; x += h)
    {
        double y = f(x);
        double sum = pow(1 + x, 1 / 3);
        double tmp;
        int k = 2;

        while (fabs(sum - y) > eps)
        {
            k++;
            tmp /= pow(3, k) * k * k;
            sum += pow(-1, k - 1) * (3 * k - 4) * pow(x, k) * tmp;
        }

        cout << "x=" << x << " y(x)=" << y << " S(x)=" << sum << " k=" << k << endl;
    }
...
}

А выводится, кстати, должно. Но суть алгоритма туманна, так как "y" и "sum" считают одно и то же. И цикл "while" вообще не работает.

[Вад]

alina999, факториал можно организовать двумя способами:
1. просто: сделать отдельную функцию, которая только его и вычисляет. В цикле или рекурсивно (полагаю, для учебных нужд глубина рекурсии не будет такой уж большой):

Код: (C++)

double factorial(int k) {
  if (k <= 1) {
    return 1.0;
  }
  else {
    return k * factorial(k-1);
  }
}

2. эффективно - для этого уже нужно понимать, почему плохо вычислять по новой факториал для каждого нового члена, т.е. иметь представление, что такое за штука O(n^2). Но код уже будет не такой простой и понятный.