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

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

Восстановить пароль Регистрация
Другие темы раздела
C++ Описать логическую функцию equal (l1,l2), проверяющую на равенство списки l1 и l2 http://www.cyberforum.ru/cpp-beginners/thread581551.html
Добрый день, прошу помощи в написании кода (работаю в QT Creator 2.0.1). Задание: Описать логическую функцию equal (l1,l2), проверяющую на равенство списки l1 и l2.
C++ Вывод на экран предложения, содержащие максимальное количество знаков пунктуации Здравстуйте. Задание такое: Написать программу, которая считывает текст из файла и выводит на экран предложения, содержащие максимальное количество знаков пунктуации. Подразумевается, что уже имеется текстовый файл. Пытался сделать программу, но она не работает. Скорее всего неправильно сделал перебор предложений. #include <iostream> #include <fstream> #include <string> using namespace... http://www.cyberforum.ru/cpp-beginners/thread581537.html
C++ Составить программу вычисления с применением условного оператора
u = \begin{cases} {x}^{2}+|x-8|, & x < 8 \\ {x}^{3}+4x-5, & x \geq 8 \end{cases}
Решение по примера C++
ребят не могу написать код под уравнения. помогите плиз... Числа: а, m, n задаются через int
C++ Составить программу вычисления с применением условного оператора http://www.cyberforum.ru/cpp-beginners/thread581518.html
u=\begin{cases} {x}^{2}-3, & x < 0 \\ x, & x \geq 0 \end{cases}
C++ Подсчет суммы с заданной точностью Подсчет суммы с заданной точностью "e>0" #include "stdafx.h" #include <iostream> #include <cmath> using namespace std; int main() { int n, i=1; подробнее

Показать сообщение отдельно
unt
0 / 0 / 0
Регистрация: 21.05.2012
Сообщений: 8
21.05.2012, 17:59     Процедура вычисления кратного интеграла методом Монте-Карло
Создать функцию для вычисления кратных интегралов методом Монте Карло.
Помогите написать рабочий код. Укажите на ошибки.

Процедура вычисления кратного интеграла методом Монте-Карло
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
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
#include <stdio.h>
#include <time.h>
#include <math.h>
#include <conio.h>
#include <stdlib.h>
 
double func(double a, double b, double c)
{
    return (a*b*c)/sqrt(pow(a,2)+pow(b,2)+pow(c,2));
}
 
int debug(int x)
{
    if((x>-32768) && (x<32767));
    else
   {
    printf("\n \t \t Vy vveli nevernyj simvol. \n \t \t Neobhodimo bylo vvesti celoe chislo. \n \n");
    getch();
    exit(1);
    }
    return 0;
}
 
void main()
{
    double x, y, z;
   double a, b, c;
   unsigned int n;
   int i;
   double sum=0;
   double S=0;
   double min=0;
 
   srand(time(NULL));
    printf("\n \t Vvediteniz. granicu integrirovanija (x): ");
    scanf("%lf",&x);
   debug(x);
 
    printf("\t Vvedite granicu integrirovanija (y): ");
    scanf("%lf",&y);
   debug(y);
 
   printf("\t Vvedite  granicu integrirovanija (z): ");
    scanf("%lf",&z);
   debug(z);
 
 
    inv: printf("\n \t Vvedite chislio ispytanij \n");
    printf("\t N = ");
    scanf("%d",&n);
   if((n<0) || (n>65535))
   {
    printf("\n \t Vy vveli nepravilnoe chislo, kolichestvo mozhet byt ot 0 do 65535 \n");
    goto inv;
   }
 
 
    printf("\n \t / / / \n");
    printf("\t | | | (x*y*z)/sqr(pow(x,2)+pow(y,2)+pow(z,2))dxdydz \n");
    printf("\t / / / \n");
   printf("\t   D   \n");
 
   // вычисления двойного интеграла
   for(i=0;i<n;i++)
   {
 
    a=x*rand()/RAND_MAX;//генерируем x в интервале
      b=y*rand()/RAND_MAX;//генерируем y в интервале
      c=z*rand()/RAND_MAX;//генерируем z в интервале
 
      if (sqrt(pow(a,2)+pow(b,2)+pow(c,2))<=1)
    {
           sum+=func(a, b, c);
        }
         printf("\n \n \t a,b,c = %8.3lf %8.3lf %8.3lf\n", a,b,c);
   }
   S=(((x-0)*(y-0)*(z-0))/n)*sum;
   printf("\n \n \t Integral = %8.3lf \n", S);
 
   getch();
}
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 00:33. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru