double precision d1,d2
      x=log10(5.5)
      print *,x
      y=-4.4e3
      z=log10(y)
      stop
      end