ইনক্রিমেন্টাল রেট অফ রিটার্ন অ্যানালাইসিস

# ইনক্রিমেন্টাল রেট অফ রিটার্ন অ্যানালাইসিস

def npv(cashflows, i):
    return sum(cf / (1 + i) ** t for t, cf in cashflows)

def find_irr(cashflows, lo=-0.5, hi=1.0, tol=1e-7, max_iter=200):
    f_lo = npv(cashflows, lo)
    f_hi = npv(cashflows, hi)
    if f_lo * f_hi > 0:
        raise ValueError("bracket-এ NPV চিহ্ন বদলাচ্ছে না")
    for _ in range(max_iter):
        mid = (lo + hi) / 2
        f_mid = npv(cashflows, mid)
        if abs(f_mid) < tol or (hi - lo) / 2 < tol:
            return mid
        if f_lo * f_mid < 0:
            hi = mid
        else:
            lo = mid
            f_lo = f_mid
    return (lo + hi) / 2

MARR = 0.10

# Alt A: খরচ $10,000, বার্ষিক আয় $4,000, ৪ বছর
cf_A = [(0, -10000), (1, 4000), (2, 4000), (3, 4000), (4, 4000)]
# Alt B: খরচ $20,000, বার্ষিক আয় $7,000, ৪ বছর
cf_B = [(0, -20000), (1, 7000), (2, 7000), (3, 7000), (4, 7000)]
# ইনক্রিমেন্টাল ক্যাশ ফ্লো -- বড় (B) বিয়োগ ছোট (A), বিনিয়োগ অনুযায়ী ছোট থেকে বড় সাজানো
cf_incremental = [(t, b - a) for (t, a), (_, b) in zip(cf_A, cf_B)]

irr_A = find_irr(cf_A)
irr_B = find_irr(cf_B)
irr_inc = find_irr(cf_incremental)

print(f"Alt A: IRR_A = {irr_A*100:.2f}%   (NPV চেক = ${npv(cf_A, irr_A):,.6f})")
print(f"Alt B: IRR_B = {irr_B*100:.2f}%   (NPV চেক = ${npv(cf_B, irr_B):,.6f})")
print("ইনক্রিমেন্টাল (B - A) ক্যাশ ফ্লো:", cf_incremental)
print(f"ইনক্রিমেন্টাল IRR = {irr_inc*100:.2f}%   (NPV চেক = ${npv(cf_incremental, irr_inc):,.6f})")
print()
print(f"MARR = {MARR*100:.0f}%")

if irr_inc >= MARR:
    print("সিদ্ধান্ত: ইনক্রিমেন্টাল IRR MARR-এর চেয়ে বেশি -- বাড়তি বিনিয়োগ যুক্তিসঙ্গত, তাই Alt B বাছাই করা উচিত।")
else:
    print("সিদ্ধান্ত: ইনক্রিমেন্টাল IRR MARR-এর চেয়ে কম -- বাড়তি বিনিয়োগ যুক্তিসঙ্গত নয়, তাই Alt A-তেই থাকা উচিত।")