# 5.12.py
# created on 2016-12-31.
# $Id: 5.12.py 1.2 2016/12/31 11:38:16 s Exp $
#
# df : discount factor
# d : dividend
# r : interest rate
# g : growth rate
# v : value of the company

import numpy as np

# input
n = 5
r = 0.15
g1 = 0.10 # for the next 5 years
g2 = 0.05 # thereafter
d0 = 10000
i  = np.arange(n)

# calculate
d1 = d0 * (1 + g1)
df = 1 / (1 + r) ** (i + 1)
v5 = sum(d1 * (1 + g1) ** i * df)
d5 = d1 * (1 + g1) ** 4
v6 = (1 + g2) * d5 / (r - g2) / (1 + r) ** 5
v0 = v5 + v6

# output
print('5.12(a)')
print("{:5}{:10.3f}".format("r", r))
print("{:5}{:10.3f}".format("g1", g1))
print("{:5}{:10.3f}".format("g2", g2))
print("{:5}{:10.0f}".format("d0", d0))
print("{:5}{:10.0f}".format("d1", d1))
print("df")
for i in range(len(df)):
  print("{:5d}{:10.3f}".format(i+1, df[i]))
print("{:5}{:10.1f}".format("v5", v5))
print("{:5}{:10.1f}".format("d5", d5))
print("{:5}{:10.1f}".format("v6", v6))
print("{:5}{:10.1f}".format("v0", v0))
# end of program