রেট অফ রিটার্ন ও ব্রেকইভেনে ইনফ্লেশনের প্রভাব

# নমিনাল IRR (বিসেকশন) -> ফিশার সম্পর্ক দিয়ে রিয়েল IRR

def npv(rate, cash_flows):
    return sum(cf / (1 + rate) ** t for t, cf in enumerate(cash_flows))

def find_irr(cash_flows, lo=-0.5, hi=1.0, tol=1e-9, max_iter=200):
    npv_lo = npv(lo, cash_flows)
    npv_hi = npv(hi, cash_flows)
    assert npv_lo * npv_hi < 0, "bracket-এ সাইন বদল নেই -- রেঞ্জ বাড়ান"
    mid, npv_mid = lo, npv_lo
    for _ in range(max_iter):
        mid = (lo + hi) / 2
        npv_mid = npv(mid, cash_flows)
        if abs(npv_mid) < tol:
            break
        if npv_lo * npv_mid < 0:
            hi = mid
        else:
            lo, npv_lo = mid, npv_mid
    return mid, npv_mid

def real_rate(i_combined, f):
    return (i_combined - f) / (1 + f)

# --- প্রজেক্ট: আজ $20,000 খরচ, ৫ বছর ধরে বছরে $5,500 (actual dollars) আয় ---
cash_flows = [-20000] + [5500] * 5

nominal_irr, converged_npv = find_irr(cash_flows)
print(f"নমিনাল/actual-dollar IRR = {nominal_irr*100:.4f}%")
print(f"কনভার্জড NPV @ IRR       = ${converged_npv:.6f}  (প্রায় শূন্য -- বিসেকশন সঠিকভাবে কনভার্জ করেছে)\n")

f = 0.06   # প্রত্যাশিত ইনফ্লেশন
real_irr = real_rate(nominal_irr, f)
print(f"ইনফ্লেশন f            = {f*100:.2f}%")
print(f"ফিশার সম্পর্ক থেকে রিয়েল IRR = {real_irr*100:.4f}%")
print(f"-> নমিনাল IRR দেখতে ভালো মনে হলেও, ইনফ্লেশন বাদ দিলে প্রকৃত মার্জিন অনেক পাতলা")

# ডিফারেনশিয়াল ইনফ্লেশন -- SP বছরে ৩% বাড়ে, VC বছরে ৭% বাড়ে -- ব্রেকইভেন শিফট

FC = 50000       # ফিক্সড কস্ট (নমিনাল, স্থির ধরা হয়েছে -- মূল ফোকাস মার্জিন-শিফটে)
SP0 = 40.0       # আজকের প্রতি-ইউনিট বিক্রয়মূল্য
VC0 = 25.0       # আজকের প্রতি-ইউনিট ভেরিয়েবল কস্ট
sp_inflation = 0.03
vc_inflation = 0.07

print(f"FC = ${FC:,.0f}, SP0 = ${SP0:.2f}, VC0 = ${VC0:.2f}")
print(f"SP ইনফ্লেশন = {sp_inflation*100:.0f}%/yr, VC ইনফ্লেশন = {vc_inflation*100:.0f}%/yr\n")
print(f"{'বছর':>4} | {'SP':>8} | {'VC':>8} | {'মার্জিন (SP-VC)':>16} | {'Q_BE':>10}")
print("-" * 56)
for t in range(0, 6):
    sp_t = SP0 * (1 + sp_inflation) ** t
    vc_t = VC0 * (1 + vc_inflation) ** t
    margin_t = sp_t - vc_t
    q_be_t = FC / margin_t
    print(f"{t:>4} | ${sp_t:>6.2f} | ${vc_t:>6.2f} | ${margin_t:>14.2f} | {q_be_t:>9,.1f}")