Программа перестает работать при смене интервала

@valentinaaa · Регистрация: 05.02.2014

Author24 — интернет-сервис помощи студентам

почему, когда я меняю интервал с 0 - 6,2832 программа не работает?

C

#include <stdio.h>
#include <math.h>
 
double f1(double x) {
    return (sin(x) / (x * x - x + 1));
}
 
int main() {
    int k, 
    int n = 2;
    double a = 0, b = 6.28, i1 = 0., i2, h;
    double eps = 0.00001;
 
    i2 = (b - a) * (f1(a) + f1(b)) / 2;
 
    while (fabs(i2 - i1) > eps) {
        n = n * 2;
        h = (b - a) / n;
        i1 = i2;
        i2 = 0;
 
        for (k = 0; k < n; k++) {
            i2 = i2 + h * f1(a + (k + .5)) * h;
        }
    }
    printf("f(x)=sin(x)/(x*x-x+1) \n");
    printf("Интервал с %f по %f \n", a, b);
    printf("Значение определенного интеграла it %.2f \n", i2);
    
}

07.03.2014, 22:42

что конкретно не так?

@valentinaaa · 07.03.2014, 22:46 **[ТС]**

когда в 11 строчке меняю с 6,28 на 6,2832 программа перестает считать интегралл

@Black Fregat · 07.03.2014, 22:58

Считает:

Код

$ ./integral
f(x)=sin(x)/(x*x-x+1) 
Интервал с 0.000000 по 6.283200 
Значение определенного интеграла it 0.00

Что именно у Вас происходит? Какая ошибка?

@valentinaaa · 07.03.2014, 23:01 **[ТС]**

значение интеграла должно быть 1,1..., не могу понять почему, когда я меняю интервал он его не считает и как исправить, чтоб работал при 6,2832? ведь при 6,28 он считает верно

@Black Fregat · 08.03.2014, 00:09

Просто у Вас неудачное начальное значение i1.
При изменении пределов начальное значение i2 становится слишком близким к нулю
и цикл вылетает до начала первой итерации.

Поставьте, например, i1 = -1;

А меня ввела в заблуждение ошибка в строке 23.
У Вас стоит

C++

1	i2 = i2 + h * f1(a + (k + .5)) * h;

а нужно

C++

1	i2 = i2 + h * f1(a + (k + .5) * h);

А с этой ошибкой получается нуль в обоих случаях.
Вот я и не понял, что именно сломалось..

@valentinaaa · 08.03.2014, 11:38 **[ТС]**

спасибо огромное, теперь все работает

id_slk
	07.03.2014, 22:42	2
	что конкретно не так?

@valentinaaa 0 / 0 / 0 Регистрация: 05.02.2014 Сообщений: 6
	07.03.2014, 22:46 [ТС]	3
	когда в 11 строчке меняю с 6,28 на 6,2832 программа перестает считать интегралл 0

@Black Fregat Фрилансер 3705 / 2077 / 567 Регистрация: 31.05.2009 Сообщений: 6,683
	07.03.2014, 22:58	4
	Считает: Код $ ./integral f(x)=sin(x)/(x*x-x+1) Интервал с 0.000000 по 6.283200 Значение определенного интеграла it 0.00 Что именно у Вас происходит? Какая ошибка? 0

@valentinaaa 0 / 0 / 0 Регистрация: 05.02.2014 Сообщений: 6
	07.03.2014, 23:01 [ТС]	5
	значение интеграла должно быть 1,1..., не могу понять почему, когда я меняю интервал он его не считает и как исправить, чтоб работал при 6,2832? ведь при 6,28 он считает верно 0

@valentinaaa 0 / 0 / 0 Регистрация: 05.02.2014 Сообщений: 6
	08.03.2014, 11:38 [ТС]	7
	спасибо огромное, теперь все работает 0

Программа перестает работать при смене интервала

Решение