Saturday, November 20, 2021

This Is How It's Done - Basic Code - Full Rydberg Equation Roots, part 3

Basic Program  <-- Load this basic program into the online basic interpreter at: http://www.calormen.com/jsbasic/ or use your own basic interpreter.  Will use this code to create the MS Excel macro/VBA document later.  It's a BASIC numerical iteration program that uses the sign-flip algorithm to solve for the roots of the Full Rydberg Polynomial. Erdos number TBD (IEEE references, Patents?).

Click the "Show output" button for more readable and total output log.



Full Rydberg Polynomial

For now, the results for the roots for our most recent re-derivation:
(mks units)
elementary charge, e=1.6022514387454(393)e-19  
Planck's constant, h=6.6257738721130(28)e-34  
electron mass, m_e=9.1097907336274(77)e-31  
Proton radius, r_p=8.4119972191308(31)e-16  
Rydberg constant, R_H=10973241.071596(634)  
Permittivity of free space, epsilon_0=8.8537920581144(05)e-12  
Speed of Light, c=299779058.06636(304)  
Fine-structure constant, alpha=0.0072976787382677(79)  

Digits resolution 14
Calc'd proton mass= 1.6726933172773(2)e-27
NIST   proton mass= 1.67262192369(51) x 10-27 kg
Final  fine-structure constant= 0.0072976787382677(79)
Input  fine-structure constant= 0.00729735256  (compare to above line)
Proton/electron mass ratio=1836.1490029653(637)

The Surfer, OM-IV
©2021 Mark Eric Rohrbaugh & Lyz Starwalker © 2021


3 comments:

  1. The results are slightly different from The Oracle TOPPCG2 runs due to previously only solving for 7 of the 8 roots.

    The previous attempts used m_e in the "new" term (the term that is dropped by mainstream using the reduced mass approximation). Therefore, by adding this extra degree of freedom/variable, the final convergent result is different. This is similar to holding any of the constants to a desired or known GOLDEN value, and letting the algorith optimize the others based upon choice for a specific constant.

    This can be done by setting the sign{?} = 0 instead of 1 or -1 (defaults to sign = 1).

    When all of the constants are plotted versus iteration #, as well as the sign{?} variable, one can see if things converge smoothly or if there is random chaotic behavior. I've only zoomed into one or two, so once the MS Excel Spreadsheet is complete, all can be observed.

    ReplyDelete
  2. I could always program the loop using Python which is very BASIC like, so an easy port. Gives more more motivation to progress past the "Hello World!" stage.

    ReplyDelete
  3. For Part 4 likely will simply compare the NIST/CODATA values to The Oracle TOPPCG 2.1, that'll save me wondering as well as readers. The NIST/CODATA is only accurate to 4-6 decimal places due to measurement errors and fitting to incomplete Rydberg constant equation (instead of the Full Rydberg Polynomial)

    ReplyDelete

Watch the water = Lake πŸ‘© πŸŒŠπŸ¦†