В общем облазил и-нет, форум, но в своём алгоритме ошибку так и не вижу, может кто увидит и прекратит мои страдания? Проблема в самом алгоритме, извращения с классами из-за задания, так что не судите по ним. Решение интеграла методом Монте-Карло.
| 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 <iostream>
#include <time.h>
int const top=6;
int const bottom=3;
using namespace std;
class A
{
private:
int a;
public:
A() {a;}
void seta() {a=bottom;}
int geta() {return a;}
};
class B
{
private:
int b;
public:
B() {b;}
void setb() {b=top;}
int getb() {return b;}
};
class Integr:public B,public A
{
private:
int N,count,x,y,max,min;
public:
Integr() {N=200;}
void solve()
{
count=0;
setb();
seta();
if (geta()>getb()) {max=func(geta()); min=func(getb());}
else {max=func(getb()); min=func(geta());}
for (int i=0;i<N;i++)
{
x=geta()+rand()%(getb()-geta()+1);
y=min+rand()%(max-min+1);
if (y<=func(x)) {count++;}
}
cout << (max-min)*(getb()-geta())*count/N << endl;
}
int func(int handle)
{
return handle*handle*handle/3;
}
};
void main()
{
srand(time(NULL));
Integr t;
t.solve();
system("pause");
} |
|
(
