Не могу найти ошибку в программе, вычисляющей определенный интеграл по методу трапеций

@Миха_ил · Регистрация: 24.08.2015

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

Написал программу для вычисления интеграла от x^(1/3) в пределах от -1 до 8. На выходе выдает нули и бесконечности

Подскажите плиз в чем дело. Использовал тот факт, что f(x)=f_n+((f_n+1 - f_n)*(x-x_n))/(x_n+1-x_n) А интеграл считал как J=dx*(f₀+f_n)/2+dx*(f₁+f₂+...+f_n-1])

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
31
32
33
34
35
36
37
#include<iostream>
#include<conio.h>
#include<cmath>
#include<fstream>
using namespace std;
 
int main()
{
    ofstream F("C:\\Users\\user\\Documents\\Visual Studio 2013\\Projects\\Методы мат_моделирования\\Second_lab\\2.rtf", ios::in);
    setlocale(LC_ALL,"rus");
 
        double x = 0, y = 0, J, dx ,x1, y0, yn, k;
    int a = 0, b = 9, n = 100;
    dx = (b - a) / 100;
    x1 = dx;
    y0 = pow(dx,1/3);
    yn = pow(9,1/3);
    k = (y0 + yn)*dx / 2;
    J = k;
    for (int i = 0; i < n; i++)
    {
        x += dx;
        x1 += x;
        y += pow(x, 1 / 3) + ((pow(x1, 1 / 3) - pow(x, 1 / 3))*(9 - x)) / (x1 - x);
    }
 
    for (int i = 1; i < n; i++)
    {
        J += pow(x,1/3)*dx;
        printf("%5f %5f %5f\n", x, y, J);
        F << x << "   " << y << "   " << J << endl;
    }
 
 
    _getch();
    return 0;
}

@IrineK · 15.09.2015, 02:28

В строке 14 dx будет равен 0.
Сделайте так:

C++
1
dx = (b - a) / 100.0;

Добавлено через 2 минуты
Похожий сюжет с 1/3
1.0/3 лучше.

Добавлено через 5 минут
Что-то в строке 23 не того.

Новые блоги и статьи Все статьи Все блоги /
Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .	SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .
Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .	Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.