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

#include <iostream>
#include <iomanip>
#include <cmath>
 
int main()
{
    using namespace std;
 
    cout << "Введите x [в радианах]: ";
    double x;
    cin >> x;
 
    cout << "Задайте точность вычислений: ";
    double e;
    cin >> e;
    if (e < 0) e = -e;
 
    unsigned k = 1;
    short minus_one_power = 1;
    long double x_power = x * x * x;
    unsigned long long three_power = 27ull, denominator = 6ull;
    long double numerator = (three_power - 3) * x_power, term = minus_one_power * numerator / denominator, sum = term;
    while (abs(term) > e)
    {
        ++k;
        minus_one_power = -minus_one_power;
        three_power *= 3 * 3;
        x_power *= x * x;
        numerator = (three_power - 3) * x_power;
        denominator *= 2 * k * (2 * k + 1);
        term = minus_one_power * numerator / denominator;
        sum += term;
    }
 
    cout << "Правая часть выражения: " << setprecision(20) << sum / 4 << ", количество слагаемых: " << k << endl;
    cout << "Левая часть выражения: " << pow(sin(x), 3) << endl;
 
    return 0;
}