Провести программную реализацию расчета приблизительного значения функции y, разложив функцию в ряд Тейлора

#include <iostream>
#include <cmath>
#include <iomanip>
#include <cstdlib>
#include <clocale>
 
using namespace std;
 
double S(double x, double e)
{
    double   k = 1;
    double sum = 0;
    double ak  = x - 1;//рассчитываем лишь самый первый слен разложения
    while(e < fabs(ak))//Вот оно условие окончания итерационного процесса
    {
        sum += ak;
        //вот рекуррентная формула для ln
        ak  *=((-1)*(x - 1)*k/(k + 1));
        k   =  k + 1;
    }
    return sum;
}
 
int main()
{
    setlocale(LC_ALL, "");
    double xstart, xend, xdelta, e;
    cout << "Начало (xstart):";
    cin >> xstart;
    cout << "Конец (xend): ";
    cin >> xend;
    cout << "Шаг(xdelta): ";
    cin >> xdelta;
    cout << "Точность(e, eps): ";
    cin >> e;
    if (xstart <= 0)
    {
        cout << "Error!" << endl;
    return 0;
    }
    for(double x = xstart; x <= xend; x += xdelta)
    {
        if (x <= 0)
            continue;
        else
            cout << "f(" << x << ")="  << setw (9) << log(x) << " == " << S(x, e) << endl;
    }
 
    system("pause");
    return 0;
}