LISTING OF PROGRAM HWSP

1 PROGRAM

2 $CONSTANTS

3 PI=ACOS(-1)

4 SIGNIFIC 7

5 $

6 $DATA

7 E = 39683.163, 39760.285, 39955.053, 40991.884

8 $

9 $SLATER AND SO PARAMETERS

10 F0 = (E(4)+E(2))/2 + (E(3)-E(1))/6

11 G1 = (E(4)+E(2)-E(3)-E(1))/2

12 ZETA = (E(3)-E(1))*2/3

13 PRINT F0, G1, ZETA

14 $

15 $CHECK SPLITTING

16 EXPT = E(4)-E(2)

17 CALC = 2*SQRT((G1+ZETA/4)**2 + ZETA**2/2)

18 PRINT EXPT CALC

19 $

20 $COMPUTE MIXING ANGLE

21 COT2T=(4*G1+ZETA)/(SQRT(8)*ZETA)

22 THETA2=ACOT(COT2T)

23 THETA=(THETA2/2)*180/PI

24 PRINT THETA

25 END

:_EXECUTE HWSP

EXECUTION STARTED

F0 = 40421.4 G1 = 556.9765 ZETA = 181.26

EXPT = 1231.599 CALC = 1231.556

THETA = 6.006788