-- CODATA Internationally recommended 2022 values -- https://physics.nist.gov/cuu/Constants/ -- c, e, h, and k are exact values a0 = 5.29177210544 10^(-11) meter alpha = 7.2973525643 10^(-3) c = 299792458.0 meter / second e = 1.602176634 10^(-19) coulomb epsilon0 = 8.8541878188 10^(-12) farad / meter h = 6.62607015 10^(-34) joule second hbar = h / float(2 pi) k = 1.380649 10^(-23) joule / kelvin me = 9.1093837139 10^(-31) kilogram mp = 1.67262192595 10^(-27) kilogram mu0 = 1.25663706127 10^(-6) newton / ampere^2 coulomb = ampere second farad = coulomb / volt joule = kilogram meter^2 / second^2 newton = kilogram meter / second^2 tesla = kilogram / second^2 / ampere volt = joule / coulomb ampere = "ampere" kelvin = "kelvin" kilogram = "kilogram" meter = "meter" second = "second" alpha e^2 / (4 float(pi) epsilon0 hbar c) err(a,b) = 2 abs((a - b) / (a + b)) -- relative error err(alpha, e^2 / (4 float(pi) epsilon0 hbar c)) Run