Программа по вычислению Ряда Тейлора

double a = (x-1);
double S  = 0;
for(n=2; fabs(a) > E; n++) {
  S += a;
  a = - a*(x-1)/n
}

#include <sstream>  
#include <iostream> 
#include <iomanip>  
#include <cmath>    
 
using namespace std;
 
int main()
{
    double Dx;
    double xBeg;
    double xEnd;
    double E,n;
 
    cout << "Enter x begin > ";   cin >> xBeg;
    cout << "Enter x end > ";     cin >> xEnd;
    cout << "Enter Dx > ";        cin >> Dx;
    cout << "Enter precision > "; cin >> E;
 
    string part;
    cout << setw(6) << "x" << setw(10) << "S" << setw(20) << "n" << endl;
    cout << "--------------------------------------------------" << endl;
 
    if (!cin || xBeg <= 0 || xEnd > 2) {
        cerr << "Wrong input\n";
        return 1;
    }
 
    errno = 0;
    for (double x = xBeg; x <= xEnd; x += Dx) {
 
        double a = (x - 1);
        double S = 0;
        for (n = 2; fabs(a) > E; n++) 
        {
            S += a;
            a = -a*(x - 1) / n;
        }
 
        if (errno > 0) {
            cerr << "Fail calculate\n";
            return 1;
        }
 
        cout << setw(6) << x << setw(10) << S << setw(20) << n << endl;
    }
    system("pause");
    return 0;
}

for(n=2; fabs(a/n) > E; n++) {