Mercurial > octave-nkf
view libcruft/misc/d1mach.f @ 5103:e2ed74b9bfa0 after-gnuplot-split
[project @ 2004-12-28 02:43:01 by jwe]
author | jwe |
---|---|
date | Tue, 28 Dec 2004 02:43:01 +0000 |
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