Си ++ сравнение предыдущего интеграла с текущим и если разница больше выбранной точности - увеличить число разбиений

[lana_ar]

Помогите подправить программу, пожалуйста

#include <iostream>
#include <iomanip>
#include <math.h>

double f(double x)
{
	return pow(x, 2) / pow(pow(x, 3) + 1, 1. / 2);
}

using namespace std;
int main()
{
	cout.fixed;
	cout.precision(6);
	setlocale(LC_ALL, "RUS");
	double a, b, n, h, tr, x[11];
	int i;
	cout << "Введите a"; cin >> a;
	cout << "Введите b"; cin >> b;
	cout << "Введите число разбиений n"; cin >> n;
	h = (b - a) / n;
	x[0] = a;
	i = 0;

	cout << "n" << setw(20) << "x" << endl;
	cout << i << setw(20) << x[0] << setw(20) << f(x[i]) << endl;

	for (i = 0;i < n;i++)
	{
		x[i + 1] = x[i] + h;
		cout << i + 1 << setw(20) << x[i + 1] << endl;
	}
	tr = ((f(x[0]) + f(x[10])) / 2 + f(x[1]) + f(x[2]) + f(x[3]) + f(x[4]) + f(x[5]) + f(x[6]) + f(x[7]) + f(x[8]) + f(x[9]))*h;
	cout << "tr=" << tr << endl;

	for (i = 0;i < n;i++)
	{
		cout << "f(x[i])=" << f(x[i]) << endl;
	}
	return 0;

}

[digitalis]

	double a, b, n, h, tr, tr0, eps;
...........................................
        h = (b-a)/n;
        tr0 = 0 ;
lab:   tr = f(a) + f(b) ;
         x = a;
         for (i = 1;i<n;i++)
    {
          x = a + h ;
          tr = tr + f(x)*2 ;
      }
           tr = tr * h/2 ;
           if (abs(tr-tr0)>eps) 
      { tr0 = tr ; n := n + n ; h = h/2;  goto lab }

Да, применять goto некошерно, но это я оставляю для думанья ТСу на предмет заменить на while. Можно и for с выходом по условию достижения точности, это как кто любит, а мне думать лень, бо в сон клонит. Хотя боюсь, что после кошмара исходника засну не скоро
И коммент неправильный: "Интеграл ревёт" или "Интеграл рванёт" ?

[digitalis]

Поскольку ТСу это уже 100 лет не надо (привычка задать вопрос и уйти на фиг уже стала нормой), напишем это с while

	double a, b, n, h, tr, tr0, eps,d;
...........................................
        h = (b-a)/n;
        d = 1.0 ;
    while (d>eps)
  {  
        tr = f(a) + f(b) ;
         x = a;
         for (i = 1;i<n;i++)
     {
          x = a + h ;
          tr = tr + f(x)*2 ;
      }
           tr = tr * h/2 ;
           d = abs(tr-tr0) ;
           tr0 = tr ;
           n := n + n ;
           h = h/2;
   }

[lana_ar]

Да нет, не пропала! спасибо большое за помощь!

[Alar]

На сколько сложно вопрос написать текстом? зачем форуму нужны ваши картинки? незачем!

[digitalis]

Это был не вопрос, а рецензия.

[Alar]

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

Это был не вопрос, а рецензия.

А как мне брать текст на переименование темы? пришлось вручную перепечатывать, а поисковик так ничего кроме моего заголовка не возмёт никогда.

а в теме с вопросом с картинки вообще нет смысла именно для форума для индексации поисковиками, для поиска по форуму.

А если форумчанину нужна только индивидуальная помощь, то можно и денег за неё заплатить.

[digitalis]

За такой простой вопрос "гусары с да'мав денег не берут" А такая, как теперь, формулировка вопроса делает сам вопрос ненужным: ответ содержится уже в заголовке (дык я и имею в виду - в исправленном).

[Alar]

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

ответ содержится уже в заголовке.

В каком заголовке? название темы я сделал. было СИ++ и всё