===================================================================
RCS file: /cvs/ext/extinction_risk_1.R,v
retrieving revision 1.8
retrieving revision 1.9
diff -u -p -r1.8 -r1.9
--- ext/extinction_risk_1.R	2015/08/02 12:41:02	1.8
+++ ext/extinction_risk_1.R	2015/08/04 13:24:42	1.9
@@ -1,5 +1,5 @@
 # extinction_risk_1.R
-# $Id: extinction_risk_1.R,v 1.8 2015/08/02 12:41:02 hako Exp $
+# $Id: extinction_risk_1.R,v 1.9 2015/08/04 13:24:42 hako Exp $
 #
 # Author: Hiroshi Hakoyama <hako@affrc.go.jp>
 # Copyright (c) 2013-2015 Hiroshi Hakoyama <hako@affrc.go.jp>, All rights reserved.
@@ -95,7 +95,7 @@ ext1 <- function(dat, t = 100, ne = 1, alpha = 0.05, v
     if(z < 35) {
       pnorm(-w) + exp((z^2 - w^2) / 2) * pnorm(-z)
     } else {
-      pnorm(-w) + exp(- w^2 / 2) * (1 / sqrt(pi)) * (1 / z - 1 / z^3 + 3 / z^5 - 15 / z^7 + 15 * 7 / z^9 - 15 * 7 * 9 / z^11 + 15 * 7 * 9 * 11 / z^13)
+      pnorm(-w) + exp(- w^2 / 2) * (sqrt(2) / (2 * sqrt(pi))) * (1 / z - 1 / z^3 + 3 / z^5 - 15 / z^7 + 15 * 7 / z^9 - 15 * 7 * 9 / z^11 + 15 * 7 * 9 * 11 / z^13)
     }
   }