Интеграл методом Симпсона

@Mee · Регистрация: 05.07.2010

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

Вообщем на форуме есть такая програма, спасибо её автору:

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
ouble Simpson(double (*Fx)(double))
{
        double I,x,h,a,b;
        int m;
        cout << "\n Input integral atributes:\n" << endl;
        cout << "\t-> Enter a = ";
        cin >> a;
        cout << "\t-> Enter b = ";
        cin >> b;
        cout << "\t-> Enter m = ";
        cin >> m;
        h=(b-a)/(m-1);
        x = a;
        I = Fx(a) + Fx(b);
        int n=0;
        while (n < m-2)
        {
                x = x+h;
                if (n%2==0) I=I+4*Fx(x);
                        else  I = I+2*Fx(x);
                n++;
        }
        return I*h/3;
}

Но здесь m - число отрезков, на которые разбивается основой, если я правильно понял, она влияет на точность, а вот так задавать точность мне кажется неправельным. Помогите переделать так, что бы интеграл находился не по количеству разбитых отрезков, а ориентировался на точность eps например 10е-15, которая например передается аргументом. Или как в случае, если есть точность, нати m?

Евгений М. · 30.11.2010, 08:54

Думаю переменная h - то что нужно.

@Mee · 30.11.2010, 19:57 **[ТС]**

Сообщение от Евгений М.

Думаю переменная h - то что нужно.

А как найти h если мы не знаем m? я что-то не понял к чему это

Евгений М. · 01.12.2010, 03:38

Мы его не находим, а задаем. Далее либо мы через h, a, b находим m и не будем переделывать часть кода с циклом, либо переделать часть кода с циклом, чтобы не использовать m. Например сделать цикл от a+h до b-h.

@impuls__s · 13.02.2011, 18:04

а подскажите пожалуйста где можно найти эту программу целиком ? очень надо ...

Евгений М. · 13.02.2011, 21:51

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

где можно найти

В интернете.

Новые блоги и статьи Все статьи Все блоги /
Управление камерой с помощью скрипта 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-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .	Реалии 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. . .

Евгений М. 1080 / 1007 / 107 Регистрация: 28.02.2010 Сообщений: 2,889
	30.11.2010, 08:54
	Думаю переменная h - то что нужно. 0

Евгений М. 1080 / 1007 / 107 Регистрация: 28.02.2010 Сообщений: 2,889
	01.12.2010, 03:38
	Мы его не находим, а задаем. Далее либо мы через h, a, b находим m и не будем переделывать часть кода с циклом, либо переделать часть кода с циклом, чтобы не использовать m. Например сделать цикл от a+h до b-h. 0

@impuls__s 0 / 0 / 0 Регистрация: 10.02.2010 Сообщений: 37
	13.02.2011, 18:04
	а подскажите пожалуйста где можно найти эту программу целиком ? очень надо ... 0