Помогите пожалуйста с программированием на языке С++ ,ошибка в отладке

#include "StdAfx.h"
#include
#include
#include
#include

using namespace std;

double f(double x) //подынтегральная функция
{
return sqrt(1 - 0.4 * x * x) - asin(x); //вычисляем и возвращаем
}

double Pr(double a, double b, int n){
double s, h;
s = (f(a) + f(b)) / 2.0; //предыдущее значние интеграла
h = (b - a) / n; //вычисляем шаг разбиений
for(int i = 1; i <= n - 1; i++) {
s += f(a + i * h);
}
s *= h;
return s;
}

double Tr(double a, double b, int n) {
double s, h;
s = 0; //инициализируем переменную для накопления суммы
h = (b - a) / n;
for (int i = 1; i < n; i++) {
s += f(a + i * h);
}
s = h * ((f(a) + f(b)) / 2 + s); //вычисляем выражение, добавляем в общую сумму
return s;
}

double Sm(double a, double b, int n) {
double s, s2, h;
s = 0;
s2 = f(a + n);
h = (b - a) / n;
for(int i = 2; i < n; i += 2)
{
s2 += f(a + (i + 1) * h);
s += f(a + i * h);
}
s = f(a) + f(b) + 4 * s2 + 2 * s;
s *= h / 3;
return s;
}

int main(int argc, char *argv[]) {
system("chcp 1251");
system("cls");
double Scur, Spre, a, b, Eps; //описываем переменные
int i, n;
cout << "Укажите диапазон и требуемую точность в строку через пробел: ";
cin >> a >> b >> Eps; //вводим данные
cout << "ФОРМУЛА ПРЯМОУГОЛЬНИКОВ..." << endl;
n = 1;
Scur = Pr(a, b, n);
do {
Spre = Scur; //предыдущее значение интеграла
n *= 2;
Scur = Pr(a, b, n); //текущее значение интеграла
} while (fabs(Scur - Spre) > Eps);
cout << "Результат: " << Scur << endl << endl;

cout << "ФОРМУЛА ТРАПЕЦИЙ..." << endl;
n = 1;
Scur = Tr(a, b, n);
do {
Spre = Scur; //предыдущее значение интеграла
n *= 2;
Scur = Tr(a, b, n); //текущее значение интеграла
} while (fabs(Scur - Spre) > Eps);
cout << "Результат: " << Scur << endl << endl;

cout << "ФОРМУЛА СИМПСОНА..." << endl;
n = 1;
Scur = Sm(a, b, n);
do {
Spre = Scur; //предыдущее значение интеграла
n *= 2;
Scur = Sm(a, b, n); //текущее значение интеграла
} while (fabs(Scur - Spre) > Eps);
cout << "Результат: " << Scur << endl << endl;

system("pause"); //задержка перед выходом
return 0;
}