c calculates 3-body energies from the He3 fit (Victor's 35 + 60 new points)
      implicit real*8 (a-h,o-z)
      data E2K /315774.65d0/

c      unit=1.d6  ! microhartrees
      unit=E2K   ! kelvins

      write (*,'(a)')'     R1        R2        R3         SCF'
     >// '           E1            E2            E3            E4'
     >// '            TOT           TOT+eps3'
      !read (*,*) n
      OPEN(Unit=20,File="He3Output.txt")
      n=2000
      r2=100.0
      r3=100.0
      rmin=.1
      rmax=10.0
      dr=(rmax-rmin)/(n-1)
      do 10 i=1,n
      r1=rmin+(i-1)*dr
      r2=r1
      r3=r1
      !read (*,*) R1,R2,R3
      call He3 (R1,R2,R3, scf,e1,e2,e3,eps3,e4)
      tot=scf+e1+e2+e3+e4
      write (20,'(3f10.6,7E21.12)') R1,R2,R3, scf*unit, e1*unit, e2*unit,
     >       e3*unit, e4*unit, tot*unit, (tot+eps3)*unit
c      write (*,'(3f10.6,6g14.6)') R1,R2,R3,scf
c      write (*,'(3f10.6,6g14.6)') R1,R2,R3,tot,tot+eps3
 10   continue
      end