Вычисление интеграла методом левых треугольников с заданной точностью

#include "stdafx.h"
#include <iostream>
#include <cmath>
 
double func(const double &x)
{
    return sin(x);
}
 
double integr_lefttriangle(const int &n, const double &a, const double &b)
{
    double sum = 0;
    double step = (b - a) / n;//шаг
    for (double i = a; i<b; i = i + step)
        sum += func(i);//суммируем значения функции в узловых точках
    return sum*step;//множим на величину шага и возвращаем в вызывающую функцию
}
 
using namespace std;
 
 
int _tmain(int argc, _TCHAR* argv[])
{
    setlocale (LC_ALL, "Russian");
 
    int n = 1;
    double a, b, eps, s, s1;
 
    cout << "Введите верхнюю границу интегрирования a = ";
    cin >> a;
    cout << "\nВведите нижнюю границу интегрирования b = ";
    cin >> b;
    cout << "\nВведите требуемую точность eps = ";
    cin >> eps;
    
    s1 = integr_lefttriangle(n,a,b);  //первое приближение
    do {
        s = s1;     //второе приближение
        n = 2 * n;  //увеличение числа шагов в два раза, 
        //т.е. уменьшение значения шага в два раза
        s1 = integr_lefttriangle (a, b, n);
    } while (fabs(s1 - s) > eps);  //сравнение приближений с заданной точностью
    cout << "\nИнтеграл = " << s1 << endl;
 
    system("pause");
 
    return 0;
}