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

#include "stdafx.h"
#include <iostream>
#include <math.h>
#include <iomanip>
typedef double(*uf)(double, double, int &);
void tabl(double, double, double, double, uf);
double y(double, double, int &);
double s(double, double, int &);
 
long double fact(int k)
{
    if (k < 0) 
        return 0; 
    if (k == 0)
        return 1; 
    else 
        return k * fact(k - 1);
}
using namespace std;
int main()
{
    cout << setw(8) << "x" << setw(15) << "y(x)" << setw(10) << "k" << endl;
    tabl(-3, 3, 0.5, 0.0001, y);
    cout << endl;
    cout << setw(8) << "x" << setw(15) << "s(x)" << setw(10) << "k" << endl;
    tabl(-3, 3, 0.5, 0.0001, s);
    system("pause");
}
void tabl(double a, double b, double h, double eps, uf fun)
{
 
    int k = 0;
    double sum;
    for (double x = a; x<b + h / 2; x += h)
    {
        sum = fun(x, eps, k);
        cout << setw(8) << x << setw(15) << sum << setw(10) << k << endl;
    }
}
double y(double x, double eps, int &k)
{
    return pow(9, x);
}
double s(double x, double eps, int &k)
{
    double a, c, sum;
    sum = a = c = x;
    k = 1;
    while (fabs(c) > eps)
    {
        c = pow (log(9.), k)*pow(x,k)/fact(k);
        a *= c;
        sum += a;
        k++;
    }
    return sum;
}