      program DPTest

      implicit none

      double precision R, T, term1, term2, term3, term4, easy
	double precision pi, phi, C, RT, M, P, Mp
      parameter(pi = 3.14159265358979)
c     --------------------------------------------------------
c     Input
c     --------------------------------------------------------
	T = 923.15
      R = 8.314
	C =0.005
	phi = 0.5
c	--------------------------------------------------------
c	calculation
c	--------------------------------------------------------
      RT = R*T 
      term1 = 0.02*1.e-4*exp(-10115.0/T)*exp(0.5898
     #     *(1.0+2.0/pi*atan(14.985-15309.0/T)))
      term2 = 4.529e-7*exp(-(1.0/T-2.221e-4)*17767.0)
      term3 = term2/term1
      term4 = (1.0/T-2.221e-4)*26436.0
      M  = term3*exp(C*term4)
      p = phi**3*(10.0 - 15.0*phi + 6.0*phi**2)
      Mp = M**p
      write(*,*) 'M**p = ', Mp


	easy = 5.11*5.234*4.987
	write(*,*) 'easy =', easy


	write(*,*), 'atan term = ', atan(14.985-15309.0/T)
      end
 
