То ли баг, то ли я говнокод написал

#include <math.h>
#include <iostream>
using namespace std;
 
const double e = 0.0001, n = 5, step = 1/n;
int    fact             (int); // Вычисление факториала
double rowsum        (double); // Вычисление суммы ряда в точке х
double rowsum_epsilon(double); // Вычисление суммы ряда в точке х для точности e
double truefunc      (double); // Вычисление компьютером значения функции в точке х 
 
int main()
{
    cout<<n+1<<" steps: "<<endl;
    for (double i = 0; i <= 1; i += step) 
        cout<<"x = "<<i<<": "<<rowsum(i)<<endl;
    cout<<endl<<e<<" accuracy: "<<endl;
    for (double i = 0; i <= 1; i += step) 
        cout<<"x = "<<i<<": "<<rowsum_epsilon(i)<<endl;
    cout<<endl<<"Computer result: "<<endl;
    for (double i = 0; i <= 1; i += step) 
        cout<<"x = "<<i<<": "<<truefunc(i)<<endl;
    cin.get();
}
 
int fact(int number)
{
    return (number == 0) ? 1 : number * fact(number - 1); 
}
 
double rowsum(double x)
{
    double sum = 0;
    for (int i = 0; i <= n; i++) sum += pow(x, 2*i) * (2*i + 1) / fact(i);
    return sum;
}
 
double rowsum_epsilon(double x)
{
    double sum = 0, el = 0;
    int i=0;
    do sum += el = pow(x, 2*i) * (2*i + 1) / fact(i);
    while ((el > e)&&(++i));
    return sum;
}
 
double truefunc(double x)
{
    return (1 + 2*pow(x, 2))*exp(pow(x,2));
}

    do sum += el = pow(x, 2*i) * (2*i + 1) / fact(i);
    while ((el > e)&&(++i));

    do sum += el = pow(x, 2*i) * (2*i + 1) / fact(i++);
    while (el > e);

el = pow(x, 2*i); 
el = (2*i + 1); 
el = fact(i);
el = pow(x, 2*i) * (2*i + 1) / fact(i++);

for (int i = 0; i <= n; i++) 
{
el = pow(x, 2*i)
el = (2*i + 1); 
el = fact(i);
sum += el = pow(x, 2*i) * (2*i + 1) / fact(i);
}

pow(x, 2*1) * (2*1 + 1) / fact(1);
pow(x, 2*2) * (2*2 + 1) / fact(1);
pow(x, 2*2) * (2*1 + 1) / fact(1);
pow(x, 2*1) * (2*2 + 1) / fact(1);

if((el>e)&&(++i))

if((el>e)
  if(++i)

if(++i)
  if(el>e)