用重要性抽样法计算积分02/5dxexIx
令xexg,则inifnfEdxxgxfI101,其中2/5xxf,是
以xg分布的随机数,其分布情况可用直接抽样法求得:ln,为0~1上的随机数。
#include
#include
#include
#define N 1000000
void main()
{
double g,I=0;
long int i;
srand((unsigned)time(NULL));
for(i=0;i
g=1.0*rand()/RAND_MAX;
if(g>0) I+=pow(-log(g),2.5)/N;
}
printf("I=%.12f",I);
}
计算结果:I=3.326815065081