Интегрирование методом Монте-Карло

@tazer · Регистрация: 25.04.2013

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

доброго времени суток форумчане, была поставлена задача написать программу которая будет вычислять интеграл методом Монте-Карло
вот что я написал

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 <math.h>
#include <iostream>
#include <conio.h>
 
using namespace std;
 
double F(double x, double y)
{
    return exp(x + y) + exp(x - y)*sin(x);
}
 
double DInt(int n)
{
    int Ax = 0, Bx = 4, Ay = 1, By = 5;
    double S = 0, x, y;
 
    for (int i = 0; i < n; i++)
    {
        x = Ax + (Bx - Ax) * rand() % 2;
        y = Ay + (By - Ay) * rand() % 2;
 
        if ((x >= 0) && (x <= 4) && (y >= 1) && (y <= 5))
            S = F(x, y) + S;
    }
    return (Bx - Ax)*(By - Ay)*S / n;
}
 
void main()
{
    setlocale(LC_ALL, "rus");
    int n;
    cout << "Введите n ";
    cin >> n;
    cout << DInt(n);
    //cout << res;
    system("pause");
}

в итоге ответ получается приблизительно 43 с копейками, начал проверять этот интеграл в Maple и онлайн решалках - ответ приблизительно 7718

где допущена ошибка? подскажите пожалуйста

@zss · 30.11.2015, 21:34

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

x = Ax + (Bx - Ax) * rand() % 2;

C++
1
x = Ax + (Bx - Ax) * (double)rand() /RAND_MAX;

Аналогично y

Новые блоги и статьи Все статьи Все блоги /
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 https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11680&d=1772460536 Одним из. . .	Реалии 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. На борту пять. . .