আফটার-ট্যাক্স ক্যাশ ফ্লো অ্যানালাইসিস

# আফটার-ট্যাক্স ক্যাশ ফ্লো (ATCF) টেবিল -- BTCF, D, taxable income, tax, ATCF
# P/A ফ্যাক্টর M3-এ প্রতিষ্ঠিত হয়েছে -- এখানে আফটার-ট্যাক্স PW বের করতে পুনর্ব্যবহার করা হলো

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

cost = 40_000          # ইকুইপমেন্ট কস্ট, বছর ০ (মূলধনী বিনিয়োগ -- সরাসরি deductible নয়)
salvage = 0
life = 5
annual_net = 15_000    # প্রতি বছর রেভিনিউ - খরচ (BTCF, বছর ১-৫)
tax_rate = 0.25
marr = 0.08             # আফটার-ট্যাক্স MARR

D = (cost - salvage) / life   # স্ট্রেইট-লাইন ডেপ্রিসিয়েশন (L30 রিভিশন)

print(f"{'বছর':>4} | {'BTCF':>12} | {'D':>10} | {'ট্যাক্সেবল ইনকাম':>18} | {'ট্যাক্স':>10} | {'ATCF':>12}")
print("-" * 80)

rows = [(0, -cost, 0.0, 0.0, 0.0, -cost)]
print(f"{0:>4} | {-cost:>12,.2f} | {0.0:>10,.2f} | {0.0:>18,.2f} | {0.0:>10,.2f} | {-cost:>12,.2f}")

for t in range(1, life + 1):
    btcf = annual_net
    ti = btcf - D
    tax = tax_rate * ti
    atcf = btcf - tax
    rows.append((t, btcf, D, ti, tax, atcf))
    print(f"{t:>4} | {btcf:>12,.2f} | {D:>10,.2f} | {ti:>18,.2f} | {tax:>10,.2f} | {atcf:>12,.2f}")

total_d = sum(r[2] for r in rows)
assert abs(total_d - (cost - salvage)) < 1e-9, "ডেপ্রিসিয়েশন টোটাল মিলছে না"
print(f"\nচেক: মোট ডেপ্রিসিয়েশন = ${total_d:,.2f}  (cost - salvage = ${cost - salvage:,.2f}) -- মিলেছে")

atcf_annual = rows[1][5]
pw = -cost + atcf_annual * pa_factor(marr, life)

print(f"\nবার্ষিক ATCF (বছর ১-৫, uniform) : ${atcf_annual:,.2f}")
print(f"আফটার-ট্যাক্স PW @ {marr*100:.0f}% MARR       : ${pw:,.2f}")
verdict = "গ্রহণযোগ্য (PW > 0)" if pw > 0 else "প্রত্যাখ্যান (PW <= 0)"
print("সিদ্ধান্ত:", verdict)