সাম-অফ-ইয়ার্স-ডিজিট ও ইউনিটস-অফ-প্রোডাকশন

# সাম-অফ-ইয়ার্স-ডিজিটস (SYD) -- SYD = N(N+1)/2, D_t = [(N-t+1)/SYD] * (B-S)

def syd_depreciation(B, S, N):
    SYD = N * (N + 1) / 2
    schedule = []
    for t in range(1, N + 1):
        D_t = ((N - t + 1) / SYD) * (B - S)
        schedule.append(D_t)
    return schedule

B, S, N = 36000, 6000, 5
syd_schedule = syd_depreciation(B, S, N)
SYD = N * (N + 1) / 2

print(f"B = ${B:,.2f}, S = ${S:,.2f}, N = {N}")
print(f"SYD = N(N+1)/2 = {N}({N}+1)/2 = {SYD:.0f}\n")
print(f"{'বছর':>4} | {'ভগ্নাংশ':>10} | {'D_t':>10}")
print("-" * 30)
for t, D_t in enumerate(syd_schedule, start=1):
    frac = (N - t + 1) / SYD
    print(f"{t:>4} | {frac:>9.4f} | ${D_t:>8,.2f}")

total_syd = sum(syd_schedule)
check_syd = abs(total_syd - (B - S)) < 1e-9
print("-" * 30)
print(f"যাচাই -- সঞ্চিত D == B - S ?  ${total_syd:,.2f} == ${B - S:,.2f} -> {check_syd}")

# ইউনিটস-অফ-প্রোডাকশন -- D_t = (units_t / total_units) * (B - S)

def units_of_production_depreciation(B, S, unit_amounts_per_year):
    total_units = sum(unit_amounts_per_year)
    schedule = []
    for units_t in unit_amounts_per_year:
        D_t = (units_t / total_units) * (B - S)
        schedule.append(D_t)
    return schedule

B, S = 100000, 10000
units_per_year = [50000, 45000, 40000, 35000, 30000]   # মোট আয়ুষ্কাল উৎপাদন = 200,000 ইউনিট
uop_schedule = units_of_production_depreciation(B, S, units_per_year)

print(f"B = ${B:,.2f}, S = ${S:,.2f}, মোট আয়ুষ্কাল উৎপাদন = {sum(units_per_year):,} ইউনিট\n")
print(f"{'বছর':>4} | {'উৎপাদন':>10} | {'D_t':>10}")
print("-" * 32)
for t, (units_t, D_t) in enumerate(zip(units_per_year, uop_schedule), start=1):
    print(f"{t:>4} | {units_t:>10,} | ${D_t:>8,.2f}")

total_uop = sum(uop_schedule)
check_uop = abs(total_uop - (B - S)) < 1e-9
print("-" * 32)
print(f"যাচাই -- সঞ্চিত D == B - S ?  ${total_uop:,.2f} == ${B - S:,.2f} -> {check_uop}")