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

#include <iostream>
#include <cmath>
using namespace std;
 
int main(){
    int k = 0;
    double x, e, sum = 0.;
    double s0, s1; // члены ряда на предыдущей и текущей итерациях
    
    cout << "input x: ";
        cin >> x;
    cout << "input (0 < e < 1): ";
        cin >> e;
 
    // вычислим первый член ряда
    double fact = k + 1; // будем отдельно вычислять часть с факториалом из знаменателя
    s1 = pow(-1, k) * pow(x, k + 2) / fact / (4 * k + 3);
    sum += s1;
    
    // начинаем вычислять сумму ряда с заданной точностью
    do{
        k++;
        s0 = s1;
        fact *= (k + 1);
        s1 = pow(-1, k) * pow(x, k + 2) / fact / (4 * k + 3);
        sum += s1;
    } while(fabs(s1 - s0) > e);
    
    cout << "сумма из " << k << " членов ряда равна " << sum << endl;
    
    return 0;
}