how-dense-are-neutron-stars

import%20marimo%0A%0A__generated_with%20%3D%20%220.17.6%22%0Aapp%20%3D%20marimo.App(width%3D%22medium%22%2C%20app_title%3D%22how-dense-are-neutron-stars%22)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_()%3A%0A%20%20%20%20%23%20Created%20by%20Shanika%20Galaudage%0A%0A%20%20%20%20%23%20The%20production%20of%20this%20notebook%20was%20aided%20by%0A%20%20%20%20%23%20Claude%20AI%20for%20setting%20this%20up%20with%20marimo.%0A%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20import%20numpy%20as%20np%0A%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%0A%20%20%20%20plt.rcParams.update(%7B%0A%20%20%20%20%20%20%20%20%22font.family%22%3A%20%22serif%22%2C%0A%20%20%20%20%20%20%20%20%22font.serif%22%3A%20%5B%22Computer%20Modern%22%2C%20%22CMU%20Serif%22%2C%20%22DejaVu%20Serif%22%5D%2C%0A%20%20%20%20%20%20%20%20%22font.size%22%3A%209.0%2C%0A%20%20%20%20%20%20%20%20%22mathtext.fontset%22%3A%20%22cm%22%2C%0A%20%20%20%20%20%20%20%20%22axes.formatter.use_mathtext%22%3A%20False%2C%0A%20%20%20%20%20%20%20%20%22axes.labelsize%22%3A%2011%2C%0A%20%20%20%20%20%20%20%20%22axes.titlesize%22%3A%2011%2C%0A%20%20%20%20%20%20%20%20%22axes.labelpad%22%3A%203.0%2C%0A%20%20%20%20%20%20%20%20%22axes.grid%22%3A%20True%2C%0A%20%20%20%20%20%20%20%20%22grid.color%22%3A%20%22black%22%2C%0A%20%20%20%20%20%20%20%20%22grid.linestyle%22%3A%20(0%2C%20(1%2C%204))%2C%0A%20%20%20%20%20%20%20%20%22grid.linewidth%22%3A%200.55%2C%0A%20%20%20%20%20%20%20%20%22grid.alpha%22%3A%200.4%2C%0A%20%20%20%20%20%20%20%20%22figure.subplot.left%22%3A%200.17%2C%0A%20%20%20%20%20%20%20%20%22figure.subplot.right%22%3A%200.97%2C%0A%20%20%20%20%20%20%20%20%22figure.subplot.bottom%22%3A%200.14%2C%0A%20%20%20%20%20%20%20%20%22figure.subplot.top%22%3A%200.98%2C%0A%20%20%20%20%20%20%20%20%22legend.frameon%22%3A%20True%2C%0A%20%20%20%20%20%20%20%20%22legend.framealpha%22%3A%200.75%2C%0A%20%20%20%20%20%20%20%20%22legend.edgecolor%22%3A%20%22none%22%2C%0A%20%20%20%20%20%20%20%20%22legend.borderpad%22%3A%200.15%2C%0A%20%20%20%20%20%20%20%20%22legend.labelspacing%22%3A%200.3%2C%0A%20%20%20%20%20%20%20%20%22legend.handletextpad%22%3A%200.5%2C%0A%20%20%20%20%20%20%20%20%22legend.handlelength%22%3A%201.6%2C%0A%20%20%20%20%20%20%20%20%22legend.handleheight%22%3A%200.25%2C%0A%20%20%20%20%20%20%20%20%22legend.fontsize%22%3A%209%2C%0A%20%20%20%20%20%20%20%20%22lines.linewidth%22%3A%201.0%2C%0A%20%20%20%20%20%20%20%20%22xtick.major.pad%22%3A%204%2C%0A%20%20%20%20%20%20%20%20%22ytick.major.pad%22%3A%202.5%2C%0A%20%20%20%20%20%20%20%20%22xtick.labelsize%22%3A%2010%2C%0A%20%20%20%20%20%20%20%20%22ytick.labelsize%22%3A%2010%2C%0A%20%20%20%20%20%20%20%20%22xtick.major.size%22%3A%203.0%2C%0A%20%20%20%20%20%20%20%20%22ytick.major.size%22%3A%203.0%2C%0A%20%20%20%20%20%20%20%20%22xtick.minor.size%22%3A%201.5%2C%0A%20%20%20%20%20%20%20%20%22ytick.minor.size%22%3A%201.5%2C%0A%20%20%20%20%20%20%20%20%22xtick.direction%22%3A%20%22in%22%2C%0A%20%20%20%20%20%20%20%20%22ytick.direction%22%3A%20%22in%22%2C%0A%20%20%20%20%20%20%20%20%22savefig.dpi%22%3A150%2C%0A%20%20%20%20%7D)%0A%20%20%20%20return%20mo%2C%20np%2C%20plt%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%20How%20much%20*space*%20is%20all%20of%20humanity%3F%0A%0A%20%20%20%20Neutron%20stars%20pack%20an%20enormous%20amount%20of%20mass%20into%20a%20sphere%20only%0A%20%20%20%20about%20the%20size%20of%20a%20city.%0A%0A%20%20%20%20So%20here's%20a%20fun%20question%3A%0A%0A%20%20%20%20%3E%20If%20you%20compressed%20**all%20of%20humanity**%20to%20the%20**average%20density%20of%20a%0A%20%20%20%20%3E%20neutron%20star**%2C%20how%20much%20space%20would%20we%20take%20up%3F%0A%0A%20%20%20%20We'll%20look%20at%20this%20two%20ways%3A%20first%20for%20**every%20human%20who%20has%20ever%20lived**%2C%0A%20%20%20%20and%20then%20for%20**everyone%20alive%20today**.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20%23%20The%20mass%20of%20the%20Sun%2C%20in%20kilograms.%0A%20%20%20%20%23%20Neutron%20star%20masses%20are%20quoted%20as%20multiples%20of%20this%20(%22solar%20masses%22).%0A%20%20%20%20M_SUN%20%3D%201.989e30%20%20%23%20kg%0A%0A%20%20%20%20%23%20Preset%20stars%3A%20(mass%20in%20solar%20masses%2C%20equatorial%20radius%20in%20km).%0A%20%20%20%20%23%20SOURCE%3A%20https%3A%2F%2Fui.adsabs.harvard.edu%2Fabs%2F2026ApJ...1005..201K%2Fabstract%0A%20%20%20%20CUSTOM%20%3D%20%22Custom%20(drag%20a%20slider%20or%20type%20a%20value)%22%0A%20%20%20%20presets%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22PSR%20J0030%2B0451%20%E2%80%94%20Kini%20et%20al.%202026%20(median)%22%3A%20(1.43%2C%2012.68)%2C%0A%20%20%20%20%20%20%20%20%22Canonical%20Mass%20and%20Classic%20Radius%20Value%22%3A%20(1.40%2C%2010.00)%2C%0A%20%20%20%20%20%20%20%20CUSTOM%3A%20None%2C%0A%20%20%20%20%7D%0A%20%20%20%20star_dropdown%20%3D%20mo.ui.dropdown(%0A%20%20%20%20%20%20%20%20options%3Dlist(presets)%2C%20value%3Dlist(presets)%5B0%5D%2C%20label%3D%22Neutron%20Star%3A%22%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20Shared%20state%3A%20each%20mass%20widget%20(and%20each%20radius%20widget)%20reads%20and%20writes%0A%20%20%20%20%23%20the%20same%20value%2C%20so%20dragging%20a%20slider%20updates%20its%20box%20and%20vice-versa.%0A%20%20%20%20get_mass%2C%20set_mass%20%3D%20mo.state(1.43)%0A%20%20%20%20get_radius%2C%20set_radius%20%3D%20mo.state(12.68)%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20CUSTOM%2C%0A%20%20%20%20%20%20%20%20M_SUN%2C%0A%20%20%20%20%20%20%20%20get_mass%2C%0A%20%20%20%20%20%20%20%20get_radius%2C%0A%20%20%20%20%20%20%20%20presets%2C%0A%20%20%20%20%20%20%20%20set_mass%2C%0A%20%20%20%20%20%20%20%20set_radius%2C%0A%20%20%20%20%20%20%20%20star_dropdown%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(get_mass%2C%20mo%2C%20set_mass)%3A%0A%20%20%20%20mass_slider%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%201.10%2C%202.20%2C%20step%3D0.01%2C%20value%3Dget_mass()%2C%20on_change%3Dset_mass%2C%0A%20%20%20%20%20%20%20%20label%3D%22Mass%20(solar%20masses%2C%20M_sun)%3A%22%2C%20show_value%3DTrue%2C%0A%20%20%20%20)%0A%20%20%20%20return%20(mass_slider%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(get_mass%2C%20mo%2C%20set_mass)%3A%0A%20%20%20%20mass_box%20%3D%20mo.ui.number(%0A%20%20%20%20%20%20%20%201.10%2C%202.20%2C%20step%3D0.01%2C%20value%3Dget_mass()%2C%20on_change%3Dset_mass%2C%0A%20%20%20%20%20%20%20%20label%3D%22or%20type%3A%22%2C%0A%20%20%20%20)%0A%20%20%20%20return%20(mass_box%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(get_radius%2C%20mo%2C%20set_radius)%3A%0A%20%20%20%20radius_slider%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%209.00%2C%2015.00%2C%20step%3D0.01%2C%20value%3Dget_radius()%2C%20on_change%3Dset_radius%2C%0A%20%20%20%20%20%20%20%20label%3D%22Radius%20(km)%3A%22%2C%20show_value%3DTrue%2C%0A%20%20%20%20)%0A%20%20%20%20return%20(radius_slider%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(get_radius%2C%20mo%2C%20set_radius)%3A%0A%20%20%20%20radius_box%20%3D%20mo.ui.number(%0A%20%20%20%20%20%20%20%209.00%2C%2015.00%2C%20step%3D0.01%2C%20value%3Dget_radius()%2C%20on_change%3Dset_radius%2C%0A%20%20%20%20%20%20%20%20label%3D%22or%20type%3A%22%2C%0A%20%20%20%20)%0A%20%20%20%20return%20(radius_box%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20---%0A%0A%20%20%20%20%23%23%201.%20Choose%20a%20neutron%20star%0A%0A%20%20%20%20First%2C%20we%20need%20a%20neutron%20star%20whose%20average%20density%20we%20can%20calculate.%0A%0A%20%20%20%20Pick%20one%20below%2C%20or%20choose%20your%20own%20mass%20and%20radius.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20CUSTOM%2C%0A%20%20%20%20mass_box%2C%0A%20%20%20%20mass_slider%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20presets%2C%0A%20%20%20%20radius_box%2C%0A%20%20%20%20radius_slider%2C%0A%20%20%20%20star_dropdown%2C%0A)%3A%0A%20%20%20%20controls%20%3D%20%5Bstar_dropdown%5D%0A%20%20%20%20if%20star_dropdown.value%20%3D%3D%20CUSTOM%3A%0A%20%20%20%20%20%20%20%20controls%20%2B%3D%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.hstack(%5Bmass_slider%2C%20mass_box%5D%2C%20justify%3D%22start%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.hstack(%5Bradius_slider%2C%20radius_box%5D%2C%20justify%3D%22start%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20m%2C%20r%20%3D%20presets%5Bstar_dropdown.value%5D%0A%20%20%20%20%20%20%20%20controls.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(r%22%24M%20%3D%20%25.2f%5C%20M_%5Codot%20%5Cqquad%20r%20%3D%20%25.2f%5C%20%5Ctext%7Bkm%7D%24%22%20%25%20(m%2C%20r))%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20mo.vstack(controls)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(CUSTOM%2C%20M_SUN%2C%20get_mass%2C%20get_radius%2C%20np%2C%20presets%2C%20star_dropdown)%3A%0A%20%20%20%20import%20math%0A%0A%20%20%20%20if%20star_dropdown.value%20%3D%3D%20CUSTOM%3A%0A%20%20%20%20%20%20%20%20star_solar_masses%20%3D%20get_mass()%0A%20%20%20%20%20%20%20%20star_radius_km%20%3D%20get_radius()%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20star_solar_masses%2C%20star_radius_km%20%3D%20presets%5Bstar_dropdown.value%5D%0A%0A%20%20%20%20%23%20Convert%20units%20and%20compute%20the%20average%20density.%0A%20%20%20%20star_mass_kg%20%3D%20star_solar_masses%20*%20M_SUN%20%20%20%20%20%20%20%20%20%20%20%20%23%20solar%20masses%20-%3E%20kg%0A%20%20%20%20star_radius_cm%20%3D%20star_radius_km%20*%201e5%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20km%20-%3E%20cm%0A%20%20%20%20star_volume_cm3%20%3D%20(4.%20%2F%203.)%20*%20np.pi%20*%20star_radius_cm**3%0A%20%20%20%20density_kg_per_mL%20%3D%20star_mass_kg%20%2F%20star_volume_cm3%20%20%23%201%20cm%5E3%20%3D%201%20mL%0A%0A%20%20%20%20def%20sci(x%2C%20sig%3D4)%3A%0A%20%20%20%20%20%20%20%20%23%20Format%20a%20number%20as%20LaTeX%20scientific%20notation%0A%20%20%20%20%20%20%20%20exp%20%3D%20int(math.floor(math.log10(abs(x))))%0A%20%20%20%20%20%20%20%20mant%20%3D%20x%20%2F%2010.0**exp%0A%20%20%20%20%20%20%20%20return%20(r%22%25.*f%20%5Ctimes%2010%5E%7B%25d%7D%22)%20%25%20(sig%20-%201%2C%20mant%2C%20exp)%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20density_kg_per_mL%2C%0A%20%20%20%20%20%20%20%20sci%2C%0A%20%20%20%20%20%20%20%20star_mass_kg%2C%0A%20%20%20%20%20%20%20%20star_radius_cm%2C%0A%20%20%20%20%20%20%20%20star_radius_km%2C%0A%20%20%20%20%20%20%20%20star_solar_masses%2C%0A%20%20%20%20%20%20%20%20star_volume_cm3%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20M_SUN%2C%0A%20%20%20%20density_kg_per_mL%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20sci%2C%0A%20%20%20%20star_mass_kg%2C%0A%20%20%20%20star_radius_cm%2C%0A%20%20%20%20star_radius_km%2C%0A%20%20%20%20star_solar_masses%2C%0A%20%20%20%20star_volume_cm3%2C%0A)%3A%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20r%22%22%22%0A%20%20%20%20%20%20%20%20---%0A%0A%20%20%20%20%20%20%20%20%23%23%202.%20How%20dense%20is%20it%3F%0A%0A%20%20%20%20%20%20%20%20We'll%20treat%20the%20neutron%20star%20as%20a%20sphere%20and%20calculate%20its%0A%20%20%20%20%20%20%20%20**average%20density**%20from%20its%20mass%20and%20radius.%0A%0A%20%20%20%20%20%20%20%20To%20make%20the%20final%20density%20easier%20to%20compare%20with%20everyday%20amounts%2C%0A%20%20%20%20%20%20%20%20we'll%20work%20in%20**kilograms**%20for%20mass%20and%20**cubic%20centimetres**%20for%20volume.%0A%0A%20%20%20%20%20%20%20%20First%2C%20convert%20the%20mass%20from%20solar%20masses%20to%20kilograms%3A%0A%0A%20%20%20%20%20%20%20%20%24M%20%3D%20%25.2f%5C%2CM_%5Codot%0A%20%20%20%20%20%20%20%20%3D%20%25.2f%20%5Ctimes%20%25s%5C%20%5Ctext%7Bkg%7D%0A%20%20%20%20%20%20%20%20%3D%20%25s%5C%20%5Ctext%7Bkg%7D%24%0A%0A%20%20%20%20%20%20%20%20Next%2C%20convert%20the%20radius%20from%20kilometres%20to%20centimetres%3A%0A%0A%20%20%20%20%20%20%20%20%24r%20%3D%20%25.2f%5C%20%5Ctext%7Bkm%7D%0A%20%20%20%20%20%20%20%20%3D%20%25.2f%20%5Ctimes%2010%5E%7B5%7D%5C%20%5Ctext%7Bcm%7D%0A%20%20%20%20%20%20%20%20%3D%20%25s%5C%20%5Ctext%7Bcm%7D%24%0A%0A%20%20%20%20%20%20%20%20Assuming%20the%20neutron%20star%20is%20spherical%2C%20its%20volume%20is%0A%0A%20%20%20%20%20%20%20%20%24V%20%3D%20%5Cfrac%7B4%7D%7B3%7D%5Cpi%20r%5E3%0A%20%20%20%20%20%20%20%20%3D%20%25s%5C%20%5Ctext%7Bcm%7D%5E3%24%0A%0A%20%20%20%20%20%20%20%20Now%20we%20can%20calculate%20its%20average%20density%3A%0A%0A%20%20%20%20%20%20%20%20%24%5Cdisplaystyle%5Crho_%7B%5Ctext%7BNS%7D%7D%20%3D%20%5Cfrac%7BM%7D%7BV%7D%0A%20%20%20%20%20%20%20%20%3D%20%5Cfrac%7B%25s%5C%20%5Ctext%7Bkg%7D%7D%7B%25s%5C%20%5Ctext%7Bcm%7D%5E3%7D%24%0A%0A%20%20%20%20%20%20%20%20Because%20**1%20cm%24%5E3%24%20%3D%201%20mL**%2C%20we%20can%20write%20this%20as%0A%0A%20%20%20%20%20%20%20%20%24%5Crho_%7B%5Ctext%7BNS%7D%7D%20%3D%20%25s%5C%20%5Ctext%7Bkg%2FmL%7D%24%0A%0A%20%20%20%20%20%20%20%20In%20other%20words%2C%20**just%201%20mL%20of%20material%20at%20this%20average%20density%0A%20%20%20%20%20%20%20%20would%20have%20a%20mass%20of%20about%20%24%25s%24%20kg!**%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20%25%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20star_solar_masses%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20star_solar_masses%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(M_SUN)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(star_mass_kg)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20star_radius_km%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20star_radius_km%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(star_radius_cm)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(star_volume_cm3)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(star_mass_kg)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(star_volume_cm3)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(density_kg_per_mL)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(density_kg_per_mL)%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20---%0A%0A%20%20%20%20%23%20Every%20human%20who%20has%20ever%20lived%0A%0A%20%20%20%20Let's%20start%20with%20the%20biggest%20version%20of%20the%20question%3A%20**what%20if%20we%0A%20%20%20%20compressed%20every%20human%20who%20has%20ever%20lived%3F**%0A%0A%20%20%20%20%23%23%20How%20many%20people%20have%20ever%20lived%3F%0A%0A%20%20%20%20The%20Population%20Reference%20Bureau%20(PRB)%20estimates%20the%20cumulative%20number%20of%0A%20%20%20%20people%20ever%20born%20at%20several%20points%20in%20time.%20We%20can%20interpolate%20between%0A%20%20%20%20those%20benchmark%20values%20to%20estimate%20the%20total%20today.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(np)%3A%0A%20%20%20%20from%20datetime%20import%20date%0A%0A%20%20%20%20%23%20Population%20Reference%20Bureau%20(PRB)%20estimates%20of%20the%20cumulative%0A%20%20%20%20%23%20number%20of%20humans%20ever%20born.%0A%20%20%20%20%23%0A%20%20%20%20%23%20SOURCE%3A%0A%20%20%20%20%23%20Kaneda%20%26%20Haub%20(2022)%2C%20Population%20Reference%20Bureau%0A%20%20%20%20%23%20https%3A%2F%2Fwww.prb.org%2Fnews%2Fhow-many-people-have-ever-lived-on-earth%2F%0A%20%20%20%20%23%0A%20%20%20%20%23%20Years%20BCE%20are%20represented%20as%20negative%20numbers%20for%20plotting.%0A%20%20%20%20%23%20PRB%20also%20gives%20a%20starting%20benchmark%20at%20190%2C000%20BCE%2C%20but%20does%20not%0A%20%20%20%20%23%20report%20a%20cumulative%20%22number%20ever%20born%22%20for%20that%20row%2C%20so%20it%20is%20not%0A%20%20%20%20%23%20included%20here.%0A%20%20%20%20all_benchmarks%20%3D%20%5B%0A%20%20%20%20%20%20%20%20(-50_000%2C%20%20%207_856_100_002)%2C%0A%20%20%20%20%20%20%20%20(-8_000%2C%20%20%20%208_993_889_771)%2C%0A%20%20%20%20%20%20%20%20(1%2C%20%20%20%20%20%20%20%2055_019_222_125)%2C%0A%20%20%20%20%20%20%20%20(1200%2C%20%20%20%20%2081_610_565_125)%2C%0A%20%20%20%20%20%20%20%20(1650%2C%20%20%20%20%2094_392_567_578)%2C%0A%20%20%20%20%20%20%20%20(1750%2C%20%20%20%20%2097_564_499_091)%2C%0A%20%20%20%20%20%20%20%20(1850%2C%20%20%20%20101_610_739_100)%2C%0A%20%20%20%20%20%20%20%20(1900%2C%20%20%20%20104_510_976_956)%2C%0A%20%20%20%20%20%20%20%20(1950%2C%20%20%20%20107_901_175_171)%2C%0A%20%20%20%20%20%20%20%20(2000%2C%20%20%20%20113_966_170_055)%2C%0A%20%20%20%20%20%20%20%20(2010%2C%20%20%20%20115_330_173_460)%2C%0A%20%20%20%20%20%20%20%20(2022%2C%20%20%20%20117_020_448_575)%2C%0A%20%20%20%20%20%20%20%20(2035%2C%20%20%20%20118_779_027_464)%2C%0A%20%20%20%20%20%20%20%20(2050%2C%20%20%20%20120_847_437_072)%2C%0A%20%20%20%20%5D%0A%0A%20%20%20%20all_years%20%3D%20np.array(%0A%20%20%20%20%20%20%20%20%5Bb%5B0%5D%20for%20b%20in%20all_benchmarks%5D%2C%0A%20%20%20%20%20%20%20%20dtype%3Dfloat%2C%0A%20%20%20%20)%0A%0A%20%20%20%20all_totals%20%3D%20np.array(%0A%20%20%20%20%20%20%20%20%5Bb%5B1%5D%20for%20b%20in%20all_benchmarks%5D%2C%0A%20%20%20%20%20%20%20%20dtype%3Dfloat%2C%0A%20%20%20%20)%0A%0A%20%20%20%20benchmarks%20%3D%20%5B%0A%20%20%20%20%20%20%20%20(2010%2C%20115_330_173_460)%2C%0A%20%20%20%20%20%20%20%20(2022%2C%20117_020_448_575)%2C%0A%20%20%20%20%20%20%20%20(2035%2C%20118_779_027_464)%2C%0A%20%20%20%20%20%20%20%20(2050%2C%20120_847_437_072)%2C%0A%20%20%20%20%5D%0A%0A%20%20%20%20years%20%3D%20np.array(%0A%20%20%20%20%20%20%20%20%5Bb%5B0%5D%20for%20b%20in%20benchmarks%5D%2C%0A%20%20%20%20%20%20%20%20dtype%3Dfloat%2C%0A%20%20%20%20)%0A%0A%20%20%20%20totals%20%3D%20np.array(%0A%20%20%20%20%20%20%20%20%5Bb%5B1%5D%20for%20b%20in%20benchmarks%5D%2C%0A%20%20%20%20%20%20%20%20dtype%3Dfloat%2C%0A%20%20%20%20)%0A%0A%20%20%20%20today%20%3D%20date(2026%2C%208%2C%2022)%0A%0A%20%20%20%20y0%20%3D%20date(today.year%2C%201%2C%201)%0A%20%20%20%20y1%20%3D%20date(today.year%20%2B%201%2C%201%2C%201)%0A%0A%20%20%20%20today_decimal%20%3D%20(%0A%20%20%20%20%20%20%20%20today.year%0A%20%20%20%20%20%20%20%20%2B%20(today%20-%20y0).days%20%2F%20(y1%20-%20y0).days%0A%20%20%20%20)%0A%0A%0A%20%20%20%20%23%20Estimate%20the%20cumulative%20number%20of%20humans%20ever%20born%20today%0A%20%20%20%20%23%20by%20straight-line%20interpolation%20between%20the%20neighbouring%0A%20%20%20%20%23%20PRB%20benchmarks.%0A%20%20%20%20humans_ever_born%20%3D%20float(%0A%20%20%20%20%20%20%20%20np.interp(today_decimal%2C%20years%2C%20totals)%0A%20%20%20%20)%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20all_totals%2C%0A%20%20%20%20%20%20%20%20all_years%2C%0A%20%20%20%20%20%20%20%20date%2C%0A%20%20%20%20%20%20%20%20humans_ever_born%2C%0A%20%20%20%20%20%20%20%20today%2C%0A%20%20%20%20%20%20%20%20today_decimal%2C%0A%20%20%20%20%20%20%20%20totals%2C%0A%20%20%20%20%20%20%20%20years%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(humans_ever_born%2C%20plt%2C%20today_decimal%2C%20totals%2C%20years)%3A%0A%20%20%20%20fig%2C%20ax%20%3D%20plt.subplots(figsize%3D(6%2C%204))%0A%0A%20%20%20%20%23%20Neutron%20Star%20Week%20branding%0A%20%20%20%20ns_cyan%20%3D%20%22%2300FFFF%22%0A%20%20%20%20ns_orange%20%3D%20%22%23FFA233%22%0A%0A%20%20%20%20ax.plot(%0A%20%20%20%20%20%20%20%20years%2C%20totals%20%2F%201e9%2C%20%22-o%22%2C%0A%20%20%20%20%20%20%20%20color%3Dns_cyan%2C%0A%20%20%20%20%20%20%20%20markerfacecolor%3Dns_cyan%2C%0A%20%20%20%20%20%20%20%20markeredgecolor%3D%22black%22%2C%0A%20%20%20%20%20%20%20%20linewidth%3D2.3%2C%0A%20%20%20%20%20%20%20%20markersize%3D6%2C%0A%20%20%20%20%20%20%20%20label%3D%22PRB%20benchmarks%22%0A%20%20%20%20)%0A%0A%20%20%20%20ax.scatter(%0A%20%20%20%20%20%20%20%20%5Btoday_decimal%5D%2C%0A%20%20%20%20%20%20%20%20%5Bhumans_ever_born%20%2F%201e9%5D%2C%0A%20%20%20%20%20%20%20%20color%3Dns_orange%2C%0A%20%20%20%20%20%20%20%20edgecolor%3D%22black%22%2C%0A%20%20%20%20%20%20%20%20linewidth%3D1%2C%0A%20%20%20%20%20%20%20%20zorder%3D6%2C%0A%20%20%20%20%20%20%20%20s%3D85%2C%0A%20%20%20%20%20%20%20%20label%3D%22Today%22%0A%20%20%20%20)%0A%0A%20%20%20%20ax.set_xlabel(%22Year%22)%0A%20%20%20%20ax.set_ylabel(%22Cumulative%20humans%20ever%20born%20(billions)%22)%0A%20%20%20%20ax.legend(%0A%20%20%20%20%20%20%20%20frameon%3DFalse%2C%0A%20%20%20%20%20%20%20%20handlelength%3D1.5%2C%0A%20%20%20%20%20%20%20%20handletextpad%3D1%2C%0A%20%20%20%20%20%20%20%20labelspacing%3D1%2C%0A%20%20%20%20%20%20%20%20borderaxespad%3D1%2C%0A%20%20%20%20%20%20%20%20scatteryoffsets%3D%5B1%5D%2C%0A%20%20%20%20)%0A%20%20%20%20ax.grid(True%2C%20alpha%3D0.2)%0A%0A%20%20%20%20fig.tight_layout()%0A%20%20%20%20fig%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(humans_ever_born%2C%20mo%2C%20today)%3A%0A%20%20%20%20mo.md(f%22%22%22%0A%20%20%20%20Interpolating%20to%20**%7Btoday%3A%25B%20%25d%2C%20%25Y%7D**%20gives%20an%20estimated%0A%20%20%20%20**%7Bhumans_ever_born%3A%2C.0f%7D%20people%20ever%20born**.%0A%0A%20%20%20%20%23%23%20How%20much%20would%20they%20weigh%3F%0A%0A%20%20%20%20To%20estimate%20the%20combined%20mass%20of%20all%20those%20people%2C%20we%20need%20to%20choose%20an%0A%20%20%20%20average%20mass%20per%20person.%20The%20same%20choice%20will%20also%20be%20used%20later%20for%20the%0A%20%20%20%20population%20alive%20today.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20%23%20Human-mass%20presets.%20%22Adult%20average%22%20~62%20kg%20is%20the%20global%20mean%20adult%20body%0A%20%20%20%20HUMAN_CUSTOM%20%3D%20%22Custom%20(drag%20or%20type)%22%0A%20%20%20%20human_presets%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22Global%20adult%20average%20(62%20kg)%22%3A%2062%2C%0A%20%20%20%20%20%20%20%20HUMAN_CUSTOM%3A%20None%2C%0A%20%20%20%20%7D%0A%20%20%20%20human_dropdown%20%3D%20mo.ui.dropdown(%0A%20%20%20%20%20%20%20%20options%3Dlist(human_presets)%2C%20value%3Dlist(human_presets)%5B0%5D%2C%0A%20%20%20%20%20%20%20%20label%3D%22Average%20human%20mass%3A%22%2C%0A%20%20%20%20)%0A%0A%20%20%20%20get_human%2C%20set_human%20%3D%20mo.state(62)%0A%20%20%20%20return%20HUMAN_CUSTOM%2C%20get_human%2C%20human_dropdown%2C%20human_presets%2C%20set_human%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(get_human%2C%20mo%2C%20set_human)%3A%0A%20%20%20%20human_slider%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%2030%2C%2080%2C%20step%3D1%2C%20value%3Dget_human()%2C%20on_change%3Dset_human%2C%0A%20%20%20%20%20%20%20%20label%3D%22Mass%20(kg)%3A%22%2C%20show_value%3DTrue%2C%0A%20%20%20%20)%0A%20%20%20%20return%20(human_slider%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(get_human%2C%20mo%2C%20set_human)%3A%0A%20%20%20%20human_box%20%3D%20mo.ui.number(%0A%20%20%20%20%20%20%20%2030%2C%2080%2C%20step%3D1%2C%20value%3Dget_human()%2C%20on_change%3Dset_human%2C%20label%3D%22or%20type%3A%22%2C%0A%20%20%20%20)%0A%20%20%20%20return%20(human_box%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(HUMAN_CUSTOM%2C%20human_box%2C%20human_dropdown%2C%20human_slider%2C%20mo)%3A%0A%20%20%20%20human_controls%20%3D%20%5Bhuman_dropdown%5D%0A%20%20%20%20if%20human_dropdown.value%20%3D%3D%20HUMAN_CUSTOM%3A%0A%20%20%20%20%20%20%20%20human_controls.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.hstack(%5Bhuman_slider%2C%20human_box%5D%2C%20justify%3D%22start%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20mo.vstack(human_controls)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20HUMAN_CUSTOM%2C%0A%20%20%20%20get_human%2C%0A%20%20%20%20human_dropdown%2C%0A%20%20%20%20human_presets%2C%0A%20%20%20%20humans_ever_born%2C%0A)%3A%0A%20%20%20%20if%20human_dropdown.value%20%3D%3D%20HUMAN_CUSTOM%3A%0A%20%20%20%20%20%20%20%20HUMAN_KG%20%3D%20get_human()%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20HUMAN_KG%20%3D%20human_presets%5Bhuman_dropdown.value%5D%0A%0A%20%20%20%20ever_human_mass_kg%20%3D%20humans_ever_born%20*%20HUMAN_KG%0A%20%20%20%20return%20HUMAN_KG%2C%20ever_human_mass_kg%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(HUMAN_KG%2C%20ever_human_mass_kg%2C%20humans_ever_born%2C%20mo%2C%20sci)%3A%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20r%22%22%22%0A%20%20%20%20%20%20%20%20Using%20an%20average%20mass%20of%20**%25.0f%20kg%20per%20person**%3A%0A%0A%20%20%20%20%20%20%20%20%24%5Cdisplaystyle%20M_%7B%5Ctext%7Ball%20humans%20ever%7D%7D%0A%20%20%20%20%20%20%20%20%3D%20N_%7B%5Ctext%7Bever%20born%7D%7D%20%5Ctimes%20%5Clangle%20M_%7B%5Ctext%7Bhuman%7D%7D%20%5Crangle%24%0A%0A%20%20%20%20%20%20%20%20%24~~~~~~~~~~~~~~~~~~~~~~%3D%20%25s%20%5Ctimes%20%25.0f%5C%20%5Ctext%7Bkg%7D%0A%20%20%20%20%20%20%20%20%3D%20%25s%5C%20%5Ctext%7Bkg%7D%24%0A%0A%20%20%20%20%20%20%20%20So%20the%20combined%20mass%20of%20every%20human%20who%20has%20ever%20lived%20is%20about%0A%20%20%20%20%20%20%20%20**%24%25s%24%20kg**.%0A%0A%20%20%20%20%20%20%20%20%23%23%20How%20much%20space%20would%20they%20take%20up%3F%0A%0A%20%20%20%20%20%20%20%20Now%20divide%20that%20total%20mass%20by%20the%20neutron%20star's%20average%20density%3A%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20%25%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20HUMAN_KG%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(humans_ever_born)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20HUMAN_KG%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(ever_human_mass_kg)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(ever_human_mass_kg)%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(density_kg_per_mL%2C%20ever_human_mass_kg)%3A%0A%20%20%20%20ever_volume_mL%20%3D%20ever_human_mass_kg%20%2F%20density_kg_per_mL%0A%20%20%20%20return%20(ever_volume_mL%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20density_kg_per_mL%2C%0A%20%20%20%20ever_human_mass_kg%2C%0A%20%20%20%20ever_volume_mL%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20sci%2C%0A%20%20%20%20star_radius_km%2C%0A%20%20%20%20star_solar_masses%2C%0A)%3A%0A%20%20%20%20volume_units%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22tsp%20(5%20mL)%22%3A%205%2C%0A%20%20%20%20%20%20%20%20%22metric%20tbsp%20(15%20mL)%22%3A%2015%2C%0A%20%20%20%20%20%20%20%20%22AUS%20tbsp%20(20%20mL)%22%3A%2020%2C%0A%20%20%20%20%20%20%20%20%22thimbles%20(~2%20mL)%22%3A%202%2C%0A%20%20%20%20%7D%0A%0A%20%20%20%20ever_summary%20%3D%20mo.md(%0A%20%20%20%20%20%20%20%20r%22%22%22%0A%20%20%20%20%20%20%20%20%24%5Cdisplaystyle%20V_%7B%5Ctext%7Ball%20humans%20ever%7D%7D%0A%20%20%20%20%20%20%20%20%3D%20%5Cfrac%7BM_%7B%5Ctext%7Ball%20humans%20ever%7D%7D%7D%7B%5Crho_%7B%5Ctext%7BNS%7D%7D%7D%0A%20%20%20%20%20%20%20%20%3D%20%5Cfrac%7B%25s%5C%20%5Ctext%7Bkg%7D%7D%7B%25s%5C%20%5Ctext%7Bkg%2FmL%7D%7D%0A%20%20%20%20%20%20%20%20%3D%20%25.2f%5C%20%5Ctext%7BmL%7D%24%0A%0A%20%20%20%20%20%20%20%20So%2C%20compressed%20to%20the%20**average%20density%20of%20this%20neutron%20star**%0A%20%20%20%20%20%20%20%20(**%25.2f%20%24M_%5Codot%24%2C%20%25.2f%20km**)%2C%20every%20human%20who%20has%20ever%20lived%0A%20%20%20%20%20%20%20%20would%20occupy%20just%20**%25.2f%20mL**.%0A%0A%20%20%20%20%20%20%20%20**In%20everyday%20terms%2C%20that's%20about%3A**%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20%25%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(ever_human_mass_kg)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(density_kg_per_mL)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ever_volume_mL%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20star_solar_masses%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20star_radius_km%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ever_volume_mL%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%0A%20%20%20%20ever_cards_html%20%3D%20%22%22.join(%0A%20%20%20%20%20%20%20%20f%22%22%22%0A%20%20%20%20%20%20%20%20%3Cdiv%20class%3D%22volume-card%22%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20%3Cdiv%20class%3D%22volume-number%22%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7Bever_volume_mL%20%2F%20size%3A.2f%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%3C%2Fdiv%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20%3Cdiv%20class%3D%22volume-label%22%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7Bname%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%3C%2Fdiv%3E%0A%20%20%20%20%20%20%20%20%3C%2Fdiv%3E%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20for%20name%2C%20size%20in%20volume_units.items()%0A%20%20%20%20)%0A%0A%20%20%20%20ever_cards%20%3D%20mo.Html(%0A%20%20%20%20%20%20%20%20f%22%22%22%0A%20%20%20%20%20%20%20%20%3Cstyle%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20.volume-grid%20%7B%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20display%3A%20grid%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20grid-template-columns%3A%20repeat(auto-fit%2C%20minmax(120px%2C%201fr))%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20gap%3A%2012px%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20width%3A%20100%25%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20margin-top%3A%208px%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%7D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20.volume-card%20%7B%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20text-align%3A%20center%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20padding%3A%2014px%2010px%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20border-top%3A%202px%20solid%20currentColor%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20border-bottom%3A%201px%20solid%20%23bbb%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%7D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20.volume-number%20%7B%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20font-size%3A%201.6rem%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20font-weight%3A%20700%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20line-height%3A%201.2%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%7D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20.volume-label%20%7B%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20margin-top%3A%206px%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20font-size%3A%200.9rem%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20line-height%3A%201.3%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%7D%0A%20%20%20%20%20%20%20%20%3C%2Fstyle%3E%0A%0A%20%20%20%20%20%20%20%20%3Cdiv%20class%3D%22volume-grid%22%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20%7Bever_cards_html%7D%0A%20%20%20%20%20%20%20%20%3C%2Fdiv%3E%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20)%0A%0A%20%20%20%20mo.vstack(%5Bever_summary%2C%20ever_cards%5D)%0A%20%20%20%20return%20(volume_units%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(ever_volume_mL%2C%20mo)%3A%0A%20%20%20%20mo.md(f%22%22%22%0A%20%20%20%20%23%23%23%20How%20has%20this%20changed%20over%20time%3F%0A%0A%20%20%20%20Today%2C%20every%20human%20who%20has%20ever%20lived%20would%20occupy%20about%0A%20%20%20%20**%7Bever_volume_mL%3A.2f%7D%20mL**%20if%20compressed%20to%20the%20average%20density%0A%20%20%20%20of%20this%20neutron%20star.%0A%0A%20%20%20%20But%20that%20volume%20has%20been%20growing%20throughout%20human%20history%20as%20more%0A%20%20%20%20people%20have%20been%20born.%0A%0A%20%20%20%20Using%20the%20Population%20Reference%20Bureau%20estimates%20at%20each%20benchmark%20year%2C%0A%20%20%20%20we%20can%20repeat%20the%20same%20calculation%3A%0A%0A%20%20%20%20**humans%20ever%20born%20%E2%86%92%20total%20human%20mass%20%E2%86%92%20equivalent%20volume%20at%20neutron-star%20density**%0A%0A%20%20%20%20The%20cyan%20points%20show%20the%20PRB%20benchmark%20estimates.%20The%20orange%20point%20marks%0A%20%20%20%20today%2C%20at%20**%7Bever_volume_mL%3A.2f%7D%20mL**.%20The%20dashed%20horizontal%20lines%20show%0A%20%20%20%20familiar%20everyday%20volumes%20for%20comparison.%0A%0A%20%20%20%20*The%20time%20axis%20uses%20logarithmic%20spacing%20so%20that%20both%20ancient%20and%20recent%0A%20%20%20%20human%20history%20can%20be%20seen%20clearly%20on%20the%20same%20plot.*%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20all_totals%2C%0A%20%20%20%20all_years%2C%0A%20%20%20%20density_kg_per_mL%2C%0A%20%20%20%20ever_human_mass_kg%2C%0A%20%20%20%20ever_volume_mL%2C%0A%20%20%20%20humans_ever_born%2C%0A%20%20%20%20np%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20today_decimal%2C%0A)%3A%0A%20%20%20%20%23%20Average%20human%20mass%20used%20in%20the%20calculation%20above%0A%20%20%20%20plot_average_human_mass_kg%20%3D%20ever_human_mass_kg%20%2F%20humans_ever_born%0A%0A%20%20%20%20%23%20At%20each%20PRB%20benchmark%20year%3A%0A%20%20%20%20%23%20cumulative%20humans%20ever%20born%20-%3E%20total%20human%20mass%20-%3E%20volume%20at%20neutron-star%20density%0A%20%20%20%20humanity_ns_volume_mL%20%3D%20(%0A%20%20%20%20%20%20%20%20all_totals%20*%20plot_average_human_mass_kg%20%2F%20density_kg_per_mL%0A%20%20%20%20)%0A%20%20%20%20%23%20Transform%20calendar%20year%20into%20%22years%20before%202051%22.%0A%20%20%20%20%23%20This%20lets%20us%20use%20a%20log%20x-axis%20while%20keeping%202050%20finite.%0A%20%20%20%20plot_years%20%3D%202051%20-%20all_years%0A%20%20%20%20plot_today%20%3D%202051%20-%20today_decimal%0A%0A%20%20%20%20volume_fig%2C%20volume_ax%20%3D%20plt.subplots(figsize%3D(6%2C%204))%0A%0A%20%20%20%20volume_ns_cyan%20%3D%20%22%2300FFFF%22%0A%20%20%20%20volume_ns_orange%20%3D%20%22%23FFA233%22%0A%0A%20%20%20%20volume_ax.plot(%0A%20%20%20%20%20%20%20%20plot_years%2C%0A%20%20%20%20%20%20%20%20humanity_ns_volume_mL%2C%0A%20%20%20%20%20%20%20%20%22-o%22%2C%0A%20%20%20%20%20%20%20%20color%3Dvolume_ns_cyan%2C%0A%20%20%20%20%20%20%20%20markerfacecolor%3Dvolume_ns_cyan%2C%0A%20%20%20%20%20%20%20%20markeredgecolor%3D%22black%22%2C%0A%20%20%20%20%20%20%20%20linewidth%3D2.0%2C%0A%20%20%20%20%20%20%20%20markersize%3D5%2C%0A%20%20%20%20%20%20%20%20label%3D%22PRB%20benchmarks%22%2C%0A%20%20%20%20)%0A%0A%20%20%20%20volume_ax.scatter(%0A%20%20%20%20%20%20%20%20%5Bplot_today%5D%2C%0A%20%20%20%20%20%20%20%20%5Bever_volume_mL%5D%2C%0A%20%20%20%20%20%20%20%20color%3Dvolume_ns_orange%2C%0A%20%20%20%20%20%20%20%20edgecolor%3D%22black%22%2C%0A%20%20%20%20%20%20%20%20linewidth%3D1%2C%0A%20%20%20%20%20%20%20%20s%3D85%2C%0A%20%20%20%20%20%20%20%20zorder%3D6%2C%0A%20%20%20%20%20%20%20%20label%3D%22Today%22%2C%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20Logarithmic%20spacing%20in%20time%0A%20%20%20%20volume_ax.set_xscale(%22log%22)%0A%20%20%20%20volume_ax.invert_xaxis()%0A%0A%20%20%20%20%23%20Everyday-volume%20reference%20lines%0A%0A%20%20%20%20volume_ax.axhline(%0A%20%20%20%20%20%20%20%205%2C%0A%20%20%20%20%20%20%20%20color%3D%220.4%22%2C%0A%20%20%20%20%20%20%20%20linestyle%3D%22--%22%2C%0A%20%20%20%20%20%20%20%20linewidth%3D1%2C%0A%20%20%20%20%20%20%20%20alpha%3D0.7%2C%0A%20%20%20%20)%0A%0A%20%20%20%20volume_ax.axhline(%0A%20%20%20%20%20%20%20%2015%2C%0A%20%20%20%20%20%20%20%20color%3D%220.4%22%2C%0A%20%20%20%20%20%20%20%20linestyle%3D%22--%22%2C%0A%20%20%20%20%20%20%20%20linewidth%3D1%2C%0A%20%20%20%20%20%20%20%20alpha%3D0.7%2C%0A%20%20%20%20)%0A%0A%20%20%20%20volume_ax.axhline(%0A%20%20%20%20%20%20%20%2020%2C%0A%20%20%20%20%20%20%20%20color%3D%220.4%22%2C%0A%20%20%20%20%20%20%20%20linestyle%3D%22--%22%2C%0A%20%20%20%20%20%20%20%20linewidth%3D1%2C%0A%20%20%20%20%20%20%20%20alpha%3D0.7%2C%0A%20%20%20%20)%0A%0A%20%20%20%20volume_ax.text(%0A%20%20%20%20%20%20%20%20plot_years.max()%2C%0A%20%20%20%20%20%20%20%205.15%2C%0A%20%20%20%20%20%20%20%20%221%20metric%20tsp%20(5%20mL)%22%2C%0A%20%20%20%20%20%20%20%20fontsize%3D8%2C%0A%20%20%20%20%20%20%20%20va%3D%22bottom%22%2C%0A%20%20%20%20)%0A%0A%0A%20%20%20%20volume_ax.text(%0A%20%20%20%20%20%20%20%20plot_years.max()%2C%0A%20%20%20%20%20%20%20%2015.15%2C%0A%20%20%20%20%20%20%20%20%221%20metric%20tbsp%20(15%20mL)%22%2C%0A%20%20%20%20%20%20%20%20fontsize%3D8%2C%0A%20%20%20%20%20%20%20%20va%3D%22bottom%22%2C%0A%20%20%20%20)%0A%0A%20%20%20%20volume_ax.text(%0A%20%20%20%20%20%20%20%20plot_years.max()%2C%0A%20%20%20%20%20%20%20%2019.85%2C%0A%20%20%20%20%20%20%20%20%221%20AUS%20tbsp%20(20%20mL)%22%2C%0A%20%20%20%20%20%20%20%20fontsize%3D8%2C%0A%20%20%20%20%20%20%20%20va%3D%22top%22%2C%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20Selected%20calendar-year%20labels%0A%20%20%20%20tick_years%20%3D%20np.array(%5B%0A%20%20%20%20%20%20%20%20-50_000%2C%0A%20%20%20%20%20%20%20%20-8_000%2C%0A%20%20%20%20%20%20%20%201%2C%0A%20%20%20%20%20%20%20%201200%2C%0A%20%20%20%20%20%20%20%201650%2C%0A%20%20%20%20%20%20%20%201900%2C%0A%20%20%20%20%20%20%20%202000%2C%0A%20%20%20%20%20%20%20%202050%2C%0A%20%20%20%20%5D)%0A%0A%20%20%20%20tick_positions%20%3D%202051%20-%20tick_years%0A%0A%20%20%20%20tick_labels%20%3D%20%5B%0A%20%20%20%20%20%20%20%20%2250%2C000%20BCE%22%2C%0A%20%20%20%20%20%20%20%20%228%2C000%20BCE%22%2C%0A%20%20%20%20%20%20%20%20%221%20CE%22%2C%0A%20%20%20%20%20%20%20%20%221200%22%2C%0A%20%20%20%20%20%20%20%20%221650%22%2C%0A%20%20%20%20%20%20%20%20%221900%22%2C%0A%20%20%20%20%20%20%20%20%222000%22%2C%0A%20%20%20%20%20%20%20%20%222050%22%2C%0A%20%20%20%20%5D%0A%0A%20%20%20%20volume_ax.set_xticks(tick_positions)%0A%20%20%20%20volume_ax.set_xticklabels(%0A%20%20%20%20%20%20%20%20tick_labels%2C%0A%20%20%20%20%20%20%20%20rotation%3D35%2C%0A%20%20%20%20%20%20%20%20ha%3D%22right%22%2C%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20Remove%20automatic%20minor%20tick%20labels%0A%20%20%20%20volume_ax.xaxis.set_minor_formatter(plt.NullFormatter())%0A%0A%20%20%20%20volume_ax.set_xlabel(%22Year%22)%0A%20%20%20%20volume_ax.set_ylabel(%0A%20%20%20%20%20%20%20%20%22Volume%20of%20all%20humans%20ever%20born%5Cnat%20neutron-star%20density%20(mL)%22%0A%20%20%20%20)%0A%0A%20%20%20%20volume_ax.legend(%0A%20%20%20%20%20%20%20%20frameon%3DFalse%2C%0A%20%20%20%20%20%20%20%20handlelength%3D1.5%2C%0A%20%20%20%20%20%20%20%20handletextpad%3D1%2C%0A%20%20%20%20%20%20%20%20labelspacing%3D1%2C%0A%20%20%20%20%20%20%20%20borderaxespad%3D1%2C%0A%20%20%20%20%20%20%20%20scatteryoffsets%3D%5B1%5D%2C%0A%20%20%20%20)%0A%0A%20%20%20%20volume_ax.grid(True%2C%20alpha%3D0.2)%0A%0A%20%20%20%20volume_fig.tight_layout()%0A%20%20%20%20volume_fig%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20---%0A%0A%20%20%20%20%23%20Every%20human%20alive%20today%0A%0A%20%20%20%20Now%20let's%20ask%20the%20same%20question%20for%20**everyone%20alive%20today**.%0A%0A%20%20%20%20%23%23%20How%20many%20people%20are%20alive%20today%3F%0A%0A%20%20%20%20For%20this%20estimate%2C%20we'll%20interpolate%20between%20the%20UN%20World%20Population%0A%20%20%20%20Prospects%202024%20population%20estimates%20for%201%20January%202026%20and%201%20January%202027.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(date%2C%20today)%3A%0A%20%20%20%20%23%20Currently-living%20population%2C%20interpolated%20from%20UN%20WPP%202024%20(medium%20variant).%0A%20%20%20%20%23%20SOURCE%3A%20UN%20DESA%2C%20World%20Population%20Prospects%202024.%0A%20%20%20%20pop_2026%20%3D%208_266_245_291%20%20%20%23%201%20Jan%202026%0A%20%20%20%20pop_2027%20%3D%208_335_111_500%20%20%20%23%201%20Jan%202027%0A%0A%20%20%20%20frac_year%20%3D%20(today%20-%20date(today.year%2C%201%2C%201)).days%20%2F%20(%0A%20%20%20%20%20%20%20%20date(today.year%20%2B%201%2C%201%2C%201)%20-%20date(today.year%2C%201%2C%201)%0A%20%20%20%20).days%0A%0A%20%20%20%20living_now%20%3D%20pop_2026%20%2B%20frac_year%20*%20(pop_2027%20-%20pop_2026)%0A%20%20%20%20return%20(living_now%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(HUMAN_KG%2C%20living_now%2C%20mo%2C%20today)%3A%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20r%22%22%22%0A%20%20%20%20%20%20%20%20Interpolating%20to%20**%25s**%20gives%20an%20estimated%0A%20%20%20%20%20%20%20%20**%25s%20people%20alive%20today**.%0A%0A%20%20%20%20%20%20%20%20%23%23%20How%20much%20would%20they%20weigh%3F%0A%0A%20%20%20%20%20%20%20%20Using%20the%20same%20average%20mass%20as%20above%2C%20**%25.0f%20kg%20per%20person**%3A%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20%25%20(today.strftime(%22%25B%20%25d%2C%20%25Y%22)%2C%20f%22%7Bliving_now%3A%2C.0f%7D%22%2C%20HUMAN_KG)%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(HUMAN_KG%2C%20living_now)%3A%0A%20%20%20%20living_human_mass_kg%20%3D%20living_now%20*%20HUMAN_KG%0A%20%20%20%20return%20(living_human_mass_kg%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(HUMAN_KG%2C%20living_human_mass_kg%2C%20living_now%2C%20mo%2C%20sci)%3A%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20r%22%22%22%0A%20%20%20%20%20%20%20%20%24M_%7B%5Ctext%7Bhumans%20alive%7D%7D%0A%20%20%20%20%20%20%20%20%3D%20N_%7B%5Ctext%7Balive%7D%7D%20%5Ctimes%20%5Clangle%20M_%7B%5Ctext%7Bhuman%7D%7D%20%5Crangle%24%0A%0A%20%20%20%20%20%20%20%20%24%3D%20%25s%20%5Ctimes%20%25.0f%5C%20%5Ctext%7Bkg%7D%0A%20%20%20%20%20%20%20%20%3D%20%25s%5C%20%5Ctext%7Bkg%7D%24%0A%0A%20%20%20%20%20%20%20%20So%20everyone%20alive%20today%20has%20a%20combined%20mass%20of%20about%20**%24%25s%24%20kg**.%0A%0A%20%20%20%20%20%20%20%20%23%23%20How%20much%20space%20would%20they%20take%20up%3F%0A%0A%20%20%20%20%20%20%20%20Again%2C%20divide%20the%20total%20mass%20by%20the%20neutron%20star's%20average%20density%3A%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20%25%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(living_now)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20HUMAN_KG%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(living_human_mass_kg)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(living_human_mass_kg)%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(density_kg_per_mL%2C%20living_human_mass_kg)%3A%0A%20%20%20%20living_volume_mL%20%3D%20living_human_mass_kg%20%2F%20density_kg_per_mL%0A%20%20%20%20return%20(living_volume_mL%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20density_kg_per_mL%2C%0A%20%20%20%20ever_volume_mL%2C%0A%20%20%20%20living_human_mass_kg%2C%0A%20%20%20%20living_volume_mL%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20sci%2C%0A%20%20%20%20star_radius_km%2C%0A%20%20%20%20star_solar_masses%2C%0A%20%20%20%20volume_units%2C%0A)%3A%0A%20%20%20%20ever_to_living_ratio%20%3D%20ever_volume_mL%20%2F%20living_volume_mL%0A%0A%20%20%20%20_living_summary%20%3D%20mo.md(%0A%20%20%20%20%20%20%20%20r%22%22%22%0A%20%20%20%20%20%20%20%20%24%5Cdisplaystyle%20V_%7B%5Ctext%7Bhumans%20alive%7D%7D%0A%20%20%20%20%20%20%20%20%3D%20%5Cfrac%7BM_%7B%5Ctext%7Bhumans%20alive%7D%7D%7D%7B%5Crho_%7B%5Crm%20NS%7D%7D%0A%20%20%20%20%20%20%20%20%3D%20%5Cfrac%7B%25s%5C%20%5Ctext%7Bkg%7D%7D%7B%25s%5C%20%5Ctext%7Bkg%2FmL%7D%7D%0A%20%20%20%20%20%20%20%20%3D%20%25.2f%5C%20%5Ctext%7BmL%7D%24%0A%0A%20%20%20%20%20%20%20%20So%2C%20compressed%20to%20the%20**average%20density%20of%20this%20neutron%20star**%0A%20%20%20%20%20%20%20%20(**%25.2f%20%24M_%5Codot%24%2C%20%25.2f%20km**)%2C%20every%20human%20alive%20today%20would%0A%20%20%20%20%20%20%20%20occupy%20just%20**%25.2f%20mL**.%0A%0A%20%20%20%20%20%20%20%20**In%20everyday%20terms%2C%20that's%20about%3A**%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20%25%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(living_human_mass_kg)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(density_kg_per_mL)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20living_volume_mL%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20star_solar_masses%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20star_radius_km%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20living_volume_mL%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%0A%20%20%20%20_living_cards_html%20%3D%20%22%22.join(%0A%20%20%20%20%20%20%20%20f%22%22%22%0A%20%20%20%20%20%20%20%20%3Cdiv%20class%3D%22living-volume-card%22%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20%3Cdiv%20class%3D%22living-volume-number%22%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7Bliving_volume_mL%20%2F%20size%3A.2f%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%3C%2Fdiv%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20%3Cdiv%20class%3D%22living-volume-label%22%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7Bname%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%3C%2Fdiv%3E%0A%20%20%20%20%20%20%20%20%3C%2Fdiv%3E%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20for%20name%2C%20size%20in%20volume_units.items()%0A%20%20%20%20)%0A%0A%20%20%20%20_living_cards%20%3D%20mo.Html(%0A%20%20%20%20%20%20%20%20f%22%22%22%0A%20%20%20%20%20%20%20%20%3Cstyle%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20.living-volume-grid%20%7B%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20display%3A%20grid%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20grid-template-columns%3A%20repeat(auto-fit%2C%20minmax(120px%2C%201fr))%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20gap%3A%2012px%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20width%3A%20100%25%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20margin-top%3A%208px%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20margin-bottom%3A%208px%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%7D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20.living-volume-card%20%7B%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20text-align%3A%20center%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20padding%3A%2014px%2010px%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20border-top%3A%202px%20solid%20currentColor%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20border-bottom%3A%201px%20solid%20%23bbb%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%7D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20.living-volume-number%20%7B%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20font-size%3A%201.6rem%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20font-weight%3A%20700%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20line-height%3A%201.2%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%7D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20.living-volume-label%20%7B%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20margin-top%3A%206px%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20font-size%3A%200.9rem%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20line-height%3A%201.3%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%7D%0A%20%20%20%20%20%20%20%20%3C%2Fstyle%3E%0A%0A%20%20%20%20%20%20%20%20%3Cdiv%20class%3D%22living-volume-grid%22%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20%7B_living_cards_html%7D%0A%20%20%20%20%20%20%20%20%3C%2Fdiv%3E%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20)%0A%0A%20%20%20%20_living_comparison%20%3D%20mo.md(%0A%20%20%20%20%20%20%20%20r%22%22%22%0A%20%20%20%20%20%20%20%20For%20comparison%2C%20**every%20human%20who%20has%20ever%20lived%20would%20occupy%20about%0A%20%20%20%20%20%20%20%20%25.1f%20times%20as%20much%20volume%20as%20everyone%20alive%20today**%20at%20the%20same%20density.%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20%25%20ever_to_living_ratio%0A%20%20%20%20)%0A%0A%20%20%20%20mo.vstack(%5B%0A%20%20%20%20%20%20%20%20_living_summary%2C%0A%20%20%20%20%20%20%20%20_living_cards%2C%0A%20%20%20%20%20%20%20%20_living_comparison%2C%0A%20%20%20%20%5D)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20---%0A%0A%20%20%20%20%23%20Bonus%3A%20what%20about%20Mount%20Everest%3F%0A%0A%20%20%20%20Another%20way%20to%20appreciate%20the%20scale%20of%20neutron-star%20densities%20is%20to%20ask%0A%20%20%20%20the%20same%20question%20about%20something%20*really*%20massive.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(density_kg_per_mL%2C%20mo%2C%20sci%2C%20star_radius_km%2C%20star_solar_masses)%3A%0A%20%20%20%20%23%20Sense%20of%20scale%3A%20Mount%20Everest%20as%20neutron-star%20matter.%0A%20%20%20%20%23%20Everest%20mass%20%3D%20measured%20volume%20x%20rock%20density.%0A%20%20%20%20%23%0A%20%20%20%20%23%20Volume%20SOURCE%3A%20Fry%20(1987)%2C%20%22Defining%20and%20Sizing-up%20Mountains%22%2C%20Summit%20magazine%0A%20%20%20%20%23%20%20%20%E2%80%94%2018.84%20cubic%20miles%2C%20from%20planimeter%20%2B%20end-area%20integration%2C%20base%20~5170%20m.%0A%20%20%20%20%23%20Density%20SOURCE%3A%20https%3A%2F%2Fmatmake.com%2Fproperties%2Fdensity-of-rocks.html%0A%20%20%20%20%23%20%20%20%E2%80%94%20Everest's%20rocks%20(limestone%2C%20marble%2C%20granite%2C%20gneiss)%20all%20fall%20in%20the%0A%20%20%20%20%23%20%20%20~2600-2800%20kg%2Fm%5E3%20range%3B%202700%20is%20a%20good%20representative%20value.%0A%0A%20%20%20%20MI_TO_KM%20%3D%201.609344%0A%20%20%20%20EVEREST_VOLUME_KM3%20%3D%2018.84%20*%20MI_TO_KM**3%0A%20%20%20%20EVEREST_RHO%20%3D%202700%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20kg%2Fm%5E3%0A%20%20%20%20EVEREST_KG%20%3D%20EVEREST_RHO%20*%20EVEREST_VOLUME_KM3%20*%201e9%20%20%20%23%20km%5E3%20-%3E%20m%5E3%20-%3E%20kg%0A%0A%20%20%20%20everest_mL%20%3D%20EVEREST_KG%20%2F%20density_kg_per_mL%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20r%22For%20a%20neutron%20star%20of%20**%25.2f%20%24M_%5Codot%24**%20and%20radius%20**%25.2f%20km**%20%22%0A%20%20%20%20%20%20%20%20r%22(average%20density%20%24%25s%24%20kg%2FmL)%2C%20all%20of%20**Mount%20Everest**%20%22%0A%20%20%20%20%20%20%20%20r%22(%24%5Csim%20%25s%24%20kg)%20compressed%20to%20that%20same%20density%20%22%0A%20%20%20%20%20%20%20%20r%22would%20occupy%20about%20**%25.0f%20mL**%20%E2%80%94%20roughly%20%25.0f%20teaspoons%20or%20%25.1f%20metric%20cups!%22%0A%20%20%20%20%20%20%20%20%25%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20star_solar_masses%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20star_radius_km%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(density_kg_per_mL)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sci(EVEREST_KG)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20everest_mL%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20everest_mL%20%2F%205%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20everest_mL%20%2F%20250%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A