C FOR DP VERSION, REPLACE WITH IMPLICIT REAL*8 (A-H,O-Z)