# 3.7.jl
# 2024-10-11
# $Id: 3.7.jl 1.1 2024/10/31 12:23:21 s Exp $

using Printf

# input
r    = 0.1
cfsa = [-20000, -1000, -1000, -1000]
ia   = [     0,     1,     2,     3]
cfsb = [-30000, -2000, -2000, -2000, -2000, -2000]
ib   = [     0,     1,     2,     3,     4,     5]

# calculation
pva  = sum(map((x, y) -> x / (1 + r) ^ y, cfsa, ia))
pvb  = sum(map((x, y) -> x / (1 + r) ^ y, cfsb, ib))
na   = length(ia)
nb   = length(ib)
aa   = r * (1 + r) ^ na * -pva / ((1 + r) ^ na - 1)
ab   = r * (1 + r) ^ nb * -pvb / ((1 + r) ^ nb - 1)

# output
println("3.7.jl")
@printf("pv_a= %8d (\$)\n", round(pva))
@printf("A_a = %8d (\$)\n", round(aa ))
@printf("pv_b= %8d (\$)\n", round(pvb))
@printf("A_b = %8d (\$)\n", round(ab ))
# end