-- www.eigenmath.org/planck.txt
-- To run, do select all, copy, then paste into Eigenmath's script field.
"How Planck Calculated h and k"
"Results from Kurlbaum measurement"
c = 3.0 10^10 -- speed of light
pie = float(pi)
K = 7.31 10^5 -- Kurlbaum measurement
sigma = K / (373^4 - 273^4)
P = 4 sigma / c
binding(P)
Q = P 15 c^3 / (8 pie^5)
binding(Q)
"Verify Planck numerical solution (result should be close to 1.0)"
X = 4.9651
(1 - 1/5 X) exp(X)
"Result from Lummer and Pringsheim measurement"
L = 0.294 -- Lummer and Pringsheim measurement
R = X L / c
binding(R)
"Solve for k"
k = Q R^3
k
"Solve for h"
h = R k
h
"Check h and k (1=ok)"
abs(1 - 6.548 10^(-27) / h) < 0.001
abs(1 - 1.346 10^(-16) / k) < 0.001