double precision d1,d2,d3
      r1=3.3
      d1=3.3
      r2=lgamma(r1)
      r3=algama(r1)
      print *,r2,r3
      d2=lgamma(d1)
      d3=dlgama(d1)
      print *,d2,d3
      x=9.87654e75
      d1=9.8765d75
      r3=lgamma(x)
      d2=dlgama(d1)
      print *,r3,d2
      stop
      end