complex*16 dc
      double precision dp1,dp2
      dc=(1.0d0,1.0d0)
      dp1=cdabs(dc)
      dp2=dp1*dp1
      write(6,*)dp1,dp2
      write(6,'(1h ,''cdabs test passed if 2nd number=2.0d0'')')
      stop
      end