다음은 블랙-리터만 기대수익률 공식을 옮긴 것이다. 공식이 길어서 한 번에 쓰는 대신에 bl1, bl2, bl3, bl4로 나눠 작성한 후 다시 합쳐 연산하고 pi_adj에 저장한 것이다. 대부분의 연산이 내적곱이다.
공식이 길면 작성하기도 어렵지만, 결괏값이 이상해 틀린 부분을 찾으려 할 때도 불리하다. 이럴 때는 공식을 크게 나눠 각 부분의 이상 여부를 확인하는 것이 편리하다.
# 투자자 전망과 합쳐진 균형초과수익률 계산
bl1 = inv( np.dot( tau, C ) )
bl2 = np.dot( np.dot( np.transpose( P ), inv( omega ) ), P )
bl3 = np.dot( inv( np.dot( tau, C ) ), pi )
bl4 = np.dot( np.dot( np.transpose( P ), inv( omega ) ), Q )
pi_adj = np.dot( inv( bl1+bl2 ), ( bl3+bl4 ) )