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

// ихнее67.cpp: определяет точку входа для консольного приложения.
//
 
#include <stdafx.h>
#include <iostream>
#include <math.h>
#include <float.h>
#include <stdio.h>
using namespace std;
 
int _tmain(int argc, _TCHAR* argv[])
{
    setlocale (0, "");
    const int MaxIter=500;
    int n;
    double x,eps=DBL_EPSILON;
    cout<<("Введите любые значения x, где х изменяется -1<=x<1:")<<endl;
    cout<<("Введите аргумент x=");
    cin>>x;
    long double  ln=1, y=ln;
    for (n=1; fabs(ln)>eps && n<MaxIter; n++)
        ln=0;
        ln*=-(pow((double)x,(double)n)/n;
       y+=(-ln);
       cout<<"ln="<<ln<<endl;
    if(n<MaxIter)
    {
        cout<<"\nЗначение функции y="<<y<<"для x="<<x<<endl;
        cout<<"вычислено после"<<n<<"итераций"<<endl;
    }
    else
        cout<< "\nРяд расходится!"; 
}               
    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;
}

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
int main()
{
    int n;
    double xstart, xend, dx, an, S,  eps;
    do {
        cout << "xstart="; cin >> xstart;
    } while (xstart < -1. || xstart >= 1.);
    do {
        cout << "xend="; cin >> xend;
    } while (xend < -1. || xend >= 1.);
    cout << "dx="; cin >> dx;
    cout << "eps="; cin >> eps;
 
    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) >= eps)
        {
            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;
}