m********g
发帖数: 46 | 1 how to calculate F(z) numerically by writing our own subroutine without
using R/SAS/Matlab etc? | z*****u
发帖数: 322 | 2 数值积分？
【在 m********g 的大作中提到】
: how to calculate F(z) numerically by writing our own subroutine without
: using R/SAS/Matlab etc?
| m********g
发帖数: 46 | 3 怎么对一个函数做[-infinity, c] 上的数值积分？ | f***a
发帖数: 329 | 4 你积[-999999，c]不就行了
【在 m********g 的大作中提到】
: 怎么对一个函数做[-infinity, c] 上的数值积分？
| m********g
发帖数: 46 | 5 如果求F(-0.005), 假设 dx = 0.001, 这样求解也太慢了。 | s******h
发帖数: 539 | 6 显然, 这个问题等价于只要积[0, c]上面的积分即可，c >0. 比如,
F(-0.005) = 1 - F(0.005) = 1/2 - (F(0.005) - F(0)).
也就是求 \int_0^0.005 e^{-x^2/2} dx (当然不担心前面的相乘的常数1/\sqrt{2\pi}), 裸用泰勒展开就可以算得很准了. 对一般的c,也可以很容易的分成长度为1的一段一段的考虑，也就是说你可以考虑
[0, 1), [1, 2), ..., [[t], t)上面的积分。 |
|
