генерация случайных чисел с заданным распределением

[Aladin]

Здравствуйте. Нужна ваша помощь.
Необходим алгоритм или же пример. где реализована генерация случайных чисел с заданным распределением вероятности.

[Rififi]

http://www.boost.org/doc/libs/1_46_1....distributions

[Aladin]

Мне нужно для метода монте-карло, какой из них можно использовать

[MeTeOpA]

Покопайтесь по поводу Randomize() и random. Уверен, что помогут.

[Aladin]

Я делал через random программа считает правильно для простенькой функции

#include "stdafx.h"
#include "iostream"
#include "iomanip"
#include "math.h"

using namespace std;


double PI =  3.141592653589793;

 
double func(double x)
{
   return (x*x+1)/(x*x*x*x+x*x+1); 
}

double D(double min,double max)
{  
	return min + double(rand()) / RAND_MAX * (max - min); 
}


int main()
{
  double a,b,n,k,x,g,res,s=0,res2;
  a=0; b=1;
  cin>>n;
  k=b-a;
  for(int i=1; i<=n; i++)
  {
	  x=D(0,1);
	  s+=func(x);
  }
  res=1/n*k*s;

  res2=PI/(2.0*sqrt(3.0));
  cout<<setprecision(15)<<res<<endl;
  cout<<setprecision(15)<<res2<<endl;
}

[Aladin]

Но если для более сложного уравнения, считает не правильно

#include "stdafx.h" 
#include "iostream" 
#include "iomanip" 
#include "math.h" 

using namespace std; 

double e = 2.71828183; 

double func(double x,double y) 
{ 
        return pow(e,-2*x-2*y)*sin(x+y); 
} 

double DoubleRand(double min,double max) 
    { return min + double(rand()) / RAND_MAX * (max - min); } 

int main() 
{ 
  double a,b,n,k,x,g,res,s=0,res2; 
  double PI = 3.141592653589793; 
  a=0; b=100; 
  cout<<"n = "; cin>>n; 
  k=b-a; 
  double y; 
  for(int i=1; i<=n; i++) 
  { 
      x=DoubleRand(0,100); 
      y=DoubleRand(0,100); 
      s+=func(x,y); 
  } 
  res=1/n*k*s; 

  res2=4.0/25.0; 
  cout<<setprecision(15)<<res<<endl; 
  cout<<setprecision(15)<<res2<<endl; 
}

[MeTeOpA]

Aladin, распишите пожалуйста что конкретно, на ваш взгляд работает неправильно... Не все же знают метод "Монте-Карло"...