PROGRAM logsteps

c     this program illustrates how to make equal log10 steps DX per
c     logarithmic decade from XMIN to XMAX; the user need only define
c     XMIN, XMAX and DX

c     the quantity dx is the inverse of the number of equal log10 steps
c     i.e., dx=0.1 is 10 steps of equal log10 in X for each decade

      implicit none

      integer       i,nbins
      double precision xmin,xmax,dx,x,pow


      xmin = 1.0d-3
      xmax = 1.0d+3
      dx   = 0.1

      nbins = nint(log10(xmax/xmin)/dx) + 1

      WRITE(6,*) xmin,xmax
      WRITE(6,*) dx,nbins

      DO 11 i=1,nbins
         pow = float(i-1)*dx
         x = xmin*(10.0d0**pow)
C         x = xmin*(float(i-1)+(10.0d0**pow))
         WRITE(6,*) i,x
 11   CONTINUE


      STOP
      END