Нахождение суммы ряда Тейлора

#include "stdafx.h"
#include <iostream>
#include <conio.h>
#include <math.h>
using namespace std;
 
int main()
{
    double a = 0.1; // задаем параметры границ для a,b
    double b = 0.8;
    double step = (b - a) /10; // определяем шаг
    int n = 10; // вводим n
    double x = a; // даем x значение первого параметра
    double y, c, c1, sum, sum1;
    do
    {
        sum = sum1 = 0; // даем накопителям сумм для разных случаев нулевое значение
                        // случай с n:
        for (int i = 0; i<n; i++)
        {
            c = pow(-1, (i+1)+1)*pow(x,2*(i+1))/(2*(i+1)*(2*(i+1)-1)); // i-й элемент ряда
            sum += c; // суммирование ряда - к предыдущей сумме прибавляется данный элемент
        }
        // случай с e: 
        int i1 = 0; // счетчик элементов 
        do
        {
            c1 = pow(-1, (i1 + 1) + 1)*pow(x, 2 * (i1 + 1)) / (2 * (i1 + 1)*(2 * (i1 + 1) - 1)); // i-й элемент ряда
            sum1 += c1; // суммирование ряда - к предыдущему значению суммы прибавляется данный элемент
            i1++;
        } 
        while (c1 >= pow(10, -4));
        // случай с y:
        y = x*atan(x)-log(sqrt(1+x*x));
        cout<< "X="<< x << ", S(n)="<< sum<< ", S(e)="<< sum1<< ", Y=" << y<< endl;
        x += step; // изменение x на шаг
    } 
    while (x <= b);
    _getch();
    return 0;
}