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

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

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

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

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
#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

Считает:

Code
1
2
3
4
$ ./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 **[ТС]**

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

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .
Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .	SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .

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

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

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

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

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

Решение