C
C
C
c
c
c
      subroutine funcdp(r,phi,norder,x,alpha)
      implicit real*8 (a-h,o-z)
      dimension r(3),s(3),s2(3),s3(3),s4(3),s5(3)
      dimension phi(100),s6(3)
      dimension alpha(3),x(3)
C
      fak=1.0d0
      do 1 i=1,3
      s(i)=r(i)-x(i)
      s2(i)=s(i)**2
      s3(i)=s2(i)*s(i)
      s4(i)=s2(i)*s2(i)
      s5(i)=s3(i)*s2(i)
      s6(i)=s3(i)**2
      fak=fak*(1.d0-tanh(alpha(i)*s(i)))
 1    continue
C
C
      phi(1)=1.
C
      phi(2)=s(1)+s(3)
      phi(3)=s(2)
C
      phi(4)=s2(1)+s2(3)
      phi(5)=s2(2)
      phi(6)=s(2)*(s(1)+s(3))
      phi(7)=s(1)*s(3)
C
      phi(8)=s3(1)+s3(3)
      phi(9)=s3(2)
      phi(10)=s(1)*s2(3)+s2(1)*s(3)
      phi(11)=s2(2)*(s(1)+s(3))
      phi(12)=s(2)*(s2(1)+s2(3))
      phi(13)=s(1)*s(2)*s(3)
C
      phi(14)=s4(1)+s4(3)
      phi(15)=s4(2)
      phi(16)=s3(1)*s(3)+s(1)*s3(3)
      phi(17)=s2(1)*s2(3)
      phi(18)=s(2)*(s3(1)+s3(3))
      phi(19)=s2(2)*(s2(1)+s2(3))
      phi(20)=s3(2)*(s(1)+s(3))
      phi(21)=s2(2)*s(1)*s(3)
      phi(22)=s(2)*(s2(1)*s(3)+s(1)*s2(3))
C
C
      phi(23)=s5(1)+s5(3)
      phi(24)=s4(1)*s(3)+s(1)*s4(3)
      phi(25)=s3(1)*s2(3)+s2(1)*s3(3)
      phi(26)=s(2)*(s4(1)+s4(3))
      phi(27)=s(2)*(s3(1)*s(3)+s(1)*s3(3))
      phi(28)=s(2)*s2(1)*s2(3)
      phi(29)=s2(2)*(s3(1)+s3(3))
      phi(30)=s2(2)*(s2(1)*s(3)+s(1)*s2(3))
      phi(31)=s3(2)*(s2(1)+s2(3))
      phi(32)=s3(2)*s(1)*s(3)
      phi(33)=s4(2)*(s(1)+s(3))
      phi(34)=s5(2)
C
      phi(35)=s6(1)+s6(3)
      phi(36)=s5(1)*s(3)+s(1)*s5(3)
      phi(37)=s4(1)*s2(3)+s2(1)*s4(3)
      phi(38)=s3(1)*s3(3)
      phi(39)=s(2)*(s5(1)+s5(3))
      phi(40)=s(2)*(s4(1)*s(3)+s(1)*s4(3))
      phi(41)=s(2)*(s3(1)*s2(3)+s2(1)*s3(3))
      phi(42)=s2(2)*(s4(1)+s4(3))
      phi(43)=s2(2)*(s3(1)*s(3)+s(1)*s3(3))
      phi(44)=s2(1)*s2(2)*s2(3)
      phi(45)=s3(2)*(s3(1)+s3(3))
      phi(46)=s3(2)*(s2(1)*s(3)+s(1)*s2(3))
      phi(47)=s4(2)*(s2(1)+s2(3))
      phi(48)=s4(2)*s(1)*s(3)
      phi(49)=s5(2)*(s(1)+s(3))
      phi(50)=s6(2)
C
C
C
C
C
      do 2 i=1,norder
      phi(i)=fak*phi(i)
 2    continue
C
      return
      end