def radial_wave_func(n,l,r): import math from math import factorial as ff a = c.physical_constants['Bohr radius'][0] #a = float(r)**(-3/2.0) AA = 2.0 / (n*a) AA = AA**3 #print AA B = ff(n-l-1) Ca = 2*n*((ff(n+l))**3) BB= float(B)/Ca #print BB sqre = math.sqrt(AA*BB) #print sqre expp = math.exp(float(-r)/(n*a)) Lexp = ((2*float(r))/(n*a))**l pfunc = assoc_laguerre(2*l+1, n-l-1) y = pfunc(((2*r)/(n*r))) return round(sqre*expp*Lexp*y / a**(-3.0/2),5)
def angular_wave_func(m,l,theta,phi): if m >0: eps = (-1)**(m) else: eps = 1 a = (2*l + 1)/(4*math.pi) cca = ff(l-abs(m)) ccb = ff(l+abs(m)) b = float(cca)/float(ccb) sqr = math.sqrt(a*b) #print sqr expon = math.exp(1j*m*float(phi)) #print expon pfunc = assoc_legendre(m,l) y = pfunc(float(theta)) #print(y) final = eps*sqr*expon*y return complex(round(final.real,5), round(final.imag,5))
def f(N): return 10**100*ff(N-500)*ff(N-70)/ff(N-567)/ff(N)
print("you are still kind of young") else: print("enjoy your drink") age_check(23) ''' '''#for days = ("Mon","Tue","Wed","Thur","Fri") for day in days: if day == "Wed": break else: print(day) for letter in "Kimjae": print(letter) ''' #module ''' import math #math 모듈 전체를 가져옴 -비추 print(math.ceil(1.2)) ''' from math import ceil, fsum #math 모듈중 원하는 함수만 가져옴 from math import floor as ff #함수 이름 원하는대로 바꿔주기 print(ceil(2.6)) print(fsum([1, 2, 3, 4, 5, 5.9])) print(ff(4.3))
def cc(n, r): return ff(n) // ff(n - r) // ff(r)
def question_02(x,y): result = np.ones((x,y), dtype=np.int) # Takes care of first row and first col for xi in range(1,x): for yi in range(1,y): result[xi,yi] = result[xi-1,yi] + result[xi,yi-1] return result def num_paths(x,y): # list of x items, each of which are lists of length y paths = [[1 for _ in range(y+1)] for _ in range(x+1)] # Skip over first row and first col as we know those are are 1. for row in range(1,x+1): for col in range(1,y+1): paths[row][col] = paths[row-1][col]+paths[row][col-1] print(np.array(paths)) return paths[x][y] if __name__ == "__main__": #grid = question_02(2,2) #print(grid) #print("") #grid = question_02(15,10) #print(grid) print("\n\n2017\n\n") print(num_paths(2,2)) print(num_paths(15,10)) print( ff(15+10)/(ff(15)*ff(10)) )
# ex02.module.py # 방법 1 import threading import time import os import threading as th # OR import threading, time, os, sys, math # 모듈 안에 있는 모든 메소드를 메모리에 올려 놓는다. # 방법2 from math import factorial from math import factorial as fct # 모듈의 특정 메소드만 가져온다 # 방법3 from math import * # 1번과 동일 # 방법4 from math import (factorial as ff, acos as ac) n = fct(5) / fct(3) print(n) n2 = ff(5) + ac(1) print(n2)