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)