x=gamma(4.56)
      print *,x
      z=gamma(1.0e-20)
      print *,z
      end