| 28 |
C. ****************************************************************** |
C. ****************************************************************** |
| 29 |
C. |
C. |
| 30 |
IMPLICIT DOUBLE PRECISION(A-H,O-Z) |
IMPLICIT DOUBLE PRECISION(A-H,O-Z) |
| 31 |
COMMON/DELTAB/DELTA0,DELTA1 |
COMMON/DELTAB/DELTA0,DELTA1,DLT |
| 32 |
* |
* |
| 33 |
REAL VVV(3),FFF(3) |
REAL VVV(3),FFF(3) |
| 34 |
REAL*8 CHARGE, STEP, VECT(*), VOUT(*), F(4) |
REAL*8 CHARGE, STEP, VECT(*), VOUT(*), F(4) |
| 38 |
+ (XT,XYZT(1)),(YT,XYZT(2)),(ZT,XYZT(3)) |
+ (XT,XYZT(1)),(YT,XYZT(2)),(ZT,XYZT(3)) |
| 39 |
* |
* |
| 40 |
PARAMETER (MAXIT = 1992, MAXCUT = 11) |
PARAMETER (MAXIT = 1992, MAXCUT = 11) |
| 41 |
PARAMETER (EC=2.9979251D-4,DLT=1D-4,DLT32=DLT/32) |
cPP PARAMETER (EC=2.9979251D-4,DLT=1D-4,DLT32=DLT/32) |
| 42 |
PARAMETER (ZERO=0, ONE=1, TWO=2, THREE=3) |
PARAMETER (EC=2.99792458D-4) |
| 43 |
|
cPP PARAMETER (ZERO=0, ONE=1, TWO=2, THREE=3) |
| 44 |
|
PARAMETER (ZERO=0.D0, ONE=1.D0, TWO=2.D0, THREE=3.D0) |
| 45 |
PARAMETER (THIRD=ONE/THREE, HALF=ONE/TWO) |
PARAMETER (THIRD=ONE/THREE, HALF=ONE/TWO) |
| 46 |
PARAMETER (PISQUA=.986960440109D+01) |
PARAMETER (PISQUA=.986960440109D+01) |
| 47 |
PARAMETER (IX=1,IY=2,IZ=3,IPX=4,IPY=5,IPZ=6) |
PARAMETER (IX=1,IY=2,IZ=3,IPX=4,IPY=5,IPZ=6) |
| 48 |
|
|
| 49 |
REAL*8 DELTAB(3) |
REAL*8 DELTAB(3) |
| 50 |
|
REAL*8 DLT32 |
| 51 |
|
DLT32=DLT/32. |
| 52 |
|
|
| 53 |
*. |
*. |
| 54 |
*. ------------------------------------------------------------------ |
*. ------------------------------------------------------------------ |