double precision d1,d2,d3
      d1=5.4321d73
      d2=lgamma(d1)
      print *,d2
      stop
      end