Правильно ли написан код для метода трапеций на с++?

#include "stdafx.h"
#include
#include
#include

using namespace std;

double f(double x)
{
double E = 0.000000001, S = 1, R = 0, N, n = 0, y, Y = 1, F = 1;

//Разложение 1/sqrt(1-x) в ряд Маклорена и извлечение корня
for (int i = 1; i < 10; i++){
S += pow(x,i);

}

return sqrt(S)*(2 * x / sqrt ((1+x)*(1+x*x)));
}

int _tmain(int argc, _TCHAR* argv[])
{
setlocale(LC_ALL, "Russian");
cout « "Программа считает интеграл 2x/sqrt(1-x^4) с границами от 0 до 1 методом трапеции";
cout « "с точность E=10^(-9)!" « endl;
double a = 0, b = 1, n = 4, S = 0, h, integ, IntegPrev;
h = (b - a) / n;
for (int i = 1; i <= (n - 1); i++)
S += f(a + h*i);
integ = h * ((f(a) + f(b)) / 2 + S);
const double PI = 3.14159265;
//посчитали первый интеграл для n=4
do {
//IntegPrev = integ;
n += 1;
h = (b - a) / n;
S = 0;
for (int i = 1; i <= (n - 1); i++)
S += f(a + h*i);
integ = h * ((f(a) + f(b)) / 2 + S);
} while (abs(PI/2 - integ) > 0.01);
cout « "\nAnswer: " « integ « endl;
getchar();
return 0;
}