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)