Mercurial > forge

changeset 5680:a01b1d68712a octave-forge

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Replaced expm1.
author asnelt
date Sat, 23 May 2009 20:10:57 +0000
parents 2b962ee1530a
children 3876d698aca1
files main/statistics/inst/copulapdf.m
diffstat 1 files changed, 1 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/main/statistics/inst/copulapdf.m	Sat May 23 13:32:29 2009 +0000
+++ b/main/statistics/inst/copulapdf.m	Sat May 23 20:10:57 2009 +0000
@@ -137,7 +137,7 @@
       if (d != 2)
         error ("copulapdf: Frank copula PDF implemented as bivariate only");
       endif
-      p = (theta .* exp (theta .* (1 + sum (x, 2))) .* expm1 (theta))./ (exp (theta) - exp (theta + theta .* x(:, 1)) + exp (theta .* sum (x, 2)) - exp (theta + theta .* x(:, 2))) .^ 2;
+      p = (theta .* exp (theta .* (1 + sum (x, 2))) .* (exp (theta) - 1))./ (exp (theta) - exp (theta + theta .* x(:, 1)) + exp (theta .* sum (x, 2)) - exp (theta + theta .* x(:, 2))) .^ 2;
       # Product copula at columns where theta == 0
       k = find (theta == 0);
       if (any (k))