эпсилон

[Actec]

#include <iostream>
#include <cmath>
using namespace std;

double fun(double x)
{
	return (log10(x) / sqrt(x + 1));
}

void TrapezeIntegral()//метод трапеций
{
	double a = 1.0, b = 5.0, h, x, S = 0, fa, fb, f,eps;
	int i, n = 1000;

	x = a;
	h = (b - a) / n;
	fa = fun(a);
	fb = fun(b);
	f = h * (fa + fb) / 2.0;
	for (i = 1; i < n; i++)
	{
		x = a + i * h;
		f += h * fun(x);
	}
	cout << "Trap Integral = " << f << endl;
}

void RectangleIntegral()//вычисление интеграла методом правых прямоугольников
{
	double a = 1.0, b = 5.0, h, x, f = 0, eps;
	int i, n = 1000;

	h = (b - a) / n;
	for (i = 1; i <= n; i++)
	{
		x = a + i * h;
		f += fun(x);
	}
	f *= h;
	cout << "\nRect Integral = " << f << endl;
}
int main()//вызов функции
{
	TrapezeIntegral();
	RectangleIntegral();
	system("pause");
	return 0;
}

Помогите внести в код эпсилон(точность)

[Alex11223]

Сообщение от Actec

Помогите внести в код эпсилон(точность)

Ну, бери за правый край, а я за левый, и понесли!

[digitalis]

Дык это же просто, как веник:

	for (i = 1; i <= n; i++)
	{
		x = a + i * h;
                y = fun(x) ;
		f += y;
                if abs(y) < eps break ;
	}

Ну если это с for, а так while напрашивается.