untitled

import numpy as np

# Plutonium-239 half-life in years
half_life = 24100
initial_mass = 1.0 # 1 gram

# Generate time points from 0 to 5 half-lives
t = np.linspace(0, half_life * 5, 500)
mass = initial_mass * (0.5) ** (t / half_life)

# Find specific points for the table
years = [0, half_life, half_life * 2, half_life * 3, half_life * 4, half_life * 5]
remaining_mass = [initial_mass * (0.5) ** (y / half_life) for y in years]

print(list(zip(years, remaining_mass)))