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

@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

Новые блоги и статьи Все статьи Все блоги /
Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ФедосеевПавел 06.01.2026 Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ВВЕДЕНИЕ Введу сокращения: аналоговый ПИД — ПИД регулятор с управляющим выходом в виде числа в диапазоне от 0% до. . .	Модель микоризы: классовый агентный подход 2 anaschu 06.01.2026 репозиторий https:/ / github. com/ shumilovas/ fungi ветка по-частям. коммит Create переделка под биомассу. txt вход sc, но sm считается внутри мицелия. кстати, обьем тоже должен там считаться. . . .	Расчёт токов в цепи постоянного тока igorrr37 05.01.2026 / * Дана цепь постоянного тока с сопротивлениями и напряжениями. Надо найти токи в ветвях. Программа составляет систему уравнений по 1 и 2 законам Кирхгофа и решает её. Последовательность действий:. . .	Новый CodeBlocs. Версия 25.03 palva 04.01.2026 Оказывается, недавно вышла новая версия CodeBlocks за номером 25. 03. Когда-то давно я возился с только что вышедшей тогда версией 20. 03. С тех пор я давно снёс всё с компьютера и забыл. Теперь. . .
Модель микоризы: классовый агентный подход anaschu 02.01.2026 Раньше это было два гриба и бактерия. Теперь три гриба, растение. И на уровне агентов добавится между грибами или бактериями взаимодействий. До того я пробовал подход через многомерные массивы,. . .	Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Programma_Boinc 28.12.2025 Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Налог на собак: https:/ / ********/ gallery/ V06K53e Финансовый отчет в Excel: https:/ / ********/ gallery/ bKBkQFf Пост отсюда. . .	Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Programma_Boinc 26.12.2025 Нашел на реддите интересную статью под названием Anyone know where to get a free Desktop or Laptop? Ниже её машинный перевод. После долгих разбирательств я наконец-то вернула себе. . .	Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Programma_Boinc 23.12.2025 Рецензия / Мнение/ Перевод Нашел на реддите интересную статью под названием The Thinkpad X220 Tablet is the best budget school laptop period . Ниже её машинный перевод. Thinkpad X220 Tablet —. . .