ইকোনমিক সার্ভিস লাইফ

# ইকোনমিক সার্ভিস লাইফ -- প্রতিটি সম্ভাব্য রিটেনশন N-এর জন্য AW-of-total-cost গণনা করে সর্বনিম্ন খোঁজা

def a_p_factor(i, n):
    return i * (1 + i) ** n / ((1 + i) ** n - 1)

def a_f_factor(i, n):
    return i / ((1 + i) ** n - 1)

P = 40000     # ক্রয়মূল্য
i = 0.10      # সুদের হার

# ইনডেক্স 0 = বছর ১, ইনডেক্স 5 = বছর ৬
salvage = [32000, 26000, 21000, 17000, 14000, 12000]   # N বছর রাখলে স্যালভেজ ভ্যালু
om      = [3000, 3800, 4700, 5700, 6900, 8300]         # সেই বছরে হওয়া O&M খরচ

max_n = len(salvage)
aw_by_n = {}

print(f"{'N':>3} | {'Capital Recovery':>18} | {'AW of O&M':>14} | {'AW মোট':>14}")
print("-" * 58)

for n in range(1, max_n + 1):
    S_n = salvage[n - 1]

    # capital recovery with salvage -- P*(A/P,i,n) - S_n*(A/F,i,n)
    capital_recovery = P * a_p_factor(i, n) - S_n * a_f_factor(i, n)

    # n বছর পর্যন্ত প্রতিটি O&M খরচকে প্রথমে প্রেজেন্ট ওয়ার্থে আনা, তারপর n বছরের উপর AW-তে রূপান্তর
    pw_om = sum(om[t - 1] / (1 + i) ** t for t in range(1, n + 1))
    aw_om = pw_om * a_p_factor(i, n)

    aw_total = capital_recovery + aw_om
    aw_by_n[n] = aw_total

    print(f"{n:>3} | ${capital_recovery:>16,.2f} | ${aw_om:>12,.2f} | ${aw_total:>12,.2f}")

# সর্বনিম্ন AW খুঁজে বের করা -- অনুমান নয়, সার্চ করে
best_n = min(aw_by_n, key=aw_by_n.get)
print("-" * 58)
print(f"\nইকোনমিক সার্ভিস লাইফ: N = {best_n} বছর (সর্বনিম্ন AW = ${aw_by_n[best_n]:,.2f})")