def GammaDist(alpha, beta, x):
if x <= 0:
return -100000
else:
if alpha == 0 or beta == 0:
return -1000000
else:
return pow(x,alpha-1) * special.exp(-1.0*x / beta) / (pow(beta,alpha)*special.gamma(alpha))
def prob_old(c, n, j, k):
pi = k - j # how many new marbles to pick
po = c - pi # how many old marbles to pick
lnp = factln(n - c) + factln(j) + factln(n - j)
lnp -= factln(n) + factln(j - po) + factln(n - j - pi)
p = exp(lnp)
return p
def F_gauss_sine(t):
dt = 10.0
omega = 1.5
return A*exp(-(t-40)**2/dt**2)*sin(t*omega)
def F_gauss(t):
dt = 5.0
return A*exp(-(t-30)**2/dt**2)
def F_gauss_sine(t):
dt = 10.0
omega = 1.5
return A * exp(-(t - 40)**2 / dt**2) * sin(t * omega)
def F_gauss(t):
dt = 5.0
return A * exp(-(t - 30)**2 / dt**2)
