c----------------------------------------------------------------------
       subroutine lyp(rho,agr,rholap,eclyp,vclyp)
c calculation of unpolarized LYP correlation energy and potential.
c
c
c author: ? 
c     in:
c    out:
c  calls: --
       implicit real*8 (a-h, o-z)
       parameter (a=0.04918d0,b=0.132d0,c=0.2533d0,d=0.349d0,
     &         cf=2.871233994d0,thd=0.333333333d0,h=2.666666667d0,
     &         h1=1.666666667d0,h2=1.3333333333d0,h3=0.666666666d0,
     &         h4=3.666666667d0)