Почему программа зацикливается (легкая)?

Question

[tr1ck]

Почему программа зацикливается? Вроде все верно. Считается интеграл методом трапеции с точностью E=10^(-9).
#include "stdafx.h"
#include
#include
#include

using namespace std;

double f(double x)
{
return 2 * x / sqrt (1-x*x*x*x);
}

int _tmain(int argc, _TCHAR* argv[])
{
double a = 0, b = 0.9, 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);
//посчитали первый интеграл для n=4
do {
IntegPrev = integ;
n += 1;
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);
} while (abs(IntegPrev - integ) > 0,000000001);
cout <<"qq" << integ;
getchar();
return 0;
}

Comments

[tr1ck]

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

[tr1ck]

не пишет сюда библиотеки, но это и не важно.

[Миша Кринкин]

@tr1ck все он пишет, вы просто код вставляете неправильно!!!

Answer 1

[Миша Кринкин]

for (int i = 1; i <= (n - 1); i++)
S += f(a + h*i);

вы забыли обнулить S.

UPD. Вставляйте код как код, а не как текст - не разводите бардак!!!!!!

Comments

[tr1ck]

все равно циклится..

[Миша Кринкин]

@tr1ck ну покажите исправленный код, ошибка с обнулением S в любом случае есть.

[tr1ck]

А все, вдвоем помогли) спасибо) запятая и S)

Answer 2

[Алексей Павлов]

while (abs(IntegPrev - integ) > 0,000000001);
Запятая вместо точки справа от знака "больше". Нужно написать
while (abs(IntegPrev - integ) > 0.000000001);

Comments

[tr1ck]

красавчик, спасибо)

[Алексей Павлов]

Вот тут точно ошибки нет?

for (int i = 1; i <= (n - 1); i++)
S += f(a + h*i);

Может быть ошибка, что выполняется всего 3 раза, а n=4. Возможно, нужно указать, что i=0 в инициализации массива.
Уже подзабыл этот метод трапеции, но мне кажется, там делается одним циклом, вот этот цикл for не сильно нужен. Если я правильно помню, вся фигура разбивается на ряд трапеций и находится их сумма, а трапеция - это прямоугольник и прямоугольный треугольник. А вы как-то странно считаете, похоже, отдельно сумму всех прямоугольников, и прибавляете один треугольник, и так несколько раз в цикле do при разных n.

[tr1ck]

Ошибка была с обнулением S и запятой, спасибо) а вы знакомы со scilab?

[Алексей Павлов]

Про те ошибки я понял, я спрашивал, нет ли здесь логической ошибки? Ответ совпадает с расчётным значением?
Scilab я разочек запускал много лет назад. Совсем чуть-чуть в Matlab работал. Вот в MathCAD я много работал в своё время в институте и в аспирантуре.

[tr1ck]

Вроде нет) по формуле.