Line data    Source code

       1             : c.........................................................inconi
       2         332 :       SUBROUTINE inconi(l,xmj,e,csq,rr,xx1,xx2)
       3             : c initial point for income regular solution of dirac eq.
       4             : c order kap1=-L-1, kap2=L
       5             : c
       6             : C----->infinity point solutions construction
       7             : C     .. Scalar Arguments ..
       8             : 
       9             :       IMPLICIT NONE
      10             :       REAL csq,e,rr,xmj
      11             :       INTEGER l
      12             : C     ..
      13             : C     .. Array Arguments ..
      14             :       REAL xx1(4),xx2(4)
      15             : C     ..
      16             : C     .. Local Scalars ..
      17             :       REAL beta,bova
      18             : C     ..
      19             : C     .. Intrinsic Functions ..
      20             :       INTRINSIC abs,exp,sqrt
      21             : C     ..
      22         332 :       beta = sqrt(-e-e*e/csq)
      23         332 :       bova = -beta/ (1.+e/csq)
      24         332 :       IF (abs(xmj).LT.l) THEN
      25         170 :          xx1(1) = exp(-beta*rr)
      26         170 :          xx1(2) = bova*xx1(1)
      27         170 :          xx1(3) = 0.
      28         170 :          xx1(4) = 0.
      29             : c
      30         170 :          xx2(1) = 0.
      31         170 :          xx2(2) = 0.
      32         170 :          xx2(3) = exp(-beta*rr)
      33         170 :          xx2(4) = bova*xx2(3)
      34             :       ELSE
      35             : c      XX1(1) =1.D-05
      36         162 :          xx1(1) = exp(-beta*rr)
      37         162 :          xx1(2) = bova*xx1(1)
      38             : c not needed
      39         162 :          xx1(3) = 0.0
      40         162 :          xx1(4) = 0.0
      41         162 :          xx2(1) = 0.0
      42         162 :          xx2(2) = 0.0
      43         162 :          xx2(3) = 0.0
      44         162 :          xx2(4) = 0.0
      45             :       END IF
      46         332 :       RETURN
      47             :       END