Метод Монте-Карло

@Stas12 · Регистрация: 20.10.2011

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

как мне перевести этот код на с++..... помогите пожалуста((

Pascal
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
 program MonteKarlo;
 uses crt;
 Label l1,l2;
 var
 j1,j,a,b,c,n1,k,n:integer;
 I,Y,x:real;
 Begin
 randomize;
 clrscr;
 writeln('Vvod znachenii');
 write('a = ');
 Read(a);
 write('b = ');
 Read(b);
 write('n = ');
 Read(n);
 writeln('--------------------------------');
 writeln('| k | integral | vsego ispitani|');
 for j:=1 to 9 do
 begin
 I:=0;
 for j1:=1 to n do
 begin
 x:=a+(b-a)*random;
 I:=I+x*x+5*x;
 end;
 I:=I*(b-a)/n;
 writeln('--------------------------------');
 writeln('| ',j,' | ',i:2:6,' | ',n,' |');
 {writeln(' Integral = ',i:6:7,' vsego ispitani = ',n,' popalo pod function = ',n1);}
 end;
 writeln('--------------------------------');
 readkey;
 end.

MrGluck · 22.04.2012, 22:36

код кривой и не рабочий

Если "дословно" перевести

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 <iomanip>
#include <conio.h>
#include <ctime>
 
int main()
{
    system("cls");
    srand( time(NULL) );
    int j1 , j, a, b, n, n1=0, k;
    double I, Y, x;
    //std::cout<< "Translit eto ne ok";
    std::cout<< "Input values: \na = ";
    std::cin>> a;
    std::cout<< "b = ";
    std::cin>> b;
    std::cout<< "n = ";
    std::cin>> n;
    std::cout<< "------------------------\n";
    std::cout<< "| k | integral | tests |\n";
    for (j = 1; j <= 9; j++)
    {
        I = 0;
        for (j1 = 1; j1 <= n; j1++)
        {
            x = a + (b - a) * rand()%100/100;
            I += x*x + 5*x;
        }
        I *= (b - a) / n;
        std::cout<< "------------------------\n";
        std::cout<< "| "<< j<< " | "<< std::setw(2)<< std::setprecision(6)<< i<< " | "<< n<< " |\n";
        std::cout<< " Integral = "<< std::setw(6)<< std::setprecision(7)<< i<< "  tests = "<< n<< "  function satisfies = "<< n1;
    }
    std::cout<< "------------------------\n";
    getch();
    return 0;
}

@Stas12 · 22.04.2012, 23:07 **[ТС]**

нашел программу метода монте карло не пойму что означает Zmax

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
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <conio.h>
 
 
 
double f( double x, double y ) {
    return ((x*x+1)/(x*x*x*x-x*x+1));
}
 
int main() {
 
    srand( time ( NULL ) );//чтобы генерированые числа в рандоме не повторялись
   
 
    double Zmax = 0., t, x = 0., y, z;//объявление переменных
    double SumINTEGR = 0., INTEGR = 0., averINTEGR;
    int i, j, HIT;
 
    for ( i = 0; i < 100; ++i ) {
 
        y = 0.;
 
        for ( j = 0; j < 100; ++j ) {
            if( Zmax < ( t = f( x, y ) ) )
                Zmax = t;
            y += 0.01;
        }
        x +=0.01;
    }
 
 
    for( i = 0; i < 10; ++i ) {
 
        // Число попаданий точек в искомый объем под поверхностью функции
        HIT = 0;
 
        //Значение интеграла
        INTEGR = 0.;
 
        for( j = 0; j < 183065; ++j ) {
 
           x = (double) ( rand() ) / RAND_MAX;//заполняем x  в интервале 0 1
           y = (double) ( rand() ) / RAND_MAX;
           z = (double) ( rand() ) * 10.0 / RAND_MAX;
 
            if( z <= f( x, y ) )//проверка
                ++HIT;
 
        }
 
        INTEGR = (double) HIT / 183065. * 10.;
        SumINTEGR += INTEGR;//значение интеграла
    }
 
    // Среднее значение интеграла в 10 случаях
    averINTEGR = SumINTEGR / 10.;
 
    printf( "Значение интеграла   = %lf\n", averINTEGR );
    printf( "Z максимальное       = %lf\n", Zmax );
 getch ();
    return 0;
}

Добавлено через 35 секунд
что это за переменная?? и зачем она вообще нужна?

MrGluck · 22.04.2012, 23:52

Задание вам дали не для того, чтобы вы скопипастили более-менее рабочий код, а чтобы научились программировать. Будут ошибки в вашей программе - обращайтесь.

И используйте наконец теги форматирования кода, проявите уважение к участникам форума.

Новые блоги и статьи Все статьи Все блоги /
Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»	SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 12.02.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .
SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .	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 и по. . .

MrGluck Форумчанин 8216 / 5047 / 1437 Регистрация: 29.11.2010 Сообщений: 13,453
	22.04.2012, 23:52
	Задание вам дали не для того, чтобы вы скопипастили более-менее рабочий код, а чтобы научились программировать. Будут ошибки в вашей программе - обращайтесь. И используйте наконец теги форматирования кода, проявите уважение к участникам форума. 0