02.12.2012, 19:00. Просмотров 616. Ответов 5
Написал программу, которая вычисляет интеграл методом Монте-Карло. Вроде как код проверял, перепроверял и должно быть всё верно. Но при указывании кол-ва точек >7000 программа не работает, то есть зависает. А если, к примеру, указать кол-во точек 10000, то выбивает с ошибкой "Divide error". С чем это связано не пойму.
| 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
| #include <iostream.h>
#include <conio.h>
#include <stdlib.h>
#include <math.h>
#include <stdio.h>
class method{
public:
double allsumm(double *x,int n);
void createrand(double *x,int n);
};
double method::allsumm(double *x,int n){
double temp=0;
for(int i=0;i<n;i++){
temp+=sqrt(pow(x[i],4)-pow(x[i],3)+8);
}
return temp;
}
void method::createrand(double *x,int n){
for(int i=0;i<n;i++){
x[i]=0.0001*(rand()%30001);
}
}
void main(){
clrscr();
double a=0,b=3.2,temp;
unsigned long int n;
double *x;
method MonteKarlo;
cout<<"Vvedite koli4estvo to4ek ";
cin>>n;
x=new double[n];
/*cout<<"Vvedite interval\n a: ";
cin>>a;
cout<<"b: ";
cin>>b;*/
MonteKarlo.createrand(x,n);
temp=MonteKarlo.allsumm(x,n);
cout<<"Summ= "<<temp;
cout<<"\nI= "<<((b-a)*temp)/n;
getch();
} |
|
Так всё сложно?
