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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Gmails
6 / 6 / 2
Регистрация: 08.04.2014
Сообщений: 248
#1

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

25.04.2014, 17:45. Просмотров 185. Ответов 0
Метки нет (Все метки)

пожалуйста помогите найти ошибку.она заключена в 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++
Всем доброго времени суток. В универе дали задание: вычислить интеграл \int_{0}^{1}dx\int_{0}^{1-x}(1-x-y)dy геометрическим методом...

Решение двойного интеграла - C++
S S ( dxdy/(24+x^2 + y^2)^(1/2) ) G:x^2 + y^2 &lt;=25 G Добавлено через 2 часа 1 минуту Я тут написал программу, но не уверен....

Решение двойного интеграла C++ - C++
Помогите решить двойной интеграл, методом прямоугольников. По матану еще не прошли интегралы, по С++ уже задают, понятия не имею как это...

Написать программу для вычисления двойного факториала числа - C++
2) написать программу для вычисления двойного факториала числа n!!=n*(n-2)*(n-4)*(n-6)... Число n вводится с клавиатуры.

Обнуляются значения вектора - C++
Всем здравствуйте! Столкнулся с очень необычной проблемой, т.к. я до этого работал только на языках, где есть &quot;сборщик мусора&quot;. Есть...

Обнуляются поля в структуре - C++
Всем привет! Курю книжку по плюсам от Шилдта. Домашнее задание было в книге реализовать поиск по имени товара. Вот объявление структуры. ...

Все элементы вектора обнуляются - C++
Делал задачу. Надо в последовательности посчитать сумма между первым и последним нулевым элементами. Для доступа к элементам...

список (в программе обнуляются указатели на структуры) - C++
здравствуйте помогите найти ошибку, при возврате указателя на структуру из функции в функцию майн он обнуляется, это меня сбивает с толку...

Программа для нахождения определенного интеграла - C++
Помогите, написать программу для нахождения определенного интеграла с точностью Е = 10-3 для заданных границ с помощью трех методов: ...


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

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

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