 #### Here is an example of the source code for getting the 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