scipy tutorial - Baustatik-Info-Server

SciPy Reference Guide, Release 0.8.dev
Display frozen pmf:
>>> rv = boltzmann(lambda_,N)
>>> x = np.arange(0,np.min(rv.dist.b,3)+1)
>>> h = plt.plot(x,rv.pmf(x))
Check accuracy of cdf and ppf:
>>> prb = boltzmann.cdf(x,lambda_,N)
>>> h = plt.semilogy(np.abs(x-boltzmann.ppf(prb,lambda_,N))+1e-20)
Random number generation:
>>> R = boltzmann.rvs(lambda_,N,size=100)
Custom made discrete distribution:
>>> vals = [arange(7),(0.1,0.2,0.3,0.1,0.1,0.1,0.1)]
>>> custm = rv_discrete(name=’custm’,values=vals)
>>> h = plt.plot(vals[0],custm.pmf(vals[0]))
Boltzmann (Truncated Discrete Exponential)
boltzmann.pmf(k,b,N) = (1-exp(-b))*exp(-b*k)/(1-exp(-b*N)) for k=0,..,N-1
630 Chapter 3. Reference