?? poidev.cpp
字號(hào):
double poidev(int xm, int &idum)
{
double pi,g,em,t,sq,alxm,y,dum,temp;
int oldm=0;
pi = 3.141592654;
if (xm < 12)
{
if (xm != oldm)
{
oldm = xm;
g = exp(-xm);
}
em = -1.0;
t = 1.0;
do
{
em = em + 1.0;
t = t * ran1(idum);
}
while (t > g);
}
else
{
if (xm != oldm)
{
oldm = xm;
sq =sqrt(2.0 * xm);
alxm = log(xm);
g = xm * alxm - gammln(xm + 1.0);
}
do
{
do
{
y = tan(pi * ran1(idum));
em = sq * y + xm;
}
while (em < 0.0);
em = int(em);
dum = em * alxm - gammln(em + 1.0) - g;
t = 0.9 * (1.0 + y *y) * exp(dum);
}
while (ran1(idum) > t);
}
temp = em;
return temp;
}
?? 快捷鍵說(shuō)明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -