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

#include <iostream>
#include <iomanip>
using namespace std;
 
double power(double x, int k) {
double result = 1;
for (int i = 1; i <= k; i++)
result *= x;
return result;
}
 
int fact2(int k) { // двойной факториал
int result1 = 1, m;
if (k % 2 == 0) {
m = k / 2;
while (m > 0) {
result1 *= k;
k = k - 2;
m = m - 1;
}
}
else {
m = (k + 1) / 2;
while (m > 0) {
result1 *= k;
k = k - 2;
m = m - 1;
}
}
return result1;
}
 
double f1(double x, double eps)
{
int i = 0; double sum = 1, slag = 1;
while (abs(slag) > eps) {
i++;
slag = power(-1, i) * power(x, i) / fact2(i + 2);
sum += slag;
}
return sum;
}
 
double f2(double x, double eps) {
int i = 0; double sum = 1, sum2 = 1, slag = 1, slag1 = 1;
while (abs(slag) > eps && abs(slag1) > eps) {
i++;
if (i % 2 == 0) {
slag *= x / i;
sum += slag;
}
else {
if (i == 1) { slag1 = -x / 3; }
slag1 *= x / i;
sum2 += slag1;
}
}
return sum, sum2;
}
 
int main() {
setlocale(LC_ALL, "Russian");
double A, B, dx, Eps, N;
cout << "Введите отрезок [A,B] в форме A B" << '\n';
cin >> A >> B;
if (A == B) {
cout << "Вы ввели точку,а не отрезок" << '\n';
return 0;
}
cout << "Введете Eps от 0.001 до 0.00001 (Eps-точность)" << '\n';
cin >> Eps;
if (Eps > 0.001 || Eps < 0.00001) {
cout << "Неверно введён Eps" << '\n';
return 0;
}
cout << "Введите N (кол-во точек табулирования) N>0" << '\n';
cin >> N;
if (N < 1) {
cout << "Вы ввели N<0" << '\n';
return 0;
}
dx = (B - A) / (N - 1);
double m = N;
double z = A - dx;
cout « "F1(x)" << "\t" << "F2(x)" << '\n';
while (m > 0) {
z += dx;
cout << f1(z, Eps) << "\t" << f2(z, Eps) << '\n';
m--;
}
system("pause");
return 0;
}