importnumpyimportpylab# number of customersn=10000# average time between arrivalsa=1# average time to serve each customerb=1# Generate samples:# inter-arrival timesat=numpy.random.exponential(a,n)st=numpy.random.exponential(b,n)foriinrange(n-1):at[i+1]=at[i+1]+at[i]ft=range(n)ft[0]=at[0]+st[0]foriinrange(n-1):ft[i+1]=max(at[i+1]+st[i+1],ft[i]+st[i+1])