11 July 2019 AD; St Benedict and St Pius I

Here is a sample of the code for Transcendental Constant C_0:

PROGRAM ALFA C_M_0

 

IMPLICIT NONE

 

COMPLEX*16 F, EXPM, FT, EXPP, ALFA, J, L, R, S, Q, T, U

COMPLEX*16 EXPM1, EXPM2, EXPM3, C_M_0, C_0, C_87

COMPLEX*16 ALFA_MINUS_HALF, ALFA_SQUARE

REAL*8 X, Y, Z, POL, THETA, THETA_DEG, C_8

REAL*8 A, C, B, D, G,K, O, P

 

C_M_0 = ( 0.986976350384356956D+00, 0.0D+00 )

 

C C_0 = 0.986976350384356956D+00

 

C_0 = ( 0.986976350384356956D+00, 0.0D+00 )

 

C C_87 = 0.314159265358979312D+01 / 0.271828182845904507D+01

 

C_87 = ( 1.1557273497909217D+00, 0.0D+00 )

 

C_8 = 0.314159265358979312D+01

 

C C_8 = ( 0.314159265358979312D+01, 0.314159265358979312D+01 )

 

 

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

 

WRITE(11,*)'CONSTANT C_M_0 ', C_M_0

WRITE(11,*)'CONSTANT C_0 ', C_0

WRITE(11,*)'RATIO C_87 = C_8/C_7 ', C_87

WRITE(11,*)'CONSTANT C_8 ', C_8

WRITE(11,*)' '

WRITE(11,*)' '

 

 

 

 

B = -8.0D+00-(2.0D+00 *(24.0D+00/24.0D+00))

 

F = C_M_0 / B

 

D = -88.0D+00/24.0D+00

 

EXPM1 = F ** D

 

WRITE(11,*)'NEW EXPM1', EXPM1

 

 

 

J = (EXPM1) ** (-1.0D+00)

 

K = (-24.0D+00/24.0D+00)

 

L = (C_0) ** K

 

EXPM2 = J * L

 

WRITE(11,*)'NEW EXPM2', EXPM2

 

 

 

O = (-8.0D+00/24.0D+00)

 

P = ( C_M_0 ) * O

 

EXPM3 = EXPM2 ** P

 

WRITE(11,*)'NEW EXPM3',EXPM3

 

EXPM = EXPM3

 

 

 

R = (-0.0D+00 + EXPM)

 

S = (C_87) ** R

 

FT = (C_0) * S

 

WRITE(11,*)'NEW FT', FT

 

 

 

 

 

 

 

 

 

C EXPM1 = (C_M_0/(-8.0D+00-(2.0D+00 *(24.0D+00/24.0D+00))))**(-88.0D&

C &+00/24.0D+00)

 

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

 

C EXPM3 = EXPM2 ** ((C_M_0) * (-8.0D+00/24.0D+00))

 

C EXPM = EXPM3

 

 

 

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

 

C FT = (C_0)*((C_87)**(-0.0D+00 + EXPM))

 

 

 

Q = -0.0D+00 - (24.0D+00/24.0D+00)

 

EXPP = Q / EXPM

 

WRITE(11,*)'NEW EXPP', EXPP

 

 

 

T = -8.0D+00 - (2.0D+00 * (24.0D+00 /24.0D+00))

 

U = FT / T

 

ALFA_MINUS_HALF = U ** EXPP

 

WRITE(11,*)'NEW ALFA_MINUS_HALF',ALFA_MINUS_HALF

 

 

 

 

C EXPP = (-0.0D+00 - (24.0D+00/24.0D+00)) / EXPM

 

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

C &))) ** EXPP

 

ALFA_SQUARE = (ALFA_MINUS_HALF) ** 2.0D+00

 

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

 

X = (DREAL(ALFA))**2.0D+00

 

Y = (DIMAG(ALFA))**2.0D+00

 

Z = (X + Y)**0.5D+00

 

POL = ABS(ALFA)

 

THETA = ATAN2((DIMAG(ALFA)),DREAL(ALFA))

 

THETA_DEG = (((THETA) * (180.0D+00)) / (C_8))

 

 

 

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

 

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

 

WRITE(11,*)' '

 

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

 

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

 

WRITE(11,*)' '

 

write(11,*)'ALFA',ALFA

 

WRITE(11,*)' '

 

C write(11,*)'REAL PART OF ALFA ',X**0.5D+00

 

write(11,*)'REAL PART OF ALFA ',(DREAL(ALFA))

 

write(11,*)'IMAGINARY PART OF ALFA',(DIMAG(ALFA))

 

WRITE(11,*)' '

 

write(11,*)'POLAR PART OF ALFA',Z

 

WRITE(11,*) 'MODULUS (ABSOLUTE VALUE) OF ALFA',ABS(ALFA)

 

WRITE(11,*)' '

 

WRITE(11,*)'ANGLE THETA OF COMPLEX NUMBER (ARGUMENT) IN RADIANS', &

&THETA

 

WRITE(11,*)' '

 

WRITE(11,*)'ANGLE THETA OF COMPLEX NUMBER (ARGUMENT) IN DEGREES', &

&THETA_DEG

 

 

 

CLOSE(11)

 

 

 

 

open(22,file='C:/FORTRAN/INVIS VIS/Alpha C_M_0 VIS.txt')

 

WRITE(22, 100) ALFA

C WRITE(22,*) ' '

WRITE(22, 100) ABS(ALFA), THETA

 

CLOSE(22)

 

 

STOP

 

 

100 FORMAT (D30.20, D30.20)

 

 

 

END PROGRAM ALFA C_M_0

Comments powered by CComment