Форум программистов, компьютерный форум CyberForum.ru

Ошибка - HIT обнуляются. Метод Монте-Карло для двойного интеграла - C++

Восстановить пароль Регистрация
 
Gmails
5 / 5 / 2
Регистрация: 08.04.2014
Сообщений: 241
25.04.2014, 17:45     Ошибка - HIT обнуляются. Метод Монте-Карло для двойного интеграла #1
пожалуйста помогите найти ошибку.она заключена в HIT,они обнуляются.метод монте-карла для двойного интеграла.
вот код:
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
#include <math.h>
#include <conio.h>
#include <stdio.h>
#include <stdlib.h> 
#include<iostream>
#include <time.h> 
#include<clocale>
#include<math.h>
using namespace std;
 int main(){
    setlocale(LC_CTYPE, "");
    srand((unsigned)time(NULL));
    printf("Функция: z=y*y*sqrt(R*R-x*x)\n");
    float R=3;
    float x,y;
    float Zmax=0;
    for(x=-R;x<=R;x=x+0.01)
    {
    for(y=-sqrt(R*R-x*x);y<=sqrt(R*R-x*x);y=y+0.01)
    {
if(Zmax<(float)(y*y*sqrt(R*R-x*x)))
            {
                Zmax=(float) (y*y*sqrt(R*R-x*x));
            }
        }
    }
    
    float SumINTEGR=0;//Сумма интегралов за 10 опытов
    
    for(int i=1;i<=10;i++)
        {
         float X,Y,Z;// Рандомные точки 
         int HIT=0;// Число попаданий точек в искомый объем под поверхностью функции
          //Значение интеграла
           float INTEGR=0;
            for(float N=1;N<=183065;N++)
            {
                X=(float)(rand()%101)/100;
                Y=(float)(rand()%101)/100;
                Z=(float)(rand()%101)/100;
    
                if(Z<=(y*y*sqrt(R*R-x*x)))
                    {
                        ++HIT;
                    }
            }
            cout<<HIT<<endl;
            INTEGR=(float)HIT/183065*10;
            SumINTEGR=SumINTEGR+INTEGR;
        }
    
    float averINTEGR=SumINTEGR/10;// Среднее значение интеграла в 10 случаях
    printf("Значение интеграла - %f\n", averINTEGR);
    
    printf("Z максимальное - %f",Zmax);
    getch();
 
}
Добавлено через 1 час 54 минуты
up!

Добавлено через 1 час 6 минут
как записать границу ((1-x*x)^(1/2))?

Добавлено через 1 час 25 минут
up!
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
25.04.2014, 17:45     Ошибка - HIT обнуляются. Метод Монте-Карло для двойного интеграла
Посмотрите здесь:

метод Монте-Карло C++
C++ метод Монте-Карло
C++ Метод Монте-Карло
C++ Процедура вычисления кратного интеграла методом Монте-Карло
Решения кратного интеграла методом Монте Карло на С++ C++
C++ Метод монте Карло
Рассчет определенного интеграла методом Монте-Карло. C++
C++ вроде метод монте карло
C++ Вычисление интеграла методом Монте-Карло
C++ Метод Монте-Карло
Метод Монте-Карло C++
Решение определенного интеграла методом Монте Карло C++

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

Текущее время: 01:30. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru