기본연금액 = 1.2 × (A값 + B값) × (1 + 0.05n/12)

def nps_breakeven(base_amount, delay_years=5):
    """
    국민연금 연기 수령 손익분기점 계산
    base_amount: 정상 수령 시 월 연금액 (만원)
    delay_years: 연기 연수
    """
    enhanced = base_amount * (1 + 0.072 * delay_years)

    # 연기 기간 동안 못 받는 총액
    foregone = base_amount * 12 * delay_years

    # 월 추가 수령액
    monthly_extra = enhanced - base_amount

    # 손익분기점 (월)
    breakeven_months = foregone / monthly_extra
    breakeven_years = breakeven_months / 12

    print(f"정상 수령액: {base_amount}만원/월")
    print(f"연기 후 수령액: {enhanced:.1f}만원/월")
    print(f"연기 기간 포기액: {foregone}만원")
    print(f"손익분기점: {breakeven_years:.1f}년 (수령 시작 후)")
    print(f"=> 수령 시작 나이 + {breakeven_years:.1f}년 이후 생존 시 유리")

    return breakeven_years

# 월 120만원 수령, 5년 연기 가정
nps_breakeven(120, delay_years=5)

import numpy as np

def dc_pension_simulator(
    current_age,
    retirement_age,
    monthly_contribution,
    current_balance,
    annual_return_rate,
    inflation_rate=0.025
):
    """
    DC형 퇴직연금 적립금 시뮬레이션

    Parameters:
        current_age: 현재 나이
        retirement_age: 퇴직 예정 나이
        monthly_contribution: 월 납입액 (만원)
        current_balance: 현재 적립금 (만원)
        annual_return_rate: 연 수익률 (예: 0.06 = 6%)
        inflation_rate: 연 물가상승률
    """
    years = retirement_age - current_age
    monthly_rate = annual_return_rate / 12

    # 현재 잔액의 미래 가치
    future_value_current = current_balance * (1 + annual_return_rate) ** years

    # 월 납입금의 미래 가치 (연금 현가 공식)
    if monthly_rate > 0:
        n = years * 12
        future_value_contributions = monthly_contribution * (
            ((1 + monthly_rate) ** n - 1) / monthly_rate
        )
    else:
        future_value_contributions = monthly_contribution * years * 12

    total_future_value = future_value_current + future_value_contributions

    # 실질 가치 (물가 조정)
    real_value = total_future_value / (1 + inflation_rate) ** years

    print(f"=== DC형 퇴직연금 시뮬레이션 ===")
    print(f"현재 나이: {current_age}세 → 퇴직 나이: {retirement_age}세 ({years}년)")
    print(f"월 납입액: {monthly_contribution}만원")
    print(f"현재 잔액: {current_balance}만원")
    print(f"연 수익률: {annual_return_rate*100:.1f}%")
    print(f"")
    print(f"예상 퇴직 적립금 (명목): {total_future_value:,.0f}만원")
    print(f"예상 퇴직 적립금 (실질): {real_value:,.0f}만원")

    return total_future_value, real_value

# 예시: 35세, 60세 퇴직, 월 30만원 납입, 현재 2,000만원 적립
dc_pension_simulator(
    current_age=35,
    retirement_age=60,
    monthly_contribution=30,
    current_balance=2000,
    annual_return_rate=0.06
)

def irp_tax_credit(annual_contribution, annual_salary):
    """
    IRP + 연금저축 세액공제 계산

    annual_contribution: 연간 납입액 (만원), 최대 900만원 한도
    annual_salary: 연봉 (만원)
    """
    max_deduction = 900  # 만원
    eligible = min(annual_contribution, max_deduction)

    if annual_salary <= 5500:
        rate = 0.165  # 16.5%
    else:
        rate = 0.132  # 13.2%

    tax_credit = eligible * rate

    print(f"연봉: {annual_salary}만원")
    print(f"납입액: {annual_contribution}만원")
    print(f"공제 대상 금액: {eligible}만원")
    print(f"세액공제율: {rate*100}%")
    print(f"연간 세액공제 금액: {tax_credit:.1f}만원")
    print(f"실질 수익률 상당: {tax_credit/annual_contribution*100:.1f}%")

    return tax_credit

# 예시 1: 연봉 4,000만원, 900만원 납입
print("[ 예시 1: 저소득 구간 ]")
irp_tax_credit(900, 4000)

print()

# 예시 2: 연봉 7,000만원, 900만원 납입
print("[ 예시 2: 고소득 구간 ]")
irp_tax_credit(900, 7000)

연금저축: 연 600만원 한도
IRP: 연 300만원 추가 (연금저축과 합산 900만원)

def isa_to_pension_strategy(
    isa_maturity_amount,
    pension_existing_contribution,
    annual_salary
):
    """
    ISA 만기 연금 전환 세액공제 계산

    isa_maturity_amount: ISA 만기 금액 (만원)
    pension_existing_contribution: 기존 연금저축+IRP 납입액 (만원)
    annual_salary: 연봉 (만원)
    """
    # ISA 전환 추가 세액공제 (전환액의 10%, 최대 300만원)
    isa_transfer_credit_base = min(isa_maturity_amount * 0.1, 300)

    # 기존 연금저축 세액공제
    pension_credit_base = min(pension_existing_contribution, 900)

    total_credit_base = pension_credit_base + isa_transfer_credit_base

    if annual_salary <= 5500:
        rate = 0.165
    else:
        rate = 0.132

    pension_tax_credit = pension_credit_base * rate
    isa_transfer_tax_credit = isa_transfer_credit_base * rate
    total_tax_credit = total_credit_base * rate

    print(f"=== ISA 연금 전환 전략 ===")
    print(f"ISA 만기 금액: {isa_maturity_amount}만원")
    print(f"ISA 전환 추가 공제 기준: {isa_transfer_credit_base:.0f}만원")
    print(f"기존 연금 납입 공제 기준: {pension_credit_base}만원")
    print(f"합계 세액공제 기준: {total_credit_base:.0f}만원")
    print(f"")
    print(f"연금저축+IRP 세액공제: {pension_tax_credit:.1f}만원")
    print(f"ISA 전환 추가 세액공제: {isa_transfer_tax_credit:.1f}만원")
    print(f"총 세액공제 합계: {total_tax_credit:.1f}만원")

# 예시: ISA 3,000만원 만기, 기존 연금 900만원, 연봉 5,000만원
isa_to_pension_strategy(3000, 900, 5000)

def pension_tax_calculator(annual_withdrawal, age):
    """
    연금 수령 시 세금 계산 (연금소득세)

    annual_withdrawal: 연간 수령액 (만원)
    age: 수령 시 나이
    """
    # 연금소득세율 (나이에 따라 다름)
    if age < 70:
        tax_rate = 0.055  # 5.5%
    elif age < 80:
        tax_rate = 0.044  # 4.4%
    else:
        tax_rate = 0.033  # 3.3%

    # 분리과세 적용 기준: 연 1,500만원 이하
    THRESHOLD = 1500

    if annual_withdrawal <= THRESHOLD:
        # 분리과세 적용
        tax = annual_withdrawal * tax_rate
        tax_type = "분리과세"
    else:
        # 종합소득세 합산 대상 (세율 개략 추정)
        # 실제로는 다른 소득과 합산하여 누진세율 적용
        basic_tax = THRESHOLD * tax_rate
        excess_tax = (annual_withdrawal - THRESHOLD) * 0.165  # 종소세 구간 가정
        tax = basic_tax + excess_tax
        tax_type = "종합과세"

    after_tax = annual_withdrawal - tax

    print(f"연간 수령액: {annual_withdrawal}만원")
    print(f"나이: {age}세 → 연금소득세율: {tax_rate*100:.1f}%")
    print(f"과세 방식: {tax_type}")
    print(f"세금: {tax:.1f}만원")
    print(f"세후 수령액: {after_tax:.1f}만원")
    print(f"실효세율: {tax/annual_withdrawal*100:.2f}%")

# 예시: 70세, 연 1,200만원 수령
print("[ 70세, 연 1,200만원 수령 ]")
pension_tax_calculator(1200, 70)
print()

# 예시: 65세, 연 2,000만원 수령
print("[ 65세, 연 2,000만원 수령 ]")
pension_tax_calculator(2000, 65)

import numpy as np
import matplotlib.pyplot as plt

def comprehensive_retirement_simulation(
    current_age=35,
    retirement_age=60,
    life_expectancy=90,

    # 현재 자산
    nps_estimated_monthly=80,       # 국민연금 예상 월 수령액 (만원)
    dc_current_balance=1000,        # DC형 퇴직연금 현재 잔액 (만원)
    dc_monthly_contribution=25,     # DC형 월 납입액 (만원)
    pension_savings_balance=500,    # 연금저축 현재 잔액 (만원)
    pension_monthly=50,             # 연금저축+IRP 월 납입액 (만원)

    # 수익률 가정
    dc_return=0.05,                 # DC형 연 수익률
    pension_return=0.06,            # 연금저축 연 수익률
    inflation=0.025,                # 연 물가상승률

    # 은퇴 후 월 필요 생활비 (현재 가치)
    monthly_expense=250             # 만원
):
    """
    종합 노후 자금 시뮬레이션
    """
    accumulation_years = retirement_age - current_age
    distribution_years = life_expectancy - retirement_age

    # 1. DC형 퇴직연금 적립금 계산
    dc_monthly_rate = dc_return / 12
    n_dc = accumulation_years * 12
    dc_future = (dc_current_balance * (1 + dc_return) ** accumulation_years +
                 dc_monthly_contribution * ((1 + dc_monthly_rate) ** n_dc - 1) / dc_monthly_rate)

    # 2. 연금저축+IRP 적립금 계산
    pension_monthly_rate = pension_return / 12
    n_pension = accumulation_years * 12
    pension_future = (pension_savings_balance * (1 + pension_return) ** accumulation_years +
                      pension_monthly * ((1 + pension_monthly_rate) ** n_pension - 1) / pension_monthly_rate)

    # 3. 국민연금 (은퇴 시점의 명목 가치로 변환)
    nps_future_monthly = nps_estimated_monthly * (1 + inflation) ** accumulation_years
    nps_annual = nps_future_monthly * 12

    # 4. 필요 생활비 (미래 가치)
    monthly_expense_future = monthly_expense * (1 + inflation) ** accumulation_years
    annual_expense_future = monthly_expense_future * 12

    # 5. 운용 자산에서 인출 가능 금액 계산 (4% 룰 적용)
    withdrawable_dc = dc_future * 0.04
    withdrawable_pension = pension_future * 0.04

    total_annual_income = withdrawable_dc + withdrawable_pension + nps_annual
    income_gap = annual_expense_future - total_annual_income

    print("=" * 50)
    print("    종합 노후 자금 시뮬레이션 결과")
    print("=" * 50)
    print(f"현재 나이: {current_age}세 → 은퇴: {retirement_age}세 → 기대 수명: {life_expectancy}세")
    print()
    print(f"[은퇴 시점 예상 적립금]")
    print(f"  DC형 퇴직연금: {dc_future:,.0f}만원")
    print(f"  연금저축+IRP: {pension_future:,.0f}만원")
    print(f"  합계 운용 자산: {dc_future+pension_future:,.0f}만원")
    print()
    print(f"[은퇴 후 연간 소득 (4% 인출 기준)]")
    print(f"  DC형 인출: {withdrawable_dc:,.0f}만원/년")
    print(f"  연금저축 인출: {withdrawable_pension:,.0f}만원/년")
    print(f"  국민연금: {nps_annual:,.0f}만원/년")
    print(f"  총 예상 소득: {total_annual_income:,.0f}만원/년")
    print()
    print(f"[은퇴 후 필요 생활비]")
    print(f"  월 필요액 (미래가치): {monthly_expense_future:,.0f}만원")
    print(f"  연 필요액 (미래가치): {annual_expense_future:,.0f}만원")
    print()
    if income_gap > 0:
        print(f"[주의] 연간 소득 부족분: {income_gap:,.0f}만원")
        print(f"=> 추가 저축 또는 지출 조정이 필요합니다.")
    else:
        surplus = -income_gap
        print(f"[양호] 연간 소득 여유분: {surplus:,.0f}만원")
        print(f"=> 현재 계획으로 충분합니다.")

# 35세 직장인 시뮬레이션
comprehensive_retirement_simulation()