Метод средних прямоугольников и количество вычислений для вычисления определенного интеграла (Си -> C++)

@kosetik · Регистрация: 22.07.2015

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

Задание:Вычислить значение определенного интеграла методом средних прямоугольников и вычислить количество разбиений , потребовавшееся для вычисления интеграла

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
38
39
40
41
42
43
#include <stdio.h>
#include <conio.h>
 
void main()
{
 clrscr();
 float a,b;   // границы отрезка
 float dx;   // приращение аргумента
 float s;   // приближенное значение интеграла
 int n;    // количество интервалов
 float x;   // аргумент
 float y;    // значение функции в начале интервала
 int i;
 
 printf ("Приближенное вычесление интеграла! (f(x)=5x^2-x+2)\n");
 printf ("Метод прямоугольников\n");
 printf ("\nВведите:\n");
 
 printf ("Нижняя граница интервала:\n");
 printf ("=>");
 scanf  ("%f",&a);
 
 printf ("Верхняя граница интервала:\n");
 printf ("=>");
 scanf  ("%f",&b);
 
 printf ("Приращение аргумента:\n");
 printf ("=>");
 scanf  ("%f",&dx);
 
 n=(b-a)/dx+1;
 x=a;
 s=0;
 
 for (i=1; i<=n; i++)
 {
  y=x*x+2;   // значение функции в начале интервала
  s+=y*dx;
  x+=dx;
 }
 printf ("\n\nЗначение интеграла: %6.3f",s);
 getch();
}

@Серж762 · 03.04.2018, 16:05

Мой вариант:

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
#include <iostream>
#include <cmath>
 
using namespace std;
 
double f(double x)
{
    return 5 * x * x - x + 2;
}
 
int main()
{
    double a,b,h;
    int n;
    double INTGRL = 0.0;
 
    cout << "Vvedite a= ";
    cin >> a;
 
    cout << "Vvedite b= ";
    cin >> b;
 
    cout << "Vvedite h= ";
    cin >> h;
 
    n = floor(((b-a)/h));
 
    for (double x = a; x <= b; x+=h)
        INTGRL += f(x-h/2);
    INTGRL*=h;
 
    cout << "Kol-tvo otrezkov = " << n << endl;
    cout << endl;
    cout << "S = " << INTGRL;
    return 0;
}

Результат:

Кликните здесь для просмотра всего текста

Метод средних прямоугольников и количество вычислений для вычисления определенного интеграла (Си -> C++)

@chessman2 · 05.06.2024, 21:19

Я не задаю точность, но задаю количество прямоугольников, кратных 2.
Начинаю с 8, вычисляю, потом удваиваю: 16, 32, 64, 128.
Каждое удвоение дает еще одну точную цифру после запятой (это примерно, и зависит от пределов интегрирования.
128 дает хорошую точность.

Код не на си к сожалению, но он понятный.

Code
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
\ Вычисление интеграла f(x) = x*x методом прямоугольников. Вещественные числа.
 
: f_x ( F: f1 -- f2 ) f**2 ;
 
\ N - количество прямоугольников, A, B - пределы интегрирования
== N  128
== A  1
== B  5
variable H
variable X
 
: integral ( F: -- sum ) 
    B s>f A s>f f- N s>f f/ H f!    \ --                     ; h = (b - a)/n
    A s>f H f@ f# 2.0 f/ f- X f!    \ --                      ; x0 = a - h/2
    fzero N 0                             \ ( F: sum ) N 0              
    DO                            \ ( F: sum )
        X f@ H f@ f+ X f!         \ ( F: sum )           ; x = x + h
        X f@                      \ ( F: sum x )
        f_x                               \ ( F: sum y )
        H f@ f* f+                \ ( F: sum' ) 
   LOOP ;
 
: main
    cr integral ." integral = " f. cr ;;

Выхлоп:

integral = 41.333008

@Pphantom · 05.06.2024, 21:47

chessman2, забавно, в первый раз за сколько-то лет вижу на форумах код на FORTH.

@chessman2 · 05.06.2024, 22:07

Забавно, что на таком почтенном форуме нет ветки forth

Вообще - не совсем форт. Это самодельный транслятор на flat assembler,
реализующий двухстековую машину. Да, лексика 90% фортовская.

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL3_image 8Observer8 10.02.2026 Содержание блога Библиотека SDL3_image содержит инструменты для расширенной работы с изображениями. Пошагово создадим проект для загрузки изображения формата PNG с альфа-каналом (с прозрачным. . .	Установка Qt-версии Lazarus IDE в Debian Trixie Xfce volvo 10.02.2026 В общем, достали меня глюки IDE Лазаруса, собранной с использованием набора виджетов Gtk2 (конкретно: если набирать текст в редакторе и вызвать подсказку через Ctrl+Space, то после закрытия окошка. . .	SDL3 для Web (WebAssembly): Работа со звуком через SDL3_mixer 8Observer8 08.02.2026 Содержание блога Пошагово создадим проект для загрузки звукового файла и воспроизведения звука с помощью библиотеки SDL3_mixer. Звук будет воспроизводиться по клику мышки по холсту на Desktop и по. . .	SDL3 для Web (WebAssembly): Основы отладки веб-приложений на SDL3 по USB и Wi-Fi, запущенных в браузере мобильных устройств 8Observer8 07.02.2026 Содержание блога Браузер Chrome имеет средства для отладки мобильных веб-приложений по USB. В этой пошаговой инструкции ограничимся работой с консолью. Вывод в консоль - это часть процесса. . .
SDL3 для Web (WebAssembly): Обработчик клика мыши в браузере ПК и касания экрана в браузере на мобильном устройстве 8Observer8 02.02.2026 Содержание блога Для начала пошагово создадим рабочий пример для подготовки к экспериментам в браузере ПК и в браузере мобильного устройства. Потом напишем обработчик клика мыши и обработчик. . .	Философия технологии iceja 01.02.2026 На мой взгляд у человека в технических проектах остается роль генерального директора. Все остальное нейронки делают уже лучше человека. Они не могут нести предпринимательские риски, не могут. . .	SDL3 для Web (WebAssembly): Вывод текста со шрифтом TTF с помощью SDL3_ttf 8Observer8 01.02.2026 Содержание блога В этой пошаговой инструкции создадим с нуля веб-приложение, которое выводит текст в окне браузера. Запустим на Android на локальном сервере. Загрузим Release на бесплатный. . .	SDL3 для Web (WebAssembly): Сборка C/C++ проекта из консоли 8Observer8 30.01.2026 Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .

@Pphantom 2281 / 1540 / 716 Регистрация: 17.03.2022 Сообщений: 4,947
	05.06.2024, 21:47
	chessman2, забавно, в первый раз за сколько-то лет вижу на форумах код на FORTH. 0

@chessman2 161 / 141 / 10 Регистрация: 21.10.2012 Сообщений: 480
	05.06.2024, 22:07
	Забавно, что на таком почтенном форуме нет ветки forth Вообще - не совсем форт. Это самодельный транслятор на flat assembler, реализующий двухстековую машину. Да, лексика 90% фортовская. 1

Метод средних прямоугольников и количество вычислений для вычисления определенного интеграла (Си -> C++)

Решение