SUBROUTINE SOLVE ( I, FIRST )
Argument Definitions (+ indicates altered content)
INTEGER            I
LOGICAL*1          FIRST
Description
  Copyright (C) 1978 Charlotte Froese Fisher 
  Copyright (C) 2012 Jacek Kobus 
 
 
When first is .true. solve computes the potential and exchange 
function and initializes variables for the i'th radial 
equation. The vector p1 is the solution of the radial equation and 
p2 the variation of the solution with respect to the energy 
parameter e(i,i).
Source file:solve.f
I/O Operations:
Unit ID  Unit No       Access  Form   Operation
    OUC                   SEQ  FMTD           W  

Operation codes A=rewind,B=backspace,C=close,E=endfile
                I=inquire,O=open,R=read,W=write
Intrinsic Functions Called
DOUBLE PRECISION   DBLE
GENERIC*8          ABS,        SQRT
External Functions and Subroutines Called
INTEGER            MMAX,       MMIN
REAL*8             HWF,        QUAD,       ZVAR
SUBROUTINE         NMRVS,      POTL,       RELPOT,     SEARCH
SUBROUTINE         START,      XCH
Parameter Variables Used
INTEGER            MAXNO              (MAXNO = 29999)
INTEGER            MAXORB             (MAXORB = 20)
Local Variables (+ indicates altered content)
INTEGER           +J,         +JJ,         M1,         M2
INTEGER            MH
REAL*8            +A11,       +AA,        +B11,       +B3
REAL*8            +B4,        +C,         +C11,       +CD
REAL*8            +CN,        +DE1,       +DE2,       +DEL1
REAL*8             DEL2,       DELH,      +DELTA,     +DISC
REAL*8            +F1,        +FL,        +FNUM,      +HW
REAL*8             P1X,        P2X,       +PNORM,     +RL
REAL*8            +V,         +X1,        +X2,        +X3
REAL*8            +X4,        +X5,        +XP
REAL*8            +XXSOLVE(MAXNO),        +XY,        +Y1
REAL*8            +Y2,        +Y3,        +ZEROTMP(MAXNO)
Referenced Equivalenced Variables (+ indicates altered content)

EQUIV              REAL*8            +P1()       Local Var
                   REAL*8            +PDE()      COMMON /WAVE/
Referenced Common Block Variables (+ indicates altered content)
INOUT              INTEGER            OUC
LABEL              CHARACTER*4        EL(MAXORB)
MCOMMBLOCK         INTEGER            L(MAXORB)
MCOMMBLOCK         INTEGER            MAXV(MAXORB)
MCOMMBLOCK         INTEGER            N(MAXORB)
MCOMMBLOCK         REAL*8             P(MAXORB,MAXNO)
MCOMMBLOCK         REAL*8             R(MAXNO),  R2(MAXNO)
MCOMMBLOCK         REAL*8             RR(MAXNO)
MCOMMBLOCK         REAL*8            +X(MAXNO), +YK(MAXNO)
MCOMMBLOCK         REAL*8            +YR(MAXNO)
PARAM              INTEGER            ND,        NO
PARAM              REAL*8             CH,        H
PARAM0             INTEGER            IONE,      ITHREE
PARAM1             REAL*8             D0,        D1
PARAM1             REAL*8             D10,       D2
PARAM1             REAL*8             D3,        D4
PARAM1             REAL*8             D5,        D6
REL                INTEGER            IREL
REL                REAL*8             DLR(MAXORB)
REL                REAL*8             VREL(MAXNO)
TEMP2              INTEGER            KK,       +M
TEMP2              INTEGER            NJ
TEMP2              REAL*8            +DELTAE,   +EM
TEMP2              REAL*8            +EU,       +FM
TEMP2              REAL*8            +FN,       +FU
TEMP2              REAL*8            +HQ(MAXNO)
TEMP2              REAL*8            +P2(MAXNO)
TEMP2              REAL*8            +PP
WAVE               REAL*8             AZ(MAXORB)
WAVE               REAL*8            +AZD,      +E(MAXORB,MAXORB)
WAVE               REAL*8            +ED,       +PDE(MAXNO)