author: boulder colo joint institute for laboratory astrophysics

... RADSS * ( 2 * RHOP + RADIUS ** 2 / 3. * D2RHOP ) ) * ) N = 1 UG = 2. * RFUN ( ENS , O , ENS , 0 , 2 , SM ) / S ** 2 VG = UG + RRFCNIENS , O , ENS , 0,0 , SP ) + S ** 4 * RRFCN ( ENS , O , ENS , 0 , -4 , SP ) + 1 2. * S ** 6 * RRFCN ...