Mercurial > octave-nkf
view libcruft/misc/d1mach.f @ 7961:a5d1e27ee1f4 ss-3-1-51
3.1.51 snapshot
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Tue, 22 Jul 2008 11:40:48 -0400 |
parents | 337a09dd1c06 |
children | 32b15d5c3147 |
line wrap: on
line source
double precision function d1mach (i) integer i logical init double precision dmach(5) save init, dmach data init /.false./ if (.not. init) then call machar (dmach(1), dmach(2), dmach(3), dmach(4), dmach(5)) init = .true. endif if (i .lt. 1 .or. i .gt. 5) goto 999 d1mach = dmach(i) return 999 write(*,1999) i 1999 format(' d1mach - i out of bounds', i10) call xstopx (' ') d1mach = 0 end