2 March 2019; Feast of Bl. Charles the Good

 

Here is an example of the source code for getting fine structure constant in Force Fortran (double precision):

 

PROGRAM ALFA C_16

 

REAL*8 A, C, B, D, F, G, EXPM, FT, EXPP, ALFA, C_16

REAL*8 EXPM1, EXPM2, EXPM3, C_0, C_87

REAL*8 ALFA_MINUS_HALF, ALFA_SQUARE

 

C_16 = 0.999983879780488093D+01

 

C_0 = 0.986976350384356956D+00

 

C_87 = 0.314159265358979312D+01 / 0.271828182845904507D+01

 

 

open(11,file='C:/FORTRAN/Alpha C_16.txt') 

 

 

EXPM1 = (C_16/ (8.0D+00+ (2.0D+00 * (24.0D+00/24.0D+00))))**(88.0D&

&+00/24.0D+00)

 

EXPM2 = ((EXPM1)**(-1.0D+00))*((C_0)**(24.0D+00/24.0D+00))

 

EXPM3 = EXPM2 ** ((C_16) * (8.0D+00/24.0D+00))

 

EXPM = EXPM3 

 

 

write(11,*)'EXPONENT MAIN',EXPM 

 

FT = (C_0)*((C_87)**(16.0D+00 + EXPM))

 

EXPP = (16.0D+00 + (24.0D+00/24.0D+00)) / EXPM

 

ALFA_MINUS_HALF = (FT/ (8.0D+00 + (2.0D+00 * (24.0D+00 / 24.0D+00)&

&))) ** EXPP

 

ALFA_SQUARE = (ALFA_MINUS_HALF) ** 2.0D+00

 

ALFA = ( ALFA_SQUARE ) ** (-1.0D+00)

 

write(11,*)'FT VALUE AT X',FT

 

write(11,*)'EXPONENT PARTIAL',EXPP

 

write(11,*)'ALFA MINUS 1/2',ALFA_MINUS_HALF

 

write(11,*)'ALFA SQUARE',ALFA_SQUARE

 

write(11,*)'ALFA',ALFA

 

CLOSE(11)

 

STOP

 

END PROGRAM ALFA C_16

Comments powered by CComment