Вычисление значения функции заданной разложением в ряд Тейлора

#include <iostream>
#include <cmath>
using namespace std;
int factorial(int i)
{
  if (i==0) return 1;
  else return i*factorial(i-1);
}
int main()
{
setlocale(LC_ALL, "Russian");
float sum=0;
float sum1=0;
float i=0;
int  a, x, n;
float p,j;
float v; v=pow((x*log(a)),n )/factorial(n);
int k; k=pow(a,x);
float t; t=k-v;
cout << " Введите погрешность = ";
cin >> p;
cout << " Введите x = ";
cin >> x;
cout << " Введите a = ";
cin >> a;
cout << " Введите n = ";
cin >> n;
cout << factorial(n);
    while (i<=n)
    {
        j=pow((x*log(a)), i)/factorial(i);
        i= i +1;
        sum=sum+j;
    }
do
{
    t=t+1;
    sum1=sum1+v;
        }
    while (t<=p);
    cout << "Ответ с погрешность равен: " << sum1 << endl;
    return 0;
}

double foo(double x, double a, double eps, unsigned &steps)
{
  assert(a > 0.0); // аргумент логарифма
  x *= log(a);
  steps = 0;
  double sum = 0.0, m = 1.0;
  
  while (fabs(m) > eps)
  {
    sum += m;
    m *= x / ++steps;
  }
  
  return sum;
}

    int a, x, n;
    float p, j;
    float v; v = pow((x*log(a)), n) / factorial(n);