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

#include <iostream>
#include <conio.h>
#include <math.h>
 
using namespace std;
 
double func1(double x, double eps)
{
    double s = 0;
    double an;
    int n;
    n = 1;
    an = log (1-x);
 
    while ( fabs(an) > eps)
 
    {
        s += an;
        n++;
        an *= - (x * n) / n;
 
        cout << "\n" << s;
    }
    return s;
}
 
int main ()
 
{
 
  double x;
  cin >> x;
  double rezult = func1(x,1e-6);
 
  return 0;
 
}

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
int main()
{    
    int n;
    double xstart, xend, dx, an, S, e, eps=1e-9; 
    cout <<"xstart="; cin>>xstart;
    cout <<"xend="; cin>>xend;
    cout <<"dx="; cin >>dx;
    cout <<"e="; cin >>e;
    
    cout <<"   x      S(x)        Y(x)\n";
    for (double x=xstart; x<=xend+eps; x+=dx)
    { 
        if (-eps<x&&x<eps) x=0.0;
        n=2;
        an=-x;
        S=0.0;
        while(fabs(an)>=e)
        {
            S+=an;            
            an*=x*(n-1)/n;            
            n++;               
        }
    cout <<setprecision(2)<< fixed<< " "<<x<<"  "<<setprecision(7)<<S<<"  "<<log(1.-x)<<"  n="<<n-2<<endl;
    }   
    
system("pause");
return 0;
}