RUNT VERSION FEB 7<r B 17ii2 0^/23/76 



SUBROUTINE CPV( X, Y, t, Cl» Sl» A9» AlO. K ) 



C ....CAUCHY PRINCIPAL VALUE INTEGHAL. 

 13 COMMON /TWO/ N,NNW, NWAVEL. ISYM, TSKIP, NC» PIE» GAMMA. M, TK.TP 



13 C0MM0N/SIX/XN(5).CN(5) 



13 REAL K 



13 IF (Y .G£. 0.0) Y • -l.OE-08 



15 TT»ATAN2( Y,X) 



24 TH«PIE/2.*TT 



C FOR NEGATIVE X,CORK£CTIQN TO RANGE OF ATAN2. 



27 IF(X.LT.O.)TH.TH+TP 



32 AA-K+Y 



34 E=EXP(AA) 



43 BB»K*X 



45 C1=C0S(B8) 



54 S1»SIN(3B) 



63 R«K*SCIRT<X**2+Y**2) 



102 SUMl-0, 



103 SUM2«0. 



104 IF(R.GE.10.)GO TG 13 

 107 SUM11«0. 



lia SUM22»0. 



Ill FAC'1.0 



113 SUMiC=l. 



114 SUM2C-1. 



115 SDLTH'O. 



116 CDLTH-O. 



117 ASSIGN 3 TO LOC 



120 IF(X.E0.O.)ASSIGN 8 TO LDC 



122 PL-1.0 



124 00 1 L'=i.lOO 



125 DL»L 



126 FAC»FAC*DL 

 130 RL-R+RL 



132 OLFAC = FAC*l)L 



133 DLTH>OL*Tri 

 135 A1=RL/0LFAC 



137 IF(ABS(CDLTH).LE.1.£-C7)G0 TO 2 



151 SUM1C«A6$(A1/SUM1) 



157 IF(SUM1C.LE.1.E-05)GQ TO 7 



165 2 COLTH»COS(OLTH) 



171 SUM11»A1*CDLTH 



172 SUMl-SUMl+SUMll 

 174 7 GO TO L0C.(3,«) 



203 9 SOM2C'0. 



204 GO TO 5 



205 3 IF(ABS(S0LTH).LE.l.F-07)Gn TO 4 

 217 SUM2C»ABS(Al/SLiM2) 



225 IF(SUM2C.LE.l.f-C»5)GD Ttl 5 



233 4 SDLTH=5IN(DLTH) 



237 SUM22«Al*SDLTri 



240 SUM2«SUM2+SUM22 



242 5 IF(SU11C.LE.i.E-i)5.AND.SUM2C.LE.l.e-t)5)'",n Tn e 



261 1 CONTINUE 



263 6 C=GAMMA+AL0G(«I*SUM1 



C DISCONTINUITY OP 2PIE If X NEGATIVE TN ■? T FIJMCTION. 



Table D-2. Continued 



131 



