# prices를 numpy matrix로 변환
prices = np.matrix( prices )

# 수익률 행렬을 만들어 계산한다
rows, cols = prices.shape
returns = np.empty( [ rows, cols – 1 ] )
for r in range( rows ):
    for c in range( cols - 1 ):
        p0, p1 = prices[ r, c ], prices[ r, c + 1 ]
        returns[ r, c ] = ( p1 / p0 ) - 1

# 수익률 계산
expreturns = np.array( [ ] )
for r in range( rows ):
    expreturns = np.append( expreturns, np.mean( returns[ r ] ) )

# 공분산 계산
covars = np.cov( returns )
R = ( 1 + expreturns ) ** 250 - 1 # 연율화
C = covars * 250 # 연율화

# 무위험이자율
rf = .015