equal
deleted
inserted
replaced
793 { |
793 { |
794 xi=1.0-std::log(V1); |
794 xi=1.0-std::log(V1); |
795 nu=V0*std::exp(-xi); |
795 nu=V0*std::exp(-xi); |
796 } |
796 } |
797 } while(nu>std::pow(xi,delta-1.0)*std::exp(-xi)); |
797 } while(nu>std::pow(xi,delta-1.0)*std::exp(-xi)); |
798 return theta*(xi-gamma(int(std::floor(k)))); |
798 return theta*(xi+gamma(int(std::floor(k)))); |
799 } |
799 } |
800 |
800 |
801 /// Weibull distribution |
801 /// Weibull distribution |
802 |
802 |
803 /// This function generates a Weibull distribution random number. |
803 /// This function generates a Weibull distribution random number. |
817 ///\param k shape parameter (<tt>k>0</tt>) |
817 ///\param k shape parameter (<tt>k>0</tt>) |
818 ///\param x_min location parameter (<tt>x_min>0</tt>) |
818 ///\param x_min location parameter (<tt>x_min>0</tt>) |
819 /// |
819 /// |
820 double pareto(double k,double x_min) |
820 double pareto(double k,double x_min) |
821 { |
821 { |
822 return exponential(gamma(k,1.0/x_min)); |
822 return exponential(gamma(k,1.0/x_min))+x_min; |
823 } |
823 } |
824 |
824 |
825 /// Poisson distribution |
825 /// Poisson distribution |
826 |
826 |
827 /// This function generates a Poisson distribution random number with |
827 /// This function generates a Poisson distribution random number with |