double precision dd
      dd=dprod(1,2)
      print *,dd
      end