球一段现成的normcdf C/C++ code - 未名空间MITBBS历史存档

国际科技财经博客移民网络热点娱乐民生时事公众号

Redian新闻

>未名空间

>PhotoGear - 摄影器材

球一段现成的normcdf C/C++ code

球一段现成的normcdf C/C++ code# PhotoGear - 摄影器材

a92011-03-13 08:03

1 楼

一条河中有鳄鱼,专门咬男人的 JJ , 一个女人和二十个男人同时过河,要怎么才能保证
每个人的安全呢?

C*c2011-03-13 08:03

2 楼

哪位大湿有？
google上搜到的有的是要收费，有的精度不够

m*u2011-03-13 08:03

3 楼

女的第一个, 男的串糖葫芦

【在 a9 的大作中提到】

: 一条河中有鳄鱼,专门咬男人的 JJ , 一个女人和二十个男人同时过河,要怎么才能保证
: 每个人的安全呢?

S*M2011-03-13 08:03

4 楼

你搜一下numerical recipe, incomplete gamma function吧
我只有fortran的code, 你用不上

【在 C****c 的大作中提到】

: 哪位大湿有？
: google上搜到的有的是要收费，有的精度不够

e*y2011-03-13 08:03

5 楼

没有女的, 男的串个圈也行哪

【在 m*****u 的大作中提到】

: 女的第一个, 男的串糖葫芦

o*62011-03-13 08:03

6 楼

这样可以吗？
double f(double x)
{
return exp(-x*x/2);
}
double F(double b=0.0,double ep=1e-6) // calculate normal CDF using
integration
{
double h,s1=0,s2=b*(1+f(b))/2;
int n;
for(n=1;fabs(s1-s2)>ep;n*=2)
{
h=b/n;
s1 = s2;
s2 = 0;
for(int k=0;k{
s2 += h*f((k+0.5)*h);
}
s2 = (s1+s2)/2;
}
return s2*sqrt(1/(8*atan(1.0)))+0.5;
}

l*e2011-03-13 08:03

7 楼

e

【在 e***y 的大作中提到】

: 没有女的, 男的串个圈也行哪

c*y2011-03-13 08:03

8 楼

太麻烦，如果精度在10e-6的话，用这个近似就足够了
where
b0 = 0.2316419, b1 = 0.319381530, b2 = 0.356563782, b3 = 1.781477937, b4
= 1
.821255978, b5 = 1.330274429

using

【在 o*******6 的大作中提到】

: 这样可以吗？
: double f(double x)
: {
: return exp(-x*x/2);
: }
: double F(double b=0.0,double ep=1e-6) // calculate normal CDF using
: integration
: {
: double h,s1=0,s2=b*(1+f(b))/2;
: int n;

T12011-03-13 08:03

9 楼

为了避免串糖葫芦，每次由女的背一个男的过河，女的再游回来。
如果附加条件是不能单人过河，就只能串糖葫芦了
不过如果可以串葫芦，这个女的就是多此一举了

o*62011-03-13 08:03

10 楼

贪版真乃牛人也！

781477937, b4

【在 c********y 的大作中提到】

: 太麻烦，如果精度在10e-6的话，用这个近似就足够了
: where
: b0 = 0.2316419, b1 = 0.319381530, b2 = 0.356563782, b3 = 1.781477937, b4
: = 1
: .821255978, b5 = 1.330274429
:
: using

T12011-03-13 08:03

11 楼

没有女的，男的也不一定非要串个圈

【在 e***y 的大作中提到】

: 没有女的, 男的串个圈也行哪

c*y2011-03-13 08:03

12 楼

牛啥啊，
我以前组里一个a3嫌python里的慢，喜欢用这个。
不过精度要求高的时候，就歇菜了。。。

【在 o*******6 的大作中提到】

: 贪版真乃牛人也！
:
: 781477937, b4