#!/usr/bin/env wolframscript
(* ::Package:: *)

LogPlot[Table[PDF[HypoexponentialDistribution[Table[2/k,{k,1,n}]],x],{n,{5,7,10}}],{x,0,40}]


nbins = 50;
nlam = 10000;
nrealizations = 100000;
R1 = 1.
gama = 2.
ncoll = Range[1, nlam, 1]
lam = R1*ncoll^gama
data=RandomVariate[HypoexponentialDistribution[lam],nrealizations];
Export["hypoexpoentialD_W.dat", data]
Show[Histogram[data,nbins,"PDF"]]








