1.1 --- a/lemon/random.h	Tue Mar 25 16:28:06 2008 +0000
     1.2 +++ b/lemon/random.h	Thu Apr 03 12:06:48 2008 +0000
     1.3 @@ -802,7 +802,7 @@
     1.4  	    nu=V0*std::exp(-xi);
     1.5  	  }
     1.6        } while(nu>std::pow(xi,delta-1.0)*std::exp(-xi));
     1.7 -      return theta*(xi-gamma(int(std::floor(k))));
     1.8 +      return theta*(xi+gamma(int(std::floor(k))));
     1.9      }
    1.10      
    1.11