1 August 2026
7 Holy Machabees (150 BC); Saint Peter in Chains (6th Century AD); Saints Faith, Hope, Charity (2nd Century AD)
Intel FORTRAN code QUAD precision for the ‘best fit’ values – the exact value of the fine structure constant, alfa, electron magnetic moment anomaly from FERMILAB experiments 2023, and adjusted value of muon magnetic moment anomaly /slightly changed FERMILAB results/
PROGRAM PI
IMPLICIT NONE
C COMPLEX*16 RECTANGULAR, RE, IM
REAL*16 PII, EEE, SSS, FFF, SPIRITUS, FILIUS, PATER
REAL*16 ALFA, AELECTRON, AMUON, ATAU
REAL*16 PART1, PART2, PART3, PART4, CFALFA, CFALFA2
REAL*16 C_16, C_0
REAL*16 RESULT1, RESULT2
REAL*16 EXPONENT1
REAL *16 ERROREXPONENT, ERRORALFA, ERRORMAGRATIO
REAL*16 RESULTALFA, RESULTMAGRATIO
REAL*16 ELECTRONMAGMOMENT1, ELECTRONMAGMOMENT2,ELECTRONMAGMOMENTFI&
&N
REAL*16 ERRORELECTRONMAGMOMENTFIN
REAL*16 MUONMAGMOMENT1, MUONMAGMOMENT2
REAL*16 ERRORMUONMAGMOMENT, AMUONEXPERIMENT
REAL*16 FINESTRUCTURECONST
C SSS = (0.314159265358979323846D+01 / 0.271828182845904523536D+01)
C INTEGER*8 I_N
PII = 0.31415926535897932384626433832795027974790680981373Q+01
EEE = 0.27182818284590452353602874713526623143584218671935Q+01
SSS= 0.11557273497909217179100931833126962730089421974651Q+01
FFF= 0.161801828970729040507Q+01
C iNVERSE OF FINE STRUCTURE CONSTNAT ALFA
ALFA = 137.0359991817272156720641913037330100Q+00
C NEW BEST VALUE OF ALFA!
C ALFA = 0.1370359991824230700000000D+03
C NEW VALUE GUESS
C ALFA = 0.13703599919513080000000000D+03
C OFFICIAL VALUE OF ALFA
C ALFA = 0.13703599917700000000D+03
C aMOMALOUS MAGNETIC DIPOLE MOMENT OF ELECTON
AELECTRON = 0.00115965218059000000Q+00
C ANOMALOUS MAGNETIC DIPLOE MOMENT OF MUON FROM EXPERIMENT
AMUONEXPERIMENT = 0.001165920705000000000000Q+00
C ADJUSTED MUON ANOMALOUS MAGNETIC DIPOLE MOMENT
AMUON = 0.11659207061348750353739127905927Q-02
C aNOMALOUS MAGNETIC DIPOLE MOMENT OF TAU PARTICLE /FROM THEORY/
ATAU = 0.0011772100000000000000000Q+00
C EXACT VALUE OF CONSTANT C_16 AND C_0
C_16 = 0.99998387978048801818832774568126350047937555873307Q+01
C_0 = 0.98697635038435719239551219363426327409509390153844Q+00
CCCCCC FFF = FIBONACCI = (7/5)X(PI/E)
SPIRITUS = 0.26179831858273051437Q+01
CCCCCC SPIRITUS IS FIBONACCI SQUARED = [(7/5)X(PI/E)] SQUARED
FILIUS = 0.52778972135218930195Q+01
CCCCCC FILIUS = (5/7)X(E SQUARED)
PATER = 0.13817446161525102066Q+02
CCCCCC PATER = (7/5)X(PI SQUARED)
C CALCULATIONS
PART1 = ALFA / ( PII ** 2.0Q+00 )
PART2 = 100.0Q+00 / (( C_16 / 10.0Q+00 ) ** ( 3.0Q+00 / ALFA ) )
C PART2 = 100.0D+00 / (( C_16 / 10.0D+00 ) ** ( 3.0D+00 / ALFA ))
PART3 = (((( AMUON / AELECTRON ) * ( ( C_16 / 10.0Q+00 ) ** ( 1.0&
&0Q+00 / 3.0Q+00 ) ) - 1.0Q+00 ) / 10.0Q+00 ) + 1.0Q+00 )
PART4 = ( ( PII / EEE ) ** ( 50.0Q+00 ) )
RESULT1 = PART1 * PART2 * PART3
RESULT2 = PART4
EXPONENT1 = ( LOG ( RESULT1 ) ) / ( LOG ( PII / EEE ) )
ERROREXPONENT = ( ( EXPONENT1 - 50.0Q+00 ) / 50.0Q+00 )
C CALCULATION OF ALFA 137
RESULTALFA = RESULT2 * ( PII ** 2.0Q+00 ) / ( PART2 * PART3 )
ERRORALFA = ( ALFA - RESULTALFA ) / ALFA
C CALCULATION RATIO OF MAGNETIC MOMENT
RESULTMAGRATIO = RESULT2 / ( PART1 * PART2 * ( PART3 / ( AMUON /&
&AELECTRON )))
ERRORMAGRATIO = ( ( AMUONEXPERIMENT / AELECTRON ) - RESULTMAGRATIO&
& ) / ( AMUON / AELECTRON )
C ELECTRON MAGNETIC MOMENT
ELECTRONMAGMOMENT1 = RESULT2 / ( PART1 * PART2 )
ELECTRONMAGMOMENT2 = ((( ELECTRONMAGMOMENT1 - 1.0Q+00 ) * 10.0Q+00&
&) + 1.0Q+00 ) / ( ( ( C_16 / 10.0Q+00 ) ** ( 1.00Q+00 / 3.0Q+00 &
&) * AMUON ) )
ELECTRONMAGMOMENTFIN = ( ELECTRONMAGMOMENT2 ) ** ( -1.0Q+00 )
ERRORELECTRONMAGMOMENTFIN = ( AELECTRON - ELECTRONMAGMOMENTFIN ) /&
&AELECTRON
C CALCULATE MUON MAGNETIC MOMENT
MUONMAGMOMENT1 = RESULT2 / ( PART1 * PART2 )
MUONMAGMOMENT2 = ((( MUONMAGMOMENT1 - 1.0Q+00 ) * 10.0Q+00 ) + 1.0&
&Q+00 ) / ( ( ( C_16 / 10.0Q+00 ) ** ( 1.00Q+00 / 3.0Q+00 ))) * AE&
&LECTRON
C MUONMAGMOMENTFIN = ( MUONMAGMOMENT2 ) ** ( -1.0D+00 )
ERRORMUONMAGMOMENT = ( MUONMAGMOMENT2 - AMUONEXPERIMENT ) / AMUON
CFALFA = RESULT2 / ( 100.0Q+00 * PART1 )
CFALFA2 = ( PART2 / 100.00Q+00 ) * PART3
OPEN(11,FILE='C:/FORTRAN/ALFA 181727 ALFA AND MUON ADJUSTED QUAD 7&
&33 CF ALFA 2023.TXT')
WRITE(11,*),'...EXACT THEORETICAL VALUE OF ALFA 33 SIG DIGITS....'
WRITE(11,*)'...ALFA.....................',ALFA
WRITE(11,150), ALFA
WRITE(11,*),'................................................... '
WRITE(11,*),'...EXACT EXPERIMENTAL VALUE OF ELECTRON MAG MOMENT..'
WRITE(11,*)'...AELECTRON................',AELECTRON
WRITE(11, 150), AELECTRON
WRITE(11,*),'................................................... '
WRITE(11,*),'...EXACT ADJUSTED VALUE OF MUON MAG MOMENT..........'
WRITE(11,*)'...AMUON....................',AMUON
WRITE(11, 150), AMUON
WRITE(11,*),'................................................... '
WRITE(11,*),'...RESULT1 = LEFT AHND SIDE OF THE EQUATION.........'
WRITE(11,*)'...RESULT1..................',RESULT1
WRITE(11, 150), RESULT1
WRITE(11,*),'................................................... '
WRITE(11,*)'.....................................................'
WRITE(11,*)'...EXPONENT1 = EXPONET OF (PI/e) ....................'
WRITE(11,*),'EXPONENT1..................',EXPONENT1
WRITE(11, 150), EXPONENT1
WRITE(11,*)'.....................................................'
WRITE(11,*),'...RELATIVE ERROR OF EXPONENT.......................'
WRITE(11,*),'...RELATIVE ERROR OF EXPONENT..',ERROREXPONENT
WRITE(11, 150), ERROREXPONENT
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*),'...RESULT OF ALFA FROM CALCULATIONS.................'
WRITE(11,*),'.RESULT OF ALFA FROM CALCULATIONSALFA.',RESULTALFA
WRITE(11, 150), RESULTALFA
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*),'...RELATIVE ERROR OF ALFA...........................'
WRITE(11,*),'...RELATIVE ERROR OF ALFA......',ERRORALFA
WRITE(11, 150), ERRORALFA
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*),'.RESULT MAGNETIC RATIO OF MUON/ELECTRON.............'
WRITE(11,*),'.RESULT MAG RATIO OF MUON/ELECTRON..',RESULTMAGRATIO
WRITE(11, 150), RESULTMAGRATIO
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*),'...RESULT MAGNETIC RATIO FROME XPERIMENT............'
WRITE(11,*),'. RESULT MAGNETIC RATIO FROM EXPERIMENT.',(AMUONEXPE&
&RIMENT / AELECTRON)
WRITE(11, 150), ( AMUONEXPERIMENT / AELECTRON )
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*),'...RELATIVE ERROR OF MAGNETIC RATIO MUON/ELECTRON...'
WRITE(11,*),'...RELATIVE ERROR OF MAGNETIC RATIO...',ERRORMAGRATIO
WRITE(11, 150), ERRORMAGRATIO
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*),'...ELECTRONMAGMOMENT1 /SECONDARY CALCULATION/.......'
WRITE(11,*),'...ELECTRONMAGMOMENT1 /SECONDARY CALCULATION/.',ELECT&
&RONMAGMOMENT1
WRITE(11, 150), ELECTRONMAGMOMENT1
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*),'...ELECTRONMAGMOMENT2 /SECONDARY CALCULATION/.......'
WRITE(11,*),'...ELECTRONMAGMOMENT2 /SECONDARY CALCULATION/.',ELECT&
&RONMAGMOMENT2
WRITE(11, 150), ELECTRONMAGMOMENT2
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*),'...ELECTRON MAGNETIC MOMENT FINAL FROM CALCULATIONS.'
WRITE(11,*),'...ELECTRON MAGNETIC MOMENT FINAL FROM CALCULATIONS.'&
&,ELECTRONMAGMOMENTFIN
WRITE(11, 150), ELECTRONMAGMOMENTFIN
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*),'...RELATIVE ERROR OF ELECTRON MAGNETIC MOMENT = FINAL&
& AND EXPERIMENTAL..'
WRITE(11,*),'...RELATIVE ERROR OF ELECTRON MAGNETIC MOMENT....',ER&
&RORELECTRONMAGMOMENTFIN
WRITE(11, 150), ERRORELECTRONMAGMOMENTFIN
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*),'...MUON MAGNETIC MOMENT1 /SECONDARY CALCULATIONS/...'
WRITE(11,*),'...MUON MAGNETIC MOMENT1 /SECONDARY CALCULATIONS/....&
&...',MUONMAGMOMENT1
WRITE(11, 150), MUONMAGMOMENT1
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*),'...MUON MAGNETIC MOMENT FINAL FROM CALCULATIONS.....'
WRITE(11,*),'...MUON MAGNETIC MOMENT FINAL FROM CALCULATIONS.....'&
&,MUONMAGMOMENT2
WRITE(11, 150), MUONMAGMOMENT2
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*),'...RELATIVE ERROR OF MUON MAGNETIC MOMENT = CALCULATE&
&D AND EXPERIMENTAL........'
WRITE(11,*),'...RELATIVE ERROR OF MUON MAGNETIC MOMENT = CALCULATE&
&D AND EXPERIMENTAL......',ERRORMUONMAGMOMENT
WRITE(11, 150), ERRORMUONMAGMOMENT
FINESTRUCTURECONST = ( 1.0Q+00 / ALFA )
WRITE(11,*)'.....................................................'
WRITE(11,*),'...FINE STRUCTURE CONSTANT 33 SIG DIGITS............'
WRITE(11,*),'...FINE STRUCTURE CONSTANT..',FINESTRUCTURECONST
WRITE(11, 150), FINESTRUCTURECONST
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*),'...CORRECTION FACTOR ALPHA 33 SIG DIGITS............'
WRITE(11,*),'...CORRECTION FACTOR ALPHA..',CFALFA
WRITE(11, 150), CFALFA
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*),'...CORRECTION FACTOR ALPHA 33 SIG DIGITS............'
WRITE(11,*),'...CORRECTION FACTOR ALPHA..',CFALFA2
WRITE(11, 150), CFALFA2
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
WRITE(11,*)'.....................................................'
C WRITE(22, 100) ALFA
C WRITE(22, 100) ABS(ALFA), THETA
CLOSE(11)
C CLOSE(22)
C PAUSE
STOP
100 format(D30.15,D30.15)
150 format(E65.40)
C 100 format(I5,D30.15,D30.15)
200 format(D100.15)
END PROGRAM PI
Comments powered by CComment