program dexptest
        implicit double precision (o-z, a-h)


        f = 1.d0
        do n=1,20
           write(*,fmt='(1p2e10.3,3x,0pf9.6)') 
     &          f, dexp(f)-1.d0, (dexp(f)-1.d0)/f
           f = f / 10.d0
        end do
        
        end